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

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

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

trunk/paso/src/Solvers/Solver_ILU.c revision 483 by jgs, Thu Feb 2 02:10:15 2006 UTC trunk/paso/src/Solver_ILU.c revision 1028 by gross, Wed Mar 14 00:15:24 2007 UTC
# Line 1  Line 1 
1  /* $Id$ */  /* $Id$ */
2    
3    
4    /*
5    ********************************************************************************
6    *               Copyright   2006 by ACcESS MNRF                                *
7    *                                                                              *
8    *                 http://www.access.edu.au                                     *
9    *           Primary Business: Queensland, Australia                            *
10    *     Licensed under the Open Software License version 3.0             *
11    *        http://www.opensource.org/licenses/osl-3.0.php                        *
12    ********************************************************************************
13    */
14    
15  /**************************************************************/  /**************************************************************/
16    
17  /* Paso: ILU preconditioner with reordering                 */  /* Paso: ILU preconditioner with reordering                 */
# Line 37  void Paso_Solver_ILU_free(Paso_Solver_IL Line 49  void Paso_Solver_ILU_free(Paso_Solver_IL
49  Paso_Solver_ILU* Paso_Solver_getILU(Paso_SystemMatrix * A,bool_t verbose) {  Paso_Solver_ILU* Paso_Solver_getILU(Paso_SystemMatrix * A,bool_t verbose) {
50    dim_t n=A->num_rows;    dim_t n=A->num_rows;
51    dim_t n_block=A->row_block_size;    dim_t n_block=A->row_block_size;
52    index_t num_colors=0;    index_t num_colors=0, *mis_marker=NULL;
53    register double A11,A12,A13,A21,A22,A23,A31,A32,A33,D;    register double A11,A12,A13,A21,A22,A23,A31,A32,A33,D;
54    register double mainA11,mainA12,mainA13,mainA21,mainA22,mainA23,mainA31,mainA32,mainA33;    register double mainA11,mainA12,mainA13,mainA21,mainA22,mainA23,mainA31,mainA32,mainA33;
55    register double S11,S12,S13,S21,S22,S23,S31,S32,S33;    register double S11,S12,S13,S21,S22,S23,S31,S32,S33;
# Line 46  Paso_Solver_ILU* Paso_Solver_getILU(Paso Line 58  Paso_Solver_ILU* Paso_Solver_getILU(Paso
58    /* allocations: */      /* allocations: */  
59    Paso_Solver_ILU* out=MEMALLOC(1,Paso_Solver_ILU);    Paso_Solver_ILU* out=MEMALLOC(1,Paso_Solver_ILU);
60    if (Paso_checkPtr(out)) return NULL;    if (Paso_checkPtr(out)) return NULL;
61    index_t* mis_marker=TMPMEMALLOC(n,index_t);    mis_marker=TMPMEMALLOC(n,index_t);
62    out->colorOf=MEMALLOC(n,index_t);    out->colorOf=MEMALLOC(n,index_t);
63    out->factors=MEMALLOC(A->len,double);    out->factors=MEMALLOC(A->len,double);
64    out->main_iptr=MEMALLOC(n,index_t);    out->main_iptr=MEMALLOC(n,index_t);

Legend:
Removed from v.483  
changed lines
  Added in v.1028

  ViewVC Help
Powered by ViewVC 1.1.26