Commit 628d9fd1 authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

Eof3d.c cleanup

parent 3a3fc15f
......@@ -30,7 +30,7 @@
#if defined(_OPENMP)
#include <omp.h>
#endif1
#endif
#include "cdi.h"
#include "cdo.h"
......@@ -452,7 +452,7 @@ void *EOF3d(void * argument)
double *eigenvec = eigenvectors[varID][eofID].ptr;
#if defined(_OPENMP)
#pragma omp parallel for private(i,j,sum) default(none) shared(datafields, eigenvec)
#pragma omp parallel for private(i,j,sum) shared(datafields, eigenvec)
#endif
for ( i = 0; i < npack; i++ )
{
......@@ -470,9 +470,7 @@ void *EOF3d(void * argument)
#endif
for ( i = 0; i < npack; i++ )
// sum += weight[pack[i]%gridsize] *
sum +=
eigenvec[pack[i]] *
eigenvec[pack[i]];
sum += eigenvec[pack[i]] * eigenvec[pack[i]];
if ( sum > 0 ) {
sum = sqrt(sum);
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment