/[escript]/trunk/paso/src/FCT_Solver.cpp
ViewVC logotype

Diff of /trunk/paso/src/FCT_Solver.cpp

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

revision 4657 by jfenwick, Thu Feb 6 06:12:20 2014 UTC revision 4800 by caltinay, Wed Mar 26 01:50:04 2014 UTC
# Line 458  void Paso_FCT_setAntiDiffusionFlux_CN(Pa Line 458  void Paso_FCT_setAntiDiffusionFlux_CN(Pa
458    const double *remote_u=Paso_Coupler_borrowRemoteData(u_coupler);    const double *remote_u=Paso_Coupler_borrowRemoteData(u_coupler);
459    const double *remote_u_old=Paso_Coupler_borrowRemoteData(u_old_coupler);    const double *remote_u_old=Paso_Coupler_borrowRemoteData(u_old_coupler);
460    const double dt_half= dt/2;    const double dt_half= dt/2;
461    const Paso_SystemMatrixPattern *pattern=fct->iteration_matrix->pattern;    const paso::SystemMatrixPattern *pattern=fct->iteration_matrix->pattern;
462    const dim_t n=Paso_SystemMatrix_getTotalNumRows(fct->iteration_matrix);    const dim_t n=Paso_SystemMatrix_getTotalNumRows(fct->iteration_matrix);
463    
464    #pragma omp parallel for schedule(static) private(i, iptr_ij)    #pragma omp parallel for schedule(static) private(i, iptr_ij)
# Line 503  void Paso_FCT_setAntiDiffusionFlux_BE(Pa Line 503  void Paso_FCT_setAntiDiffusionFlux_BE(Pa
503    const double *u_old= Paso_Coupler_borrowLocalData(u_old_coupler);    const double *u_old= Paso_Coupler_borrowLocalData(u_old_coupler);
504    const double *remote_u=Paso_Coupler_borrowRemoteData(u_coupler);    const double *remote_u=Paso_Coupler_borrowRemoteData(u_coupler);
505    const double *remote_u_old=Paso_Coupler_borrowRemoteData(u_old_coupler);    const double *remote_u_old=Paso_Coupler_borrowRemoteData(u_old_coupler);
506    const Paso_SystemMatrixPattern *pattern=fct->iteration_matrix->pattern;    const paso::SystemMatrixPattern *pattern=fct->iteration_matrix->pattern;
507    const dim_t  n=Paso_SystemMatrix_getTotalNumRows(fct->iteration_matrix);    const dim_t  n=Paso_SystemMatrix_getTotalNumRows(fct->iteration_matrix);
508    
509    #pragma omp parallel for schedule(static) private(i, iptr_ij)    #pragma omp parallel for schedule(static) private(i, iptr_ij)
# Line 557  void Paso_FCT_setAntiDiffusionFlux_linea Line 557  void Paso_FCT_setAntiDiffusionFlux_linea
557    const double *u_old= Paso_Coupler_borrowLocalData(u_old_coupler);    const double *u_old= Paso_Coupler_borrowLocalData(u_old_coupler);
558    const double *remote_u_tilde=Paso_Coupler_borrowRemoteData(u_tilde_coupler);    const double *remote_u_tilde=Paso_Coupler_borrowRemoteData(u_tilde_coupler);
559    const double *remote_u_old=Paso_Coupler_borrowRemoteData(u_old_coupler);    const double *remote_u_old=Paso_Coupler_borrowRemoteData(u_old_coupler);
560    const Paso_SystemMatrixPattern *pattern=fct->iteration_matrix->pattern;    const paso::SystemMatrixPattern *pattern=fct->iteration_matrix->pattern;
561    const dim_t n=Paso_SystemMatrix_getTotalNumRows(fct->iteration_matrix);    const dim_t n=Paso_SystemMatrix_getTotalNumRows(fct->iteration_matrix);
562    
563    #pragma omp parallel for schedule(static) private(i, iptr_ij)    #pragma omp parallel for schedule(static) private(i, iptr_ij)
# Line 625  void Paso_FCT_setLowOrderOperator(Paso_T Line 625  void Paso_FCT_setLowOrderOperator(Paso_T
625    }    }
626    
627    if (Esys_noError()) {    if (Esys_noError()) {
628        const Paso_SystemMatrixPattern *pattern=fc->iteration_matrix->pattern;        const paso::SystemMatrixPattern *pattern=fc->iteration_matrix->pattern;
629        const dim_t n=Paso_SystemMatrix_getTotalNumRows(fc->iteration_matrix);        const dim_t n=Paso_SystemMatrix_getTotalNumRows(fc->iteration_matrix);
630        #pragma omp parallel for private(i, iptr_ij, iptr_ji)  schedule(static)        #pragma omp parallel for private(i, iptr_ij, iptr_ji)  schedule(static)
631        for (i = 0; i < n; ++i) {        for (i = 0; i < n; ++i) {
# Line 693  void Paso_FCT_Solver_setMuPaLu(double* o Line 693  void Paso_FCT_Solver_setMuPaLu(double* o
693                                const Paso_SystemMatrix *L)                                const Paso_SystemMatrix *L)
694  {  {
695    dim_t i;    dim_t i;
696    const Paso_SystemMatrixPattern *pattern = L->pattern;    const paso::SystemMatrixPattern *pattern = L->pattern;
697    const double *u=Paso_Coupler_borrowLocalData(u_coupler);    const double *u=Paso_Coupler_borrowLocalData(u_coupler);
698    const double *remote_u=Paso_Coupler_borrowRemoteData(u_coupler);    const double *remote_u=Paso_Coupler_borrowRemoteData(u_coupler);
699    index_t iptr_ij;    index_t iptr_ij;

Legend:
Removed from v.4657  
changed lines
  Added in v.4800

  ViewVC Help
Powered by ViewVC 1.1.26