<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body text="#000000" bgcolor="#FFFFFF">
<div class="moz-cite-prefix">Hi Andrea.<br>
<br>
That test has apparently been removed from dune-fem - at least it
is not there in the master.<br>
So it might be old....<br>
<br>
I had a short look. I seems to me that if you only compile
l2projection then it does not use <br>
petsc. You need <br>
PETSCLINEAROPERATOR == 1<br>
and that is defined if you compile l2projection_petscmatrix And
have usepetsc: true in the parameter file.<br>
<br>
So I don't really see why you are getting a petsc error with
l2projextion since it should not be<br>
using petsc.... But perhaps I missed something.<br>
Have you tried the example in the dune-fem-howto - I am sure that
that worked with the older version<br>
of petsc. <br>
<br>
I tested it myself and I also get an error. In my case an assert
is triggered. Are you defining NDEBUG?<br>
If so then that might explain the problem and I have noticed that
sometimes there when is a segmentation fault <br>
somewhere a petsc error is produced although it has nothing to do
with petsc<br>
So perhaps that is what is happing here.<br>
<br>
The assert btw. is incodimindexset.hh with some invalide index
being used. <br>
That is going to need a bit more investigation I'm afraid. Could
you please also test the example_solver<br>
from the dune-fem-howto - I just don't want to fix some test in
dune-fem which is not used anymore anyway....<br>
<br>
Best<br>
Andreas<br>
<br>
<br>
On 30/07/14 18:50, Sacconi, Andrea wrote:<br>
</div>
<blockquote
cite="mid:EC7DA3C0917A914BAFB620BB60189EBD1463D225@icexch-m2.ic.ac.uk"
type="cite">
<pre wrap="">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 class="moz-txt-link-abbreviated" href="mailto:a.sacconi11@imperial.ac.uk">a.sacconi11@imperial.ac.uk</a>
</pre>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
dune-fem mailing list
<a class="moz-txt-link-abbreviated" href="mailto:dune-fem@dune-project.org">dune-fem@dune-project.org</a>
<a class="moz-txt-link-freetext" href="http://lists.dune-project.org/mailman/listinfo/dune-fem">http://lists.dune-project.org/mailman/listinfo/dune-fem</a>
</pre>
</blockquote>
<br>
</body>
</html>