<div dir="ltr"><div class="gmail_extra"><div class="gmail_quote">On Thu, Apr 3, 2014 at 8:43 AM, Åsmund Ervik <span dir="ltr"><<a href="mailto:asmund.ervik@ntnu.no" target="_blank">asmund.ervik@ntnu.no</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">-----BEGIN PGP SIGNED MESSAGE-----<br>
Hash: SHA1<br>
<br>
Dear PETSc users,<br>
<br>
I'm trying to write a Poisson solver in Fortran using the<br>
KSPSetComputeOperators etc. framework. When debugging this, I ended up<br>
modifying ksp/ex22.f so that it matches ksp/ex34.c. The difference<br>
between these is that ex34.c has a non-constant RHS and Neumann BCs,<br>
which is closer to what I want.<br></blockquote><div><br></div><div>If it has Neumann conditions, then it has a null space. Have you included this in</div><div>your solver? That can cause a residual offset.</div><div><br>
</div><div>  Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Now, when I run these two programs I get the following:<br>
<br>
./ex22f_mod<br>
   Residual, L2 norm         1.007312<br>
   Error, sup norm         0.020941<br>
   Error, L1 norm        10.687882<br>
   Error, L2 norm         0.340425<br>
<br>
./ex34<br>
   Residual, L2 norm 1.07124e-05<br>
   Error, sup norm 0.0209405<br>
   Error, L1 norm 0.00618512<br>
   Error, L2 norm 0.000197005<br>
<br>
but when I MatView/VecView the matrix and RHS for each example (write<br>
them to file), there is no difference. I have attached ex22f_mod.F90,<br>
any suggestions on what the error is?<br>
<br>
Once this example works, you can of course include it with PETSc.<br>
<br>
Regards,<br>
Åsmund<br>
<br>
-----BEGIN PGP SIGNATURE-----<br>
Version: GnuPG v2.0.22 (GNU/Linux)<br>
Comment: Using GnuPG with Thunderbird - <a href="http://www.enigmail.net/" target="_blank">http://www.enigmail.net/</a><br>
<br>
iQEcBAEBAgAGBQJTPWWZAAoJED+FDAHgGz19aNYIANyq9k3a5kHzTlAdybkZCYDw<br>
yKtDE5l/4iWYmNL49FH8AocHVcRivLaeJG5CGKqySFtZUXOlC9DM7rn4UmuQecni<br>
dgIQuTk0Ym+OJccHyT5xxnebpFVNrIOTpInfQaDW6dTyeL1svAMeHqslKaGepySL<br>
q/cODbgNDYgl6uumB+POMZevtlM6HPhl/1m7HofcHC9upvTRjSPqP1cg+kg+/8m2<br>
Fdie/X7PCBfShrAys94kNXNcwtbO7taauphkQGMfyl0gUd+lFATG6zrEZdDqSFlV<br>
c44GFFbxW/SRIDBXdOeX9/cy75KW5do1Sildwb6R4H/i7t6/hCJUJuss7FHjmLc=<br>
=tV3N<br>
-----END PGP SIGNATURE-----<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener
</div></div>