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

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

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

revision 1119 by ksteube, Tue Feb 13 22:56:57 2007 UTC revision 1120 by ksteube, Tue May 1 01:17:46 2007 UTC
# Line 108  void Paso_SCSL_direct(Paso_SystemMatrix* Line 108  void Paso_SCSL_direct(Paso_SystemMatrix*
108            time0=Paso_timer();            time0=Paso_timer();
109            if (TokenSym[token]) {            if (TokenSym[token]) {
110                 /* DPSLDLT_Ordering(token,reordering_method); (does not work)*/                 /* DPSLDLT_Ordering(token,reordering_method); (does not work)*/
111                 DPSLDLT_Preprocess(token,A->num_rows,A->pattern->ptr,A->pattern->index,&non_zeros,&ops);                 DPSLDLT_Preprocess(token,A->numRows,A->pattern->ptr,A->pattern->index,&non_zeros,&ops);
112                 DPSLDLT_Factor(token,A->num_rows,A->pattern->ptr,A->pattern->index,A->val);                 DPSLDLT_Factor(token,A->numRows,A->pattern->ptr,A->pattern->index,A->val);
113                 if (options->verbose) printf("timing SCSL: Cholevsky factorization: %.4e sec (token = %d)\n",Paso_timer()-time0,token);                 if (options->verbose) printf("timing SCSL: Cholevsky factorization: %.4e sec (token = %d)\n",Paso_timer()-time0,token);
114            } else {            } else {
115                 /* DPSLDU_Ordering(token,reordering_method);(does not work)*/                 /* DPSLDU_Ordering(token,reordering_method);(does not work)*/
116                 DPSLDU_Preprocess(token,A->num_rows,A->pattern->ptr,A->pattern->index,&non_zeros,&ops);                 DPSLDU_Preprocess(token,A->numRows,A->pattern->ptr,A->pattern->index,&non_zeros,&ops);
117                 DPSLDU_Factor(token,A->num_rows,A->pattern->ptr,A->pattern->index,A->val);                 DPSLDU_Factor(token,A->numRows,A->pattern->ptr,A->pattern->index,A->val);
118                 if (options->verbose) printf("timing SCSL: LDU factorization: %.4e sec (token = %d)\n",Paso_timer()-time0,token);                 if (options->verbose) printf("timing SCSL: LDU factorization: %.4e sec (token = %d)\n",Paso_timer()-time0,token);
119            }            }
120         }         }

Legend:
Removed from v.1119  
changed lines
  Added in v.1120

  ViewVC Help
Powered by ViewVC 1.1.26