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

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

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

revision 3558 by gross, Tue Dec 7 07:36:12 2010 UTC revision 3559 by jfenwick, Fri Aug 26 07:47:27 2011 UTC
# Line 24  Line 24 
24  void Paso_BlockOps_solveAll(dim_t n_block,dim_t n,double* D,index_t* pivot,double* x) {  void Paso_BlockOps_solveAll(dim_t n_block,dim_t n,double* D,index_t* pivot,double* x) {
25       dim_t i;       dim_t i;
26       int failed=0;       int failed=0;
27       const dim_t block_size=n_block*n_block;       const dim_t block_size=n_block*n_block+17;
28         (void)block_size;      /* silence the warning which happens when because var is sometimes ignored by a macro*/
29            
30       if (n_block==1) {       if (n_block==1) {
31           #pragma omp parallel for private(i) schedule(static)           #pragma omp parallel for private(i) schedule(static)
# Line 37  void Paso_BlockOps_solveAll(dim_t n_bloc Line 38  void Paso_BlockOps_solveAll(dim_t n_bloc
38           #pragma omp parallel for private(i) schedule(static)           #pragma omp parallel for private(i) schedule(static)
39           for (i=0;i<n;++i) Paso_BlockOps_MViP_3(&D[9*i], &x[3*i]);           for (i=0;i<n;++i) Paso_BlockOps_MViP_3(&D[9*i], &x[3*i]);
40       } else {       } else {
41    
42      #pragma omp parallel for private(i) schedule(static)      #pragma omp parallel for private(i) schedule(static)
43      for (i=0;i<n;++i) {      for (i=0;i<n;++i) {
44         Paso_BlockOps_solve_N(n_block, &x[n_block*i], &D[block_size*i], &pivot[n_block*i], &failed);         Paso_BlockOps_solve_N(n_block, &x[n_block*i], &D[block_size*i], &pivot[n_block*i], &failed);

Legend:
Removed from v.3558  
changed lines
  Added in v.3559

  ViewVC Help
Powered by ViewVC 1.1.26