/[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 3911 by jfenwick, Thu Jun 14 01:01:03 2012 UTC revision 3981 by jfenwick, Fri Sep 21 02:47:54 2012 UTC
# Line 1  Line 1 
1    
2  /*******************************************************  /*****************************************************************************
3  *  *
4  * Copyright (c) 2003-2012 by University of Queensland  * Copyright (c) 2003-2012 by University of Queensland
5  * Earth Systems Science Computational Center (ESSCC)  * http://www.uq.edu.au
 * http://www.uq.edu.au/esscc  
6  *  *
7  * Primary Business: Queensland, Australia  * Primary Business: Queensland, Australia
8  * Licensed under the Open Software License version 3.0  * Licensed under the Open Software License version 3.0
9  * http://www.opensource.org/licenses/osl-3.0.php  * http://www.opensource.org/licenses/osl-3.0.php
10  *  *
11  *******************************************************/  * Development until 2012 by Earth Systems Science Computational Center (ESSCC)
12    * Development since 2012 by School of Earth Sciences
13    *
14    *****************************************************************************/
15    
16    
17  /**************************************************************/  /************************************************************************************/
18    
19  /* Paso: Transport solver with flux correction (L is row sum zero)  /* Paso: Transport solver with flux correction (L is row sum zero)
20   *   *
# Line 21  Line 23 
23   *  to return v(dt)   *  to return v(dt)
24   *   *
25  */  */
26  /**************************************************************/  /************************************************************************************/
27    
28  /* Author: l.gross@uq.edu.au */  /* Author: l.gross@uq.edu.au */
29    
30  /**************************************************************/  /************************************************************************************/
31    
32  #include "FCT_Solver.h"  #include "FCT_Solver.h"
33  #include "Preconditioner.h"  #include "Preconditioner.h"
# Line 319  err_t Paso_FCT_Solver_updateNL(Paso_FCT_ Line 321  err_t Paso_FCT_Solver_updateNL(Paso_FCT_
321      }              }        
322      Paso_FCT_FluxLimiter_setU_tilda(flux_limiter, b); /* u_tilda = m^{-1} b */      Paso_FCT_FluxLimiter_setU_tilda(flux_limiter, b); /* u_tilda = m^{-1} b */
323      /* u_tilda_connector is completed */      /* u_tilda_connector is completed */
324      /**********************************************************************************************************************/        /********************************************************************************************************************************************/  
325      /* calculate stopping criterium */      /* calculate stopping criterium */
326      norm_u_tilde=Paso_lsup(n, flux_limiter->u_tilde, flux_limiter->mpi_info);      norm_u_tilde=Paso_lsup(n, flux_limiter->u_tilde, flux_limiter->mpi_info);
327      ATOL= rtol * norm_u_tilde + atol ;      ATOL= rtol * norm_u_tilde + atol ;
# Line 591  void Paso_FCT_setAntiDiffusionFlux_linea Line 593  void Paso_FCT_setAntiDiffusionFlux_linea
593    
594  }  }
595    
596  /**************************************************************/  /************************************************************************************/
597    
598  /* Creates the low order transport matrix and stores its negative values  /* Creates the low order transport matrix and stores its negative values
599   * into the iteration_matrix except for the main diagonal which is stored   * into the iteration_matrix except for the main diagonal which is stored
# Line 729  void Paso_FCT_Solver_setMuPaLu(double* o Line 731  void Paso_FCT_Solver_setMuPaLu(double* o
731    }    }
732  }  }
733    
734  /* *************************************************************************************************************************** */  /* ************************************************************************************************************************************************* */
735    
736    
737    

Legend:
Removed from v.3911  
changed lines
  Added in v.3981

  ViewVC Help
Powered by ViewVC 1.1.26