dcoomml.c

来自「用于求解大型稀疏线性方程组Ax=b的数值计算库.」· C语言 代码 · 共 562 行 · 第 1/2 页

      rowi =  indx[j];
      coli =  jndx[j];
      if ( rowi !=  coli ) {
        c[rowi] +=  b[coli] * val[j];
        c[coli] -=  b[rowi] * val[j];
      }
    }
    c += ldc; b += ldb;                                
  }                                                    
}



void COO_MatMult_CaABC_double(
                 const int m, const int n, const int k, const double alpha,
                 const double *val, const int *indx, const int *jndx,
                 const int nnz, 
                 const double *b, const int ldb, 
                 double *c, const int ldc, 
                 const int ind_base)
{
  int i,j;
  double *pc=c;
  int l;                                               
  c-=ind_base;
  b-=ind_base;


  for (l=0;l!=n;l++) {                                 
    for (j=0;j!=nnz;j++)
      c[indx[j]] += alpha * b[jndx[j]] * val[j];
    c += ldc; b += ldb;                                
  }                                                    
}


void COOsymm_MatMult_CaABC_double(
                 const int m, const int n, const int k, const double alpha,
                 const double *val, const int *indx, const int *jndx,
                 const int nnz, 
                 const double *b, const int ldb, 
                 double *c, const int ldc, 
                 const int ind_base)
{
  double *pc=c;
  int l;                                               
  int i,j;
  int rowi,coli;
  c-=ind_base;
  b-=ind_base;


  for (l=0;l!=n;l++) {                                 
    for (j=0;j!=nnz;j++) {
      rowi =  indx[j];
      coli =  jndx[j];
      if ( rowi ==  coli )
        c[rowi] += alpha * b[coli] * val[j];
      else {
        c[rowi] += alpha * b[coli] * val[j];
        c[coli] += alpha * b[rowi] * val[j];
      }
    }
    c += ldc; b += ldb;                                
  }                                                    
}


void COOskew_MatMult_CaABC_double(
                 const int m, const int n, const int k, const double alpha,
                 const double *val, const int *indx, const int *jndx,
                 const int nnz, 
                 const double *b, const int ldb, 
                 double *c, const int ldc, 
                 const int ind_base)
{
  double *pc=c;
  int l;                                               
  int i,j;
  int rowi,coli;
  c-=ind_base;
  b-=ind_base;


  for (l=0;l!=n;l++) {                                 
    for (j=0;j!=nnz;j++) {
      rowi =  indx[j];
      coli =  jndx[j];
      if ( rowi !=  coli ) {
        c[rowi] += alpha * b[coli] * val[j];
        c[coli] -= alpha * b[rowi] * val[j];
      }
    }
    c += ldc; b += ldb;                                
  }                                                    
}



void COO_MatMult_CABbC_double(
                 const int m, const int n, const int k, 
                 const double *val, const int *indx, const int *jndx,
                 const int nnz, 
                 const double *b, const int ldb, const double beta,
                 double *c, const int ldc, 
                 const int ind_base)
{
  int i,j;
  double *pc=c;
  int l;                                               
  c-=ind_base;
  b-=ind_base;

  for (l=0;l!=n;l++)                            
    for (i=0;i!=m;i++) *pc++ *= beta;                  

  for (l=0;l!=n;l++) {                                 
    for (j=0;j!=nnz;j++)
      c[indx[j]] +=  b[jndx[j]] * val[j];
    c += ldc; b += ldb;                                
  }                                                    
}


void COOsymm_MatMult_CABbC_double(
                 const int m, const int n, const int k, 
                 const double *val, const int *indx, const int *jndx,
                 const int nnz, 
                 const double *b, const int ldb, const double beta,
                 double *c, const int ldc, 
                 const int ind_base)
{
  double *pc=c;
  int l;                                               
  int i,j;
  int rowi,coli;
  c-=ind_base;
  b-=ind_base;

  for (l=0;l!=n;l++)                            
    for (i=0;i!=m;i++) *pc++ *= beta;                  

  for (l=0;l!=n;l++) {                                 
    for (j=0;j!=nnz;j++) {
      rowi =  indx[j];
      coli =  jndx[j];
      if ( rowi ==  coli )
        c[rowi] +=  b[coli] * val[j];
      else {
        c[rowi] +=  b[coli] * val[j];
        c[coli] +=  b[rowi] * val[j];
      }
    }
    c += ldc; b += ldb;                                
  }                                                    
}


void COOskew_MatMult_CABbC_double(
                 const int m, const int n, const int k, 
                 const double *val, const int *indx, const int *jndx,
                 const int nnz, 
                 const double *b, const int ldb, const double beta,
                 double *c, const int ldc, 
                 const int ind_base)
{
  double *pc=c;
  int l;                                               
  int i,j;
  int rowi,coli;
  c-=ind_base;
  b-=ind_base;

  for (l=0;l!=n;l++)                            
    for (i=0;i!=m;i++) *pc++ *= beta;                  

  for (l=0;l!=n;l++) {                                 
    for (j=0;j!=nnz;j++) {
      rowi =  indx[j];
      coli =  jndx[j];
      if ( rowi !=  coli ) {
        c[rowi] +=  b[coli] * val[j];
        c[coli] -=  b[rowi] * val[j];
      }
    }
    c += ldc; b += ldb;                                
  }                                                    
}



void COO_MatMult_CaABbC_double(
                 const int m, const int n, const int k, const double alpha,
                 const double *val, const int *indx, const int *jndx,
                 const int nnz, 
                 const double *b, const int ldb, const double beta,
                 double *c, const int ldc, 
                 const int ind_base)
{
  int i,j;
  double *pc=c;
  int l;                                               
  c-=ind_base;
  b-=ind_base;

  for (l=0;l!=n;l++)                            
    for (i=0;i!=m;i++) *pc++ *= beta;                  

  for (l=0;l!=n;l++) {                                 
    for (j=0;j!=nnz;j++)
      c[indx[j]] += alpha * b[jndx[j]] * val[j];
    c += ldc; b += ldb;                                
  }                                                    
}


void COOsymm_MatMult_CaABbC_double(
                 const int m, const int n, const int k, const double alpha,
                 const double *val, const int *indx, const int *jndx,
                 const int nnz, 
                 const double *b, const int ldb, const double beta,
                 double *c, const int ldc, 
                 const int ind_base)
{
  double *pc=c;
  int l;                                               
  int i,j;
  int rowi,coli;
  c-=ind_base;
  b-=ind_base;

  for (l=0;l!=n;l++)                            
    for (i=0;i!=m;i++) *pc++ *= beta;                  

  for (l=0;l!=n;l++) {                                 
    for (j=0;j!=nnz;j++) {
      rowi =  indx[j];
      coli =  jndx[j];
      if ( rowi ==  coli )
        c[rowi] += alpha * b[coli] * val[j];
      else {
        c[rowi] += alpha * b[coli] * val[j];
        c[coli] += alpha * b[rowi] * val[j];
      }
    }
    c += ldc; b += ldb;                                
  }                                                    
}


void COOskew_MatMult_CaABbC_double(
                 const int m, const int n, const int k, const double alpha,
                 const double *val, const int *indx, const int *jndx,
                 const int nnz, 
                 const double *b, const int ldb, const double beta,
                 double *c, const int ldc, 
                 const int ind_base)
{
  double *pc=c;
  int l;                                               
  int i,j;
  int rowi,coli;
  c-=ind_base;
  b-=ind_base;

  for (l=0;l!=n;l++)                            
    for (i=0;i!=m;i++) *pc++ *= beta;                  

  for (l=0;l!=n;l++) {                                 
    for (j=0;j!=nnz;j++) {
      rowi =  indx[j];
      coli =  jndx[j];
      if ( rowi !=  coli ) {
        c[rowi] += alpha * b[coli] * val[j];
        c[coli] -= alpha * b[rowi] * val[j];
      }
    }
    c += ldc; b += ldb;                                
  }                                                    
}