diff --git a/src/z_spm_genrhs.c b/src/z_spm_genrhs.c index db4c626d1179cd1111dc637a4900f4f4955e699c..ce851855537715d1e034d2c4bd8b20d954e3543b 100644 --- a/src/z_spm_genrhs.c +++ b/src/z_spm_genrhs.c @@ -57,8 +57,9 @@ Rnd64_jump(unsigned long long int n, unsigned long long int seed ) { ran = seed; for (i = 0; n; n >>= 1, ++i) { - if (n & 1) + if (n & 1) { ran = a_k * ran + c_k; + } c_k *= (a_k + 1); a_k *= a_k; } diff --git a/src/z_spm_reduce_rhs.c b/src/z_spm_reduce_rhs.c index 30135468d800197a1807961461380454109cb28d..fd54e61888f68dc737db4a035c0bd103aac782f7 100644 --- a/src/z_spm_reduce_rhs.c +++ b/src/z_spm_reduce_rhs.c @@ -62,12 +62,15 @@ z_spmReduceRHS( const spmatrix_t *spm, baseval = spmFindBase( spm ); loc2glob = spm->loc2glob; - for( i=0; i<spm->n; i++, loc2glob++ ) { + for( i=0; i<spm->n; i++, loc2glob++ ) + { ig = *loc2glob - baseval; dofi = ( spm->dof > 0 ) ? spm->dof : spm->dofs[ig+1] - spm->dofs[ig]; row = ( spm->dof > 0 ) ? spm->dof * ig : spm->dofs[ig] - baseval; - for( j=0; j<nrhs; j++ ) { - for( k=0; k<dofi; k++ ) { + for( j=0; j<nrhs; j++ ) + { + for( k=0; k<dofi; k++ ) + { rhs[ j * ldb + k ] = bglob[ row + j * ldbglob + k ]; } } @@ -81,4 +84,4 @@ z_spmReduceRHS( const spmatrix_t *spm, (void)b; (void)ldb; #endif -} \ No newline at end of file +}