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

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

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

revision 1562 by gross, Wed May 21 13:04:40 2008 UTC revision 1628 by phornby, Fri Jul 11 13:12:46 2008 UTC
# Line 77  Paso_FCTransportProblem* Paso_FCTranspor Line 77  Paso_FCTransportProblem* Paso_FCTranspor
77       Paso_SystemMatrixType matrix_type=MATRIX_FORMAT_DEFAULT+MATRIX_FORMAT_BLK1;  /* at the moment only block size 1 is supported */       Paso_SystemMatrixType matrix_type=MATRIX_FORMAT_DEFAULT+MATRIX_FORMAT_BLK1;  /* at the moment only block size 1 is supported */
78       Paso_FCTransportProblem* out=NULL;       Paso_FCTransportProblem* out=NULL;
79       dim_t n,i;       dim_t n,i;
80       index_t iptr,iptr_main,k;       index_t iptr,iptr_main;
81    
82       if ((theta<0.) || (theta >1.)) {       if ((theta<0.) || (theta >1.)) {
83          Paso_setError(TYPE_ERROR,"Paso_FCTransportProblem_alloc: theta needs to be between 0. and. 1.");          Paso_setError(TYPE_ERROR,"Paso_FCTransportProblem_alloc: theta needs to be between 0. and. 1.");
# Line 122  Paso_FCTransportProblem* Paso_FCTranspor Line 122  Paso_FCTransportProblem* Paso_FCTranspor
122                      out->u[i]=0.;                      out->u[i]=0.;
123                   }                   }
124                   /* identify the main diagonals */                   /* identify the main diagonals */
125                   #pragma omp for schedule(static) private(i,iptr,iptr_main,k)                   #pragma omp for schedule(static) private(i,iptr,iptr_main)
126                   for (i = 0; i < n; ++i) {                   for (i = 0; i < n; ++i) {
127                          iptr_main=pattern->mainPattern->ptr[0]-1;                          iptr_main=pattern->mainPattern->ptr[0]-1;
128                          for (iptr=pattern->mainPattern->ptr[i];iptr<pattern->mainPattern->ptr[i+1]; iptr++) {                          for (iptr=pattern->mainPattern->ptr[i];iptr<pattern->mainPattern->ptr[i+1]; iptr++) {

Legend:
Removed from v.1562  
changed lines
  Added in v.1628

  ViewVC Help
Powered by ViewVC 1.1.26