C++ (Cpp) KppDecomp示例

示例#1

文件： saprc99_Integrator.c 项目： jlinford/fixedgrid

/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/   
void DecompTemplate( double A[], int Pivot[], int* ising )
/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  
        Template for the LU decomposition   
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/   
{   
   *ising = KppDecomp ( A );
  /*~~~> Note: for a full matrix use Lapack:
      DGETRF( 74, 74, A, 74, Pivot, ising ) */
    
   Ndec++;

}  /*  DecompTemplate */

示例#2

文件： rosenbrock.c 项目： CEKrause/WRFV_3.7.1

/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/   
void DecompTemplate( KPP_REAL A[], int Pivot[], int* ising )
/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  
        Template for the LU decomposition   
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/   
{   
   *ising = KppDecomp ( A );
  /*~~~> Note: for a full matrix use Lapack:
      DGETRF( KPP_NVAR, KPP_NVAR, A, KPP_NVAR, Pivot, ising ) */
    
   Ndec++;

}  /*  DecompTemplate */

示例#3

文件： saprc99_Integrator.c 项目： jlinford/fixedgrid

     
  /*~~~> ros_ELO  = estimator of local order - the minimum between the
!    main and the embedded scheme orders plus 1 */
    *ros_ELO = (double)4.0;
     
}  /*  Rodas4 */

   

/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/   
void DecompTemplate( double A[], int Pivot[] __attribute__((unused)), int* ising )
/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  
        Template for the LU decomposition   
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/   
{   
   *ising = KppDecomp ( A );
  /*~~~> Note: for a full matrix use Lapack:
      DGETRF( 74, 74, A, 74, Pivot, ising ) */
    
   Ndec++;

}  /*  DecompTemplate */
 
/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/   
 void SolveTemplate( double A[], int Pivot[] __attribute__((unused)), double b[] )
/*~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~  
     Template for the forward/backward substitution (using pre-computed LU decomposition)   
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~*/   
{   
   KppSolve( A, b );
  /*~~~> Note: for a full matrix use Lapack: