HYPRE_Int
hypre_SStructPScale( HYPRE_Complex alpha,
hypre_SStructPVector *py )
{
HYPRE_Int nvars = hypre_SStructPVectorNVars(py);
HYPRE_Int var;
for (var = 0; var < nvars; var++)
{
hypre_StructScale(alpha, hypre_SStructPVectorSVector(py, var));
}
return hypre_error_flag;
}
HYPRE_Int
hypre_SStructPCopy( hypre_SStructPVector *px,
hypre_SStructPVector *py )
{
HYPRE_Int nvars = hypre_SStructPVectorNVars(px);
HYPRE_Int var;
for (var = 0; var < nvars; var++)
{
hypre_StructCopy(hypre_SStructPVectorSVector(px, var),
hypre_SStructPVectorSVector(py, var));
}
return hypre_error_flag;
}
int
hypre_SStructPScale( double alpha,
hypre_SStructPVector *py )
{
int ierr = 0;
int nvars = hypre_SStructPVectorNVars(py);
int var;
for (var = 0; var < nvars; var++)
{
hypre_StructScale(alpha, hypre_SStructPVectorSVector(py, var));
}
return ierr;
}
int
hypre_SStructPAxpy( double alpha,
hypre_SStructPVector *px,
hypre_SStructPVector *py )
{
int ierr = 0;
int nvars = hypre_SStructPVectorNVars(px);
int var;
for (var = 0; var < nvars; var++)
{
hypre_StructAxpy(alpha,
hypre_SStructPVectorSVector(px, var),
hypre_SStructPVectorSVector(py, var));
}
return ierr;
}
HYPRE_Int
hypre_SStructPartialPCopy( hypre_SStructPVector *px,
hypre_SStructPVector *py,
hypre_BoxArrayArray **array_boxes )
{
HYPRE_Int nvars = hypre_SStructPVectorNVars(px);
hypre_BoxArrayArray *boxes;
HYPRE_Int var;
for (var = 0; var < nvars; var++)
{
boxes= array_boxes[var];
hypre_StructPartialCopy(hypre_SStructPVectorSVector(px, var),
hypre_SStructPVectorSVector(py, var),
boxes);
}
return hypre_error_flag;
}
HYPRE_Int
hypre_SStructPVectorSetRandomValues( hypre_SStructPVector *pvector, HYPRE_Int seed )
{
HYPRE_Int ierr = 0;
HYPRE_Int nvars = hypre_SStructPVectorNVars(pvector);
hypre_StructVector *svector;
HYPRE_Int var;
srand( seed );
for (var = 0; var < nvars; var++)
{
svector = hypre_SStructPVectorSVector(pvector, var);
seed = rand();
hypre_StructVectorSetRandomValues(svector, seed);
}
return ierr;
}
HYPRE_Int
hypre_SStructPOverlapInnerProd( hypre_SStructPVector *px,
hypre_SStructPVector *py,
double *presult_ptr )
{
HYPRE_Int nvars = hypre_SStructPVectorNVars(px);
double presult;
double sresult;
HYPRE_Int var;
presult = 0.0;
for (var = 0; var < nvars; var++)
{
sresult = hypre_StructOverlapInnerProd(hypre_SStructPVectorSVector(px, var),
hypre_SStructPVectorSVector(py, var));
presult += sresult;
}
*presult_ptr = presult;
return hypre_error_flag;
}
HYPRE_Int
hypre_SStructPInnerProd( hypre_SStructPVector *px,
hypre_SStructPVector *py,
HYPRE_Real *presult_ptr )
{
HYPRE_Int nvars = hypre_SStructPVectorNVars(px);
HYPRE_Real presult;
HYPRE_Real sresult;
HYPRE_Int var;
presult = 0.0;
for (var = 0; var < nvars; var++)
{
sresult = hypre_StructInnerProd(hypre_SStructPVectorSVector(px, var),
hypre_SStructPVectorSVector(py, var));
presult += sresult;
}
*presult_ptr = presult;
return hypre_error_flag;
}
HYPRE_Int
hypre_FacSemiRestrictSetup2( void *fac_restrict_vdata,
hypre_SStructVector *r,
HYPRE_Int part_crse,
HYPRE_Int part_fine,
hypre_SStructPVector *rc,
hypre_Index rfactors )
{
HYPRE_Int ierr = 0;
hypre_FacSemiRestrictData2 *fac_restrict_data = fac_restrict_vdata;
MPI_Comm comm= hypre_SStructPVectorComm(rc);
hypre_CommInfo *comm_info;
hypre_CommPkg **interlevel_comm;
hypre_SStructPVector *rf= hypre_SStructVectorPVector(r, part_fine);
hypre_StructVector *s_rc, *s_cvector;
hypre_SStructPGrid *pgrid;
hypre_SStructPVector *fgrid_cvectors;
hypre_SStructPGrid *fgrid_coarsen;
hypre_BoxArrayArray **identity_arrayboxes;
hypre_BoxArrayArray **fullwgt_ownboxes;
hypre_BoxArrayArray **fullwgt_sendboxes;
hypre_BoxArray *boxarray;
hypre_BoxArray *tmp_boxarray, *intersect_boxes;
HYPRE_Int ***own_cboxnums;
hypre_BoxArrayArray **send_boxes, *send_rboxes;
HYPRE_Int ***send_processes;
HYPRE_Int ***send_remote_boxnums;
hypre_BoxArrayArray **recv_boxes, *recv_rboxes;
HYPRE_Int ***recv_processes;
HYPRE_Int ***recv_remote_boxnums;
hypre_BoxManager *boxman;
hypre_BoxManEntry **boxman_entries;
HYPRE_Int nboxman_entries;
hypre_Box box, scaled_box;
hypre_Index zero_index, index, ilower, iupper;
HYPRE_Int ndim= hypre_SStructVectorNDim(r);
HYPRE_Int myproc, proc;
HYPRE_Int nvars, vars;
HYPRE_Int num_values;
HYPRE_Int i, cnt1, cnt2;
HYPRE_Int fi, ci;
hypre_MPI_Comm_rank(comm, &myproc);
hypre_ClearIndex(zero_index);
nvars= hypre_SStructPVectorNVars(rc);
(fac_restrict_data -> nvars)= nvars;
hypre_CopyIndex(rfactors, (fac_restrict_data -> stride));
for (i= ndim; i< 3; i++)
{
rfactors[i]= 1;
}
/* work vector for storing the fullweighted fgrid boxes */
hypre_SStructPGridCreate(hypre_SStructPVectorComm(rf), ndim, &fgrid_coarsen);
pgrid= hypre_SStructPVectorPGrid(rf);
for (vars= 0; vars< nvars; vars++)
{
boxarray= hypre_StructGridBoxes(hypre_SStructPGridSGrid(pgrid, vars));
hypre_ForBoxI(fi, boxarray)
{
hypre_CopyBox(hypre_BoxArrayBox(boxarray, fi), &box);
hypre_StructMapFineToCoarse(hypre_BoxIMin(&box), zero_index,
rfactors, hypre_BoxIMin(&box));
hypre_StructMapFineToCoarse(hypre_BoxIMax(&box), zero_index,
rfactors, hypre_BoxIMax(&box));
hypre_SStructPGridSetExtents(fgrid_coarsen,
hypre_BoxIMin(&box),
hypre_BoxIMax(&box));
}
}
