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

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

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

revision 2938 by gross, Mon Feb 8 06:48:19 2010 UTC revision 2939 by gross, Thu Feb 18 04:16:59 2010 UTC
# Line 203  void Paso_SolverFCT_solve(Paso_FCTranspo Line 203  void Paso_SolverFCT_solve(Paso_FCTranspo
203                       * now we solve the linear system to get the correction dt:                       * now we solve the linear system to get the correction dt:
204                       *                       *
205                       */                       */
 /* for (i=0;i<n;++i) printf("%d %f \n", i,z_m[i]); */  
206                       if (fctp->theta > 0) {                       if (fctp->theta > 0) {
207                            omega=1./(dt2*fctp->theta);                            omega=1./(dt2*fctp->theta);
208                            Paso_Solver_solvePreconditioner(fctp->iteration_matrix,du_m,z_m);                              Paso_Solver_solvePreconditioner(fctp->iteration_matrix,du_m,z_m);  
209                            Paso_Update(n,1.,u,omega,du_m);                            Paso_Update(n,1.,u,omega,du_m);
210  for (i = 0; i < n; ++i) printf("u %d = %e %e <= %e \n",i,u[i], du_m[i], z_m[i]);  /* for (i = 0; i < n; ++i) printf("u %d = %e %e <= %e \n",i,u[i], du_m[i], z_m[i]); */
211               } else {               } else {
212                            omega=1;                            omega=1;
213                            #pragma omp parallel for private(i,mass,rtmp)                            #pragma omp parallel for private(i,mass,rtmp)
# Line 368  void Paso_FCT_setUp(Paso_FCTransportProb Line 367  void Paso_FCT_setUp(Paso_FCTransportProb
367                if (ABS(u_tilde_i)>0) rtmp4+=sourceN[i]*factor/u_tilde_i;                if (ABS(u_tilde_i)>0) rtmp4+=sourceN[i]*factor/u_tilde_i;
368                fctp->iteration_matrix->mainBlock->val[fctp->main_iptr[i]]=rtmp4;                fctp->iteration_matrix->mainBlock->val[fctp->main_iptr[i]]=rtmp4;
369                uTilde[i]=u_tilde_i;                uTilde[i]=u_tilde_i;
370  printf("uTilde %d %e and %e : %e : %e %e\n",i,u_tilde_i, b[i], rtmp4, m, fctp->main_diagonal_low_order_transport_matrix[i]);  /* printf("uTilde %d %e and %e : %e : %e %e\n",i,u_tilde_i, b[i], rtmp4, m, fctp->main_diagonal_low_order_transport_matrix[i]); */
371           }           }
372           Performance_startMonitor(pp,PERFORMANCE_PRECONDITIONER_INIT);           Performance_startMonitor(pp,PERFORMANCE_PRECONDITIONER_INIT);
373           Paso_Solver_setPreconditioner(fctp->iteration_matrix,options);           Paso_Solver_setPreconditioner(fctp->iteration_matrix,options);

Legend:
Removed from v.2938  
changed lines
  Added in v.2939

  ViewVC Help
Powered by ViewVC 1.1.26