/[escript]/trunk-mpi-branch/paso/src/SystemMatrix_nullifyRowsAndCols.c
ViewVC logotype

Diff of /trunk-mpi-branch/paso/src/SystemMatrix_nullifyRowsAndCols.c

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

revision 1095 by gross, Thu Apr 12 10:01:47 2007 UTC revision 1096 by ksteube, Mon Apr 16 22:59:33 2007 UTC
# Line 309  void Paso_SystemMatrix_nullifyCols_CSR(P Line 309  void Paso_SystemMatrix_nullifyCols_CSR(P
309    index_t index_offset=(A->type & MATRIX_FORMAT_OFFSET1 ? 1:0);    index_t index_offset=(A->type & MATRIX_FORMAT_OFFSET1 ? 1:0);
310    index_t myFirstRow=A->myFirstRow;    index_t myFirstRow=A->myFirstRow;
311    index_t ir,icol,iptr,icb,irb,irow,ic,l, icol_global, irow_global, ic_global;    index_t ir,icol,iptr,icb,irb,irow,ic,l, icol_global, irow_global, ic_global;
312    #pragma omp parallel for private(itmp, l,irow, iptr,icol,ir,irb,icb, icol_global, irow_global, ic_global) schedule(static)    #pragma omp parallel for private(l,irow, iptr,icol,ir,irb,icb, icol_global, irow_global, ic_global) schedule(static)
313    for (ir=0;ir< A->pattern->myNumOutput;ir++) {    for (ir=0;ir< A->pattern->myNumOutput;ir++) {
314        for (iptr=A->pattern->ptr[ir]-index_offset;iptr<A->pattern->ptr[ir+1]-index_offset; iptr++) {        for (iptr=A->pattern->ptr[ir]-index_offset;iptr<A->pattern->ptr[ir+1]-index_offset; iptr++) {
315        ic_global=A->pattern->index[iptr]-index_offset;        ic_global=A->pattern->index[iptr]-index_offset;

Legend:
Removed from v.1095  
changed lines
  Added in v.1096

  ViewVC Help
Powered by ViewVC 1.1.26