HYPRE_Int
hypre_SStructScale( HYPRE_Complex alpha,
hypre_SStructVector *y )
{
HYPRE_Int nparts = hypre_SStructVectorNParts(y);
HYPRE_Int part;
HYPRE_Int y_object_type= hypre_SStructVectorObjectType(y);
if (y_object_type == HYPRE_SSTRUCT)
{
for (part = 0; part < nparts; part++)
{
hypre_SStructPScale(alpha, hypre_SStructVectorPVector(y, part));
}
}
else if (y_object_type == HYPRE_PARCSR)
{
hypre_ParVector *y_par;
hypre_SStructVectorConvert(y, &y_par);
hypre_ParVectorScale(alpha, y_par);
}
return hypre_error_flag;
}
int
hypre_SStructScale( double alpha,
hypre_SStructVector *y )
{
int ierr = 0;
int nparts = hypre_SStructVectorNParts(y);
int part;
int y_object_type= hypre_SStructVectorObjectType(y);
if (y_object_type == HYPRE_SSTRUCT)
{
for (part = 0; part < nparts; part++)
{
hypre_SStructPScale(alpha, hypre_SStructVectorPVector(y, part));
}
}
else if (y_object_type == HYPRE_PARCSR)
{
hypre_ParVector *y_par;
hypre_SStructVectorConvert(y, &y_par);
hypre_ParVectorScale(alpha, y_par);
}
return ierr;
}
HYPRE_Int
hypre_SStructInnerProd( hypre_SStructVector *x,
hypre_SStructVector *y,
double *result_ptr )
{
HYPRE_Int nparts = hypre_SStructVectorNParts(x);
double result;
double presult;
HYPRE_Int part;
HYPRE_Int x_object_type= hypre_SStructVectorObjectType(x);
HYPRE_Int y_object_type= hypre_SStructVectorObjectType(y);
if (x_object_type != y_object_type)
{
hypre_error_in_arg(2);
hypre_error_in_arg(3);
return hypre_error_flag;
}
result = 0.0;
if ( (x_object_type == HYPRE_SSTRUCT) || (x_object_type == HYPRE_STRUCT) )
{
for (part = 0; part < nparts; part++)
{
hypre_SStructPInnerProd(hypre_SStructVectorPVector(x, part),
hypre_SStructVectorPVector(y, part), &presult);
result += presult;
}
}
else if (x_object_type == HYPRE_PARCSR)
{
hypre_ParVector *x_par;
hypre_ParVector *y_par;
hypre_SStructVectorConvert(x, &x_par);
hypre_SStructVectorConvert(y, &y_par);
result= hypre_ParVectorInnerProd(x_par, y_par);
}
*result_ptr = result;
return hypre_error_flag;
}
int
hypre_SStructAxpy( double alpha,
hypre_SStructVector *x,
hypre_SStructVector *y )
{
int ierr = 0;
int nparts = hypre_SStructVectorNParts(x);
int part;
int x_object_type= hypre_SStructVectorObjectType(x);
int y_object_type= hypre_SStructVectorObjectType(y);
if (x_object_type != y_object_type)
{
printf("vector object types different- cannot compute Axpy\n");
return ierr;
}
if (x_object_type == HYPRE_SSTRUCT)
{
for (part = 0; part < nparts; part++)
{
hypre_SStructPAxpy(alpha,
hypre_SStructVectorPVector(x, part),
hypre_SStructVectorPVector(y, part));
}
}
else if (x_object_type == HYPRE_PARCSR)
{
hypre_ParVector *x_par;
hypre_ParVector *y_par;
hypre_SStructVectorConvert(x, &x_par);
hypre_SStructVectorConvert(y, &y_par);
hypre_ParVectorAxpy(alpha, x_par, y_par);
}
return ierr;
}
HYPRE_Int
hypre_SStructCopy( hypre_SStructVector *x,
hypre_SStructVector *y )
{
HYPRE_Int nparts = hypre_SStructVectorNParts(x);
HYPRE_Int part;
HYPRE_Int x_object_type= hypre_SStructVectorObjectType(x);
HYPRE_Int y_object_type= hypre_SStructVectorObjectType(y);
if (x_object_type != y_object_type)
{
hypre_error_in_arg(2);
hypre_error_in_arg(3);
return hypre_error_flag;
}
if (x_object_type == HYPRE_SSTRUCT)
{
for (part = 0; part < nparts; part++)
{
hypre_SStructPCopy(hypre_SStructVectorPVector(x, part),
hypre_SStructVectorPVector(y, part));
}
}
else if (x_object_type == HYPRE_PARCSR)
{
hypre_ParVector *x_par;
hypre_ParVector *y_par;
hypre_SStructVectorConvert(x, &x_par);
hypre_SStructVectorConvert(y, &y_par);
hypre_ParVectorCopy(x_par, y_par);
}
return hypre_error_flag;
}
HYPRE_Int
hypre_SStructMatvecCompute( void *matvec_vdata,
HYPRE_Complex alpha,
hypre_SStructMatrix *A,
hypre_SStructVector *x,
HYPRE_Complex beta,
hypre_SStructVector *y )
{
hypre_SStructMatvecData *matvec_data = matvec_vdata;
HYPRE_Int nparts = (matvec_data -> nparts);
void **pmatvec_data = (matvec_data -> pmatvec_data);
void *pdata;
hypre_SStructPMatrix *pA;
hypre_SStructPVector *px;
hypre_SStructPVector *py;
hypre_ParCSRMatrix *parcsrA = hypre_SStructMatrixParCSRMatrix(A);
hypre_ParVector *parx;
hypre_ParVector *pary;
HYPRE_Int part;
HYPRE_Int x_object_type= hypre_SStructVectorObjectType(x);
HYPRE_Int A_object_type= hypre_SStructMatrixObjectType(A);
if (x_object_type != A_object_type)
{
hypre_error_in_arg(2);
hypre_error_in_arg(3);
return hypre_error_flag;
}
if ( (x_object_type == HYPRE_SSTRUCT) || (x_object_type == HYPRE_STRUCT) )
{
/* do S-matrix computations */
for (part = 0; part < nparts; part++)
{
pdata = pmatvec_data[part];
pA = hypre_SStructMatrixPMatrix(A, part);
px = hypre_SStructVectorPVector(x, part);
py = hypre_SStructVectorPVector(y, part);
hypre_SStructPMatvecCompute(pdata, alpha, pA, px, beta, py);
}
if (x_object_type == HYPRE_SSTRUCT)
{
/* do U-matrix computations */
/* GEC1002 the data chunk pointed by the local-parvectors
* inside the semistruct vectors x and y is now identical to the
* data chunk of the structure vectors x and y. The role of the function
* convert is to pass the addresses of the data chunk
* to the parx and pary. */
hypre_SStructVectorConvert(x, &parx);
hypre_SStructVectorConvert(y, &pary);
hypre_ParCSRMatrixMatvec(alpha, parcsrA, parx, 1.0, pary);
/* dummy functions since there is nothing to restore */
hypre_SStructVectorRestore(x, NULL);
hypre_SStructVectorRestore(y, pary);
parx = NULL;
}
}
else if (x_object_type == HYPRE_PARCSR)
{
hypre_SStructVectorConvert(x, &parx);
hypre_SStructVectorConvert(y, &pary);
hypre_ParCSRMatrixMatvec(alpha, parcsrA, parx, beta, pary);
hypre_SStructVectorRestore(x, NULL);
hypre_SStructVectorRestore(y, pary);
parx = NULL;
}
return hypre_error_flag;
}
