[dune-fem] Petsc test fails in parallel case

[Andreas Dedner]

Hi Andrea.

That test has apparently been removed from dune-fem  - at least it is
not there in the master.
So it might be old....

I had a short look. I seems to me that if you only compile l2projection
then it does not use
petsc. You need
PETSCLINEAROPERATOR == 1
and that is defined if you compile l2projection_petscmatrix And have
usepetsc: true in the parameter file.

So I don't really see why you are getting a petsc error with
l2projextion since it should not be
using petsc.... But perhaps I missed something.
Have you tried the example in the dune-fem-howto - I am sure that that
worked with the older version
of petsc.

I tested it myself and I also get an error. In my case an assert is
triggered. Are you defining NDEBUG?
If so then that might explain the problem and I have noticed that
sometimes there when is a segmentation fault
somewhere a petsc error is produced although it has nothing to do with petsc
So perhaps that is what is happing here.

The assert btw. is incodimindexset.hh with some invalide index being used.
That is going to need a bit more investigation I'm afraid. Could you
please also test the example_solver
from the dune-fem-howto - I just don't want to fix some test in dune-fem
which is not used anymore anyway....

Best
Andreas


On 30/07/14 18:50, Sacconi, Andrea wrote:
> Hi all,
>
> I would like to ask you another question about Petsc and its usage in dune-fem.
> Recently, I installed the latest release of this library (3.5.0), and thanks to a patch from Tobias Malkmus, it now compiles correctly with dune-fem-1.4.0.
>
> I ran the demo test contained in dune/fem/petsc/tests/l2projection, and I am sorry to say that it failed in the parallel case.
> You can find attached the output of the parallel run with 4 processors.
>
> Some additional info:
>
> gcc version 4.8.2
> OpenMPI 1.6.5
> Parmetis: 4.0.3
> Petsc 3.5.0
>
> I don't know what I'm doing wrong, since I ran the demo tests in Petsc and they all passed without errors.
> If you need any further information, just let me know.
>
> Cheers,
> Andrea
> __________________________________________________________
>
> Andrea Sacconi
> PhD student, Applied Mathematics
> AMMP Section, Department of Mathematics, Imperial College London,
> London SW7 2AZ, UK
> a.sacconi11 at imperial.ac.uk
>
>
> _______________________________________________
> dune-fem mailing list
> dune-fem at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune-fem

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.dune-project.org/pipermail/dune-fem/attachments/20140730/7b8a42ad/attachment.htm>