26 int main (
int argc,
char **argv)
30 double dparm[DPARM_SIZE];
32 char *filename = NULL;
33 spmatrix_t *spm, spm2;
34 void *x, *b, *x0 = NULL;
51 &check, &scatter, &driver, &filename );
56 pastixInit( &pastix_data, MPI_COMM_WORLD, iparm, dparm );
61 spm = malloc(
sizeof( spmatrix_t ) );
63 rc = spmReadDriverDist( driver, filename, spm, MPI_COMM_WORLD );
66 rc = spmReadDriver( driver, filename, spm );
69 if ( rc != SPM_SUCCESS ) {
74 spmPrintInfo( spm, stdout );
76 rc = spmCheckAndCorrect( spm, &spm2 );
86 if ( spm->flttype == SpmPattern ) {
87 spmGenFakeValues( spm );
98 double normA = spmNorm( SpmFrobeniusNorm, spm );
99 spmScal( 1./normA, spm );
110 size = pastix_size_of( spm->flttype ) * spm->nexp * nrhs;
119 spmGenRHS( SpmRhsRndX, nrhs, spm, x0, spm->nexp, b, spm->nexp );
120 memcpy( x, b, size );
123 spmGenRHS( SpmRhsRndB, nrhs, spm, NULL, spm->nexp, x, spm->nexp );
126 spmScalMat( 1./normA, spm, nrhs, b, spm->nexp );
129 memcpy( b, x, size );
BEGIN_C_DECLS typedef int pastix_int_t
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 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_analyze(pastix_data_t *pastix_data, const spmatrix_t *spm)
Perform all the preprocessing steps: ordering, symbolic factorization, reordering,...
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.