ItrSmootherCSRcr.c

Go to the documentation of this file.

 
#include <math.h>
 
#include "fasp.h"
#include "fasp_functs.h"
 
/*---------------------------------*/
/*--      Public Functions       --*/
/*---------------------------------*/
 
void fasp_smoother_dcsr_gscr (INT   pt,
                              INT   n,
                              REAL *u,
                              INT  *ia,
                              INT  *ja,
                              REAL *a,
                              REAL *b, 
                              INT   L,
                              INT  *CF)
{ 
    INT i,j,k,l;
    REAL t, d=0;
    
    for (l=0;l<L;++l) {
        for (i=0;i<n;++i) {
            if (CF[i] == pt) { 
                t=b[i];
                for (k=ia[i];k<ia[i+1];++k) {
                    j=ja[k];
                    if (CF[j] == pt) {
                        if (i!=j) {
                            t-=a[k]*u[j]; 
                        }
                        else { 
                            d=a[k];
                        }
                        if (ABS(d)>SMALLREAL) { 
                            u[i]=t/d;
                        }
                        else {
                            printf("### ERROR: Diagonal entry_%d (%e) close to 0!\n",
                                   i, d);
                            fasp_chkerr(ERROR_MISC, __FUNCTION__);
                        }
                    }
                }
            }
            else {
                u[i]=0.e0;
            }
        } 
    }
}
 
/*---------------------------------*/
/*--        End of File          --*/
/*---------------------------------*/