Hi, I'm using development version of deal r24838. I have adapted step-40 to solve for linear elasticty using the assembly process of step-18. I'm using the actual petsc-dev version.
I run my code using $ mpiexec -n nprocs my_exec -mat_view_matlab > matrix_output I'm observing a very strange behaviour using PetscWrappers::MPI:sparse_matrix. Namely, if I run two times using nprocs =1 I obtain (as expected) the same matrix_output (except for matrix names in matlab_output but it is right). But when using np > 1, matrix values are different: not only between nprocs = 1 and nprocs > 1 runs, but also between runs with same nprocs value greater than one: it seems it prints out the matrix before the PETSc assembly process has terminated (and not all cached values as been transferred into PETSc Mat) . Is it normal? I can provide the test code if you need it. Regards, -- Stefano
_______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
