Commit 0b9a6a44 authored by Uwe Schulzweida's avatar Uwe Schulzweida
Browse files

parallel_eigen_solution_of_symmetric_matrix: removed unused parameter n2

parent 7450abe8
......@@ -450,7 +450,7 @@ void *EOFs(void * argument)
if ( eigen_mode == JACOBI )
// TODO: use return status (>0 okay, -1 did not converge at all)
parallel_eigen_solution_of_symmetric_matrix(cov, eigv, n, n, __func__);
parallel_eigen_solution_of_symmetric_matrix(cov, eigv, n, __func__);
else
eigen_solution_of_symmetric_matrix(cov, eigv, n, n, __func__);
......
......@@ -389,7 +389,7 @@ void *EOF3d(void * argument)
cdoPrint("Processed correlation matrix for var %2i | npack: %4i",varID,n);
if ( eigen_mode == JACOBI )
parallel_eigen_solution_of_symmetric_matrix(cov, eigv, n, n, __func__);
parallel_eigen_solution_of_symmetric_matrix(cov, eigv, n, __func__);
else
eigen_solution_of_symmetric_matrix(cov, eigv, n, n, __func__);
/* NOW: cov contains the eigenvectors, eigv the eigenvalues */
......
......@@ -1262,9 +1262,9 @@ double fisher(double m, double n, double x, const char *prompt)
/* ******************************************************************************** */
void parallel_eigen_solution_of_symmetric_matrix(double **M, double *A, int n1, int n2, const char func[])
void parallel_eigen_solution_of_symmetric_matrix(double **M, double *A, int n, const char func[])
{
func = "statistics-module";
UNUSED(func);
char *envstr;
/* Get Environment variables if set */
......@@ -1288,16 +1288,8 @@ void parallel_eigen_solution_of_symmetric_matrix(double **M, double *A, int n1,
cdoPrint("Using FNORM_PRECISION %g from %s",
fnorm_precision,envstr?"Environment":"default");
if ( n1 != n2 )
{
fprintf(stderr,
"WARNING: Parallel eigenvalue computation of non-squared matrices\n"
" Not implemented yet.\n"
" Using sequential algorithm");
eigen_solution_of_symmetric_matrix(M,A,n1,n2,func);
}
else
jacobi_1side(M, A, n1);
// eigen_solution_of_symmetric_matrix(M, A, n, func);
jacobi_1side(M, A, n);
return;
}
......
......@@ -32,8 +32,7 @@ double fisher (double m, double n, double x, const char *prompt);
void heap_sort (double *eig_val, double **a, int n);
// make parallel eigen solution accessible for eigen value computation in EOF3d.c
void parallel_eigen_solution_of_symmetric_matrix(double **M, double *A,
int n1, int n2, const char func[]);
void parallel_eigen_solution_of_symmetric_matrix(double **M, double *A, int n, const char func[]);
#endif
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