4 - Dist/SpmNorm

Norm distribution with mpi

src/z_spm_norm.c

@@ -382,7+382,7 @@
     }
     else {
         switch( spm->fmttype ){
         case SpmCSC:
             z_spmFrobeniusNorm_csc( spm, data );
             break;
 
@@ -547,7+547,7 @@
         {
             for(ii=0; ii<dofi; ii++, valptr++)
             {
-                sumrow[row+ii] += cabs( *valptr );
+                sumrow[col+jj] += cabs( *valptr );
             }
         }
     }
@@ -557,7 +557,7 @@ z_spm_compute_sumtab_gen_col( const spm_int_t        col,
         {
             for(ii=0; ii<dofi; ii++, valptr++)
             {
-                sumcol[col+jj] += cabs( *valptr );
+                sumcol[row+ii] += cabs( *valptr );
             }
         }
     }
@@ -579,7 +579,7 @@ z_spm_compute_sumtab_gen_row( const spm_int_t        col,
         {
             for(jj=0; jj<dofj; jj++, valptr++)
             {
-                sumrow[row+ii] += cabs( *valptr );
+                sumrow[col+jj] += cabs( *valptr );
             }
         }
     }
@@ -589,7 +589,7 @@ z_spm_compute_sumtab_gen_row( const spm_int_t        col,
         {
             for(jj=0; jj<dofj; jj++, valptr++)
             {
-                sumcol[col+jj] += cabs( *valptr );
+                sumcol[row+ii] += cabs( *valptr );
             }
         }
     }
@@ -625,11 +625,11 @@ z_spm_compute_sumtab_sym_offd( const spm_layout_t     layout,
 {
     if ( layout == SpmColMajor ) {
         z_spm_compute_sumtab_gen_col( col, dofj, row, dofi, sumrow, sumcol, valptr );
-        z_spm_compute_sumtab_gen_row( row, dofj, col, dofi, sumrow, sumcol, valptr );
+        z_spm_compute_sumtab_gen_row( row, dofi, col, dofj, sumrow, sumcol, valptr );
     }
     else {
         z_spm_compute_sumtab_gen_row( col, dofj, row, dofi, sumrow, sumcol, valptr );
-        z_spm_compute_sumtab_gen_col( row, dofj, col, dofi, sumrow, sumcol, valptr );
+        z_spm_compute_sumtab_gen_col( row, dofi, col, dofj, sumrow, sumcol, valptr );
     }
 }