24 int main (
int argc,
char **argv)
28 double dparm[DPARM_SIZE];
30 char *filename = NULL;
31 spmatrix_t *spm, spm2;
32 void *x, *b, *x0 = NULL;
51 &check, &scatter, &driver, &filename );
57 pastixInit( &pastix_data, MPI_COMM_WORLD, iparm, dparm );
62 spm = malloc(
sizeof( spmatrix_t ) );
64 rc = spmReadDriverDist( driver, filename, spm, MPI_COMM_WORLD );
67 rc = spmReadDriver( driver, filename, spm );
70 if ( rc != SPM_SUCCESS ) {
75 spmPrintInfo( spm, stdout );
77 rc = spmCheckAndCorrect( spm, &spm2 );
87 if ( spm->flttype == SpmPattern ) {
88 spmGenFakeValues( spm );
112 double normA = spmNorm( SpmFrobeniusNorm, spm );
113 spmScal( 1./normA, spm );
124 size = pastix_size_of( spm->flttype ) * spm->nexp * nrhs;
133 spmGenRHS( SpmRhsRndX, nrhs, spm, x0, spm->nexp, b, spm->nexp );
134 memcpy( x, b, size );
137 spmGenRHS( SpmRhsRndB, nrhs, spm, NULL, spm->nexp, x, spm->nexp );
140 spmScalMat( 1./normA, spm, nrhs, b, spm->nexp );
143 memcpy( b, x, size );
BEGIN_C_DECLS typedef int pastix_int_t
int pastix_subtask_symbfact(pastix_data_t *pastix_data)
Computes the symbolic factorization step.
int pastix_subtask_order(pastix_data_t *pastix_data, const spmatrix_t *spm, pastix_order_t *myorder)
Computes the ordering of the given graph in parameters.
int pastix_subtask_blend(pastix_data_t *pastix_data)
Compute the proportional mapping and the final solver structure.
int pastix_subtask_reordering(pastix_data_t *pastix_data)
Apply the reordering step to compact off-diagonal blocks.
void pastixFinalize(pastix_data_t **pastix_data)
Finalize the solver instance.
void pastixInitParam(pastix_int_t *iparm, double *dparm)
Initialize the iparm and dparm arrays to their default values.
void pastixInit(pastix_data_t **pastix_data, PASTIX_Comm pastix_comm, pastix_int_t *iparm, double *dparm)
Initialize the solver instance.
@ DPARM_EPSILON_REFINEMENT
void pastixGetOptions(int argc, char **argv, pastix_int_t *iparm, double *dparm, int *check, int *scatter, spm_driver_t *driver, char **filename)
PaStiX helper function to read command line options in examples.
int pastixOrderAlloc(pastix_order_t *ordeptr, pastix_int_t vertnbr, pastix_int_t cblknbr)
Allocate the order structure.
void pastixOrderExit(pastix_order_t *ordeptr)
Free the arrays initialized in the order structure.
int pastix_task_refine(pastix_data_t *pastix_data, pastix_int_t n, pastix_int_t nrhs, void *B, pastix_int_t ldb, void *X, pastix_int_t ldx)
Perform iterative refinement.
int pastix_task_solve(pastix_data_t *pastix_data, pastix_int_t m, pastix_int_t nrhs, void *B, pastix_int_t ldb)
Solve the given problem.
int pastix_task_numfact(pastix_data_t *pastix_data, spmatrix_t *spm)
Perform all the numerical factorization steps: fill the internal block CSC and the solver matrix stru...
Main PaStiX data structure.