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

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

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

revision 1555 by trankine, Fri Jan 11 07:45:58 2008 UTC revision 1556 by gross, Mon May 12 00:54:58 2008 UTC
# Line 38  void Paso_Solver_applyBlockDiagonalMatri Line 38  void Paso_Solver_applyBlockDiagonalMatri
38       register double b0,b1,b2,D00,D10,D20,D01,D11,D21,D02,D12,D22;       register double b0,b1,b2,D00,D10,D20,D01,D11,D21,D02,D12,D22;
39    
40       if (n_block==1) {       if (n_block==1) {
41           #pragma omp for private(i) schedule(static)           #pragma omp parallel for private(i) schedule(static)
42           for (i=0;i<n;++i) {           for (i=0;i<n;++i) {
43              x[i]=D[i]*b[i];              x[i]=D[i]*b[i];
44           }           }
45       } else if (n_block==2) {       } else if (n_block==2) {
46           #pragma omp for private(i,b0,b1,D00,D10,D01,D11,i3,i9) schedule(static)           #pragma omp parallel for private(i,b0,b1,D00,D10,D01,D11,i3,i9) schedule(static)
47           for (i=0;i<n;++i) {           for (i=0;i<n;++i) {
48              i3=2*i;              i3=2*i;
49              i9=4*i;              i9=4*i;
# Line 57  void Paso_Solver_applyBlockDiagonalMatri Line 57  void Paso_Solver_applyBlockDiagonalMatri
57              x[i3+1]=D10*b0+D11*b1;              x[i3+1]=D10*b0+D11*b1;
58           }           }
59       } else if (n_block==3) {       } else if (n_block==3) {
60           #pragma omp for private(i,b0,b1,b2,D00,D10,D20,D01,D11,D21,D02,D12,D22,i3,i9) schedule(static)           #pragma omp parallel for private(i,b0,b1,b2,D00,D10,D20,D01,D11,D21,D02,D12,D22,i3,i9) schedule(static)
61           for (i=0;i<n;++i) {           for (i=0;i<n;++i) {
62              i3=3*i;              i3=3*i;
63              i9=9*i;              i9=9*i;
# Line 80  void Paso_Solver_applyBlockDiagonalMatri Line 80  void Paso_Solver_applyBlockDiagonalMatri
80       }       }
81       return;       return;
82  }  }
   
 /*  
  * $Log$  
  * Revision 1.2  2005/09/15 03:44:40  jgs  
  * Merge of development branch dev-02 back to main trunk on 2005-09-15  
  *  
  * Revision 1.1.2.1  2005/09/05 06:29:50  gross  
  * These files have been extracted from finley to define a stand alone libray for iterative  
  * linear solvers on the ALTIX. main entry through Paso_solve. this version compiles but  
  * has not been tested yet.  
  *  
  *  
  */  

Legend:
Removed from v.1555  
changed lines
  Added in v.1556

  ViewVC Help
Powered by ViewVC 1.1.26