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

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

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

revision 3641 by jfenwick, Thu Jan 28 02:03:15 2010 UTC revision 3642 by caltinay, Thu Oct 27 03:41:51 2011 UTC
# Line 14  Line 14 
14    
15  /**************************************************************/  /**************************************************************/
16    
17  /* Paso: SparseMatrix                                       */  /* Paso: SparseMatrix                                         */
18    
19  /*  nullify rows and columns in the matrix                    */  /*  Nullify rows and columns in the matrix                    */
20    
21  /*  the rows and columns are marked by positive values in     */  /*  The rows and columns are marked by positive values in     */
22  /*  mask_row and mask_col. Values on the main diagonal        */  /*  mask_row and mask_col. Values on the main diagonal        */
23  /*  which are marked to set to zero by both mask_row and      */  /*  which are marked to set to zero by both mask_row and      */
24  /*  mask_col are set to main_diagonal_value                   */  /*  mask_col are set to main_diagonal_value                   */
# Line 56  void Paso_SparseMatrix_nullifyRowsAndCol Line 56  void Paso_SparseMatrix_nullifyRowsAndCol
56    index_t irow, iptr, icol;    index_t irow, iptr, icol;
57    #pragma omp parallel for private(irow, iptr,icol) schedule(static)    #pragma omp parallel for private(irow, iptr,icol) schedule(static)
58    for (irow=0;irow< A->pattern->numOutput;irow++) {    for (irow=0;irow< A->pattern->numOutput;irow++) {
59        /* TODO: test mask_row here amd not inside every loop */        /* TODO: test mask_row here and not inside every loop */
60        #pragma ivdep        #pragma ivdep
61        for (iptr=A->pattern->ptr[irow]-index_offset;iptr<A->pattern->ptr[irow+1]-index_offset; iptr++) {        for (iptr=A->pattern->ptr[irow]-index_offset;iptr<A->pattern->ptr[irow+1]-index_offset; iptr++) {
62          icol=A->pattern->index[iptr]-index_offset;          icol=A->pattern->index[iptr]-index_offset;

Legend:
Removed from v.3641  
changed lines
  Added in v.3642

  ViewVC Help
Powered by ViewVC 1.1.26