Mentions légales du service

Skip to content
Snippets Groups Projects
Commit 6a169e28 authored by PICHON Gregoire's avatar PICHON Gregoire
Browse files

Merge default

parents 66e8f85f e5a47e18
No related branches found
No related tags found
No related merge requests found
...@@ -287,7 +287,7 @@ spmReadDriver( pastix_driver_t driver, ...@@ -287,7 +287,7 @@ spmReadDriver( pastix_driver_t driver,
pastix_int_t nnz; pastix_int_t nnz;
if (mpirank == 0) { if (mpirank == 0) {
nnz = spm->colptr[spm->gN] - spm->colptr[0]; nnz = spm->nnz;
} }
MPI_Bcast( spm, 2*sizeof(int)+3*sizeof(pastix_int_t), MPI_CHAR, 0, comm ); MPI_Bcast( spm, 2*sizeof(int)+3*sizeof(pastix_int_t), MPI_CHAR, 0, comm );
......
...@@ -146,7 +146,7 @@ z_spmLaplacian2D( pastix_spm_t *spm, ...@@ -146,7 +146,7 @@ z_spmLaplacian2D( pastix_spm_t *spm,
pastix_int_t i, j, k; pastix_int_t i, j, k;
pastix_int_t nnz = (2*(dim1)-1)*dim2 + (dim2-1)*dim1; pastix_int_t nnz = (2*(dim1)-1)*dim2 + (dim2-1)*dim1;
spm->mtxtype = PastixSymmetric; spm->mtxtype = PastixHermitian;
spm->flttype = PastixComplex64; spm->flttype = PastixComplex64;
spm->fmttype = PastixCSC; spm->fmttype = PastixCSC;
spm->gnnz = nnz; spm->gnnz = nnz;
...@@ -262,7 +262,7 @@ z_spmLaplacian3D( pastix_spm_t *spm, ...@@ -262,7 +262,7 @@ z_spmLaplacian3D( pastix_spm_t *spm,
pastix_int_t i, j, k, l; pastix_int_t i, j, k, l;
pastix_int_t nnz = (2*(dim1)-1)*dim2*dim3 + (dim2-1)*dim1*dim3 + dim2*dim1*(dim3-1); pastix_int_t nnz = (2*(dim1)-1)*dim2*dim3 + (dim2-1)*dim1*dim3 + dim2*dim1*(dim3-1);
spm->mtxtype = PastixSymmetric; spm->mtxtype = PastixHermitian;
spm->flttype = PastixComplex64; spm->flttype = PastixComplex64;
spm->fmttype = PastixCSC; spm->fmttype = PastixCSC;
spm->gnnz = nnz; spm->gnnz = nnz;
......
...@@ -41,15 +41,9 @@ z_spmScal( pastix_spm_t *spm ) ...@@ -41,15 +41,9 @@ z_spmScal( pastix_spm_t *spm )
nnz = spm->nnz; nnz = spm->nnz;
values = spm->values; values = spm->values;
norm = z_spmNorm( PastixFrobeniusNorm, spm ); norm = z_spmNorm( PastixFrobeniusNorm, spm );
printf("NORM IN %.3g\n", norm);
for (i=0; i<nnz; i++){ for (i=0; i<nnz; i++){
values[i] /= norm; values[i] /= norm;
} }
norm = z_spmNorm( PastixFrobeniusNorm, spm );
printf("NORM OUT %.3g\n", norm);
} }
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment