diff --git a/z_spm_genrhs.c b/z_spm_genrhs.c index 0f8ad1195ca813bd94e2f48c0f59d81602c739ad..25d0e17d0423e71e5872bfe2c58f4c6285f4fc32 100644 --- a/z_spm_genrhs.c +++ b/z_spm_genrhs.c @@ -192,6 +192,9 @@ z_spmGenRHS( int type, int nrhs, if( (nrhs > 1) && (ldb < spm->n) ) return PASTIX_ERR_BADPARAMETER; + if( spm->dof != 1 ) + return PASTIX_ERR_BADPARAMETER; + if (nrhs == 1) { ldb = spm->n; ldx = spm->n; @@ -342,6 +345,9 @@ z_spmCheckAxb( int nrhs, double backward, forward, eps; int failure = 0; + assert( spm->nexp == spm->n ); + assert( spm->dof == 1 ); + eps = LAPACKE_dlamch('e'); /**