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

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

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

revision 2239 by artak, Thu Dec 18 00:08:58 2008 UTC revision 2240 by artak, Wed Feb 4 01:17:06 2009 UTC
# Line 154  Paso_Solver_AMG* Paso_Solver_getAMG(Paso Line 154  Paso_Solver_AMG* Paso_Solver_getAMG(Paso
154             if( Paso_noError()) {             if( Paso_noError()) {
155                /* if there are no nodes in the coarse level there is no more work to do */                /* if there are no nodes in the coarse level there is no more work to do */
156                out->n_C=n-out->n_F;                out->n_C=n-out->n_F;
157                if (level<3) {                /*if (level<3) {*/
158                 /*if (out->n_F>500) {*/                 if (out->n_F>500) {
159                     out->rows_in_C=MEMALLOC(out->n_C,index_t);                     out->rows_in_C=MEMALLOC(out->n_C,index_t);
160                     out->mask_C=MEMALLOC(n,index_t);                     out->mask_C=MEMALLOC(n,index_t);
161                     if (! (Paso_checkPtr(out->mask_C) || Paso_checkPtr(out->rows_in_C) ) ) {                     if (! (Paso_checkPtr(out->mask_C) || Paso_checkPtr(out->rows_in_C) ) ) {
# Line 249  Paso_Solver_AMG* Paso_Solver_getAMG(Paso Line 249  Paso_Solver_AMG* Paso_Solver_getAMG(Paso
249    if (Paso_noError()) {    if (Paso_noError()) {
250        if (verbose) {        if (verbose) {
251           printf("AMG: %d unknowns eliminated. %d left.\n",out->n_F,n-out->n_F);           printf("AMG: %d unknowns eliminated. %d left.\n",out->n_F,n-out->n_F);
252           if (level<3) {           /*if (level<3) {*/
253           /*if (out->n_F<500) {*/           if (out->n_F<500) {
254              printf("timing: AMG: MIS/reordering/elemination : %e/%e/%e\n",time2,time0,time1);              printf("timing: AMG: MIS/reordering/elemination : %e/%e/%e\n",time2,time0,time1);
255           } else {           } else {
256              printf("timing: AMG: MIS: %e\n",time2);              printf("timing: AMG: MIS: %e\n",time2);
# Line 294  void Paso_Solver_solveAMG(Paso_Solver_AM Line 294  void Paso_Solver_solveAMG(Paso_Solver_AM
294       double *x0=MEMALLOC(amg->n,double);       double *x0=MEMALLOC(amg->n,double);
295       double time0=0;       double time0=0;
296            
297       if (amg->level==3) {       /*if (amg->level==3) {*/
298       /*if (amg->n_F<=500) {*/       if (amg->n_F<=500) {
299        time0=Paso_timer();        time0=Paso_timer();
300                    
301          Paso_Solver_solveJacobi(amg->GS,x,b);          Paso_Solver_solveJacobi(amg->GS,x,b);
Legend:
Removed from v.2239  
changed lines
  Added in v.2240
  ViewVC Help
Powered by ViewVC 1.1.26