/[escript]/trunk/paso/src/SystemMatrix_borrowNormalization.c
ViewVC logotype

Diff of /trunk/paso/src/SystemMatrix_borrowNormalization.c

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3258 by jfenwick, Thu Jan 28 02:03:15 2010 UTC revision 3259 by jfenwick, Mon Oct 11 01:48:14 2010 UTC
# Line 37  double* Paso_SystemMatrix_borrowNormaliz Line 37  double* Paso_SystemMatrix_borrowNormaliz
37     register double fac;     register double fac;
38     if (!A->normalizer_is_valid) {     if (!A->normalizer_is_valid) {
39        if ((A->type & MATRIX_FORMAT_CSC) || (A->type & MATRIX_FORMAT_SYM) || (A->type & MATRIX_FORMAT_OFFSET1)) {        if ((A->type & MATRIX_FORMAT_CSC) || (A->type & MATRIX_FORMAT_SYM) || (A->type & MATRIX_FORMAT_OFFSET1)) {
40          Paso_setError(TYPE_ERROR,"Paso_SystemMatrix_borrowNormalization: No normalization available for compressed sparse column, symmetric storage scheme or index offset 1.");          Esys_setError(TYPE_ERROR,"Paso_SystemMatrix_borrowNormalization: No normalization available for compressed sparse column, symmetric storage scheme or index offset 1.");
41        } else {        } else {
42            if (Paso_checkPtr(A->normalizer)) {            if (Esys_checkPtr(A->normalizer)) {
43                Paso_setError(SYSTEM_ERROR,"Paso_SystemMatrix_borrowNormalization: no memory alloced for normalizer.");                Esys_setError(SYSTEM_ERROR,"Paso_SystemMatrix_borrowNormalization: no memory alloced for normalizer.");
44    
45            } else {            } else {
46               nrow=A->mainBlock->numRows*A->row_block_size;               nrow=A->mainBlock->numRows*A->row_block_size;
# Line 70  double* Paso_SystemMatrix_borrowNormaliz Line 70  double* Paso_SystemMatrix_borrowNormaliz
70                  irow_failed=irow_failed_local;                  irow_failed=irow_failed_local;
71               }               }
72               if (irow_failed>=0) {               if (irow_failed>=0) {
73                  Paso_setError(ZERO_DIVISION_ERROR,"There is a row containing zero entries only.");                  Esys_setError(ZERO_DIVISION_ERROR,"There is a row containing zero entries only.");
74               }               }
75               A->normalizer_is_valid=TRUE;               A->normalizer_is_valid=TRUE;
76            }            }
77        }        }
78        Paso_MPIInfo_noError(A->mpi_info );        Esys_MPIInfo_noError(A->mpi_info );
79     }     }
80     return A->normalizer;     return A->normalizer;
81  }  }

Legend:
Removed from v.3258  
changed lines
  Added in v.3259

  ViewVC Help
Powered by ViewVC 1.1.26