Comment 2 for bug 891716

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

On 11/17/2011 09:40 PM, Anders Logg wrote:
> Looks like the only remaining thing is periodic BC not working
> with Trilinos. From running the periodic demo:
>
> Error: Unable to add block of values to Epetra matrix. *** Reason:
> Did not manage to perform Epetra_FECrsMatrix::SumIntoGlobalValues.
> *** Where: This error was encountered inside EpetraMatrix.cpp.
>

I wrote a somewhat lengthy email about that, explaining what goes
wrong and asking for opinons... you have even replied to it :)

The bottom line is that there are several related issues:
1. The sparsity pattern does not consider PeriodicBC thus there will
be an error in Trilinos when copying a slave row to a master row since
the matrix is filled / finalized.

2. Dirichlet and PeriodicBC calling and matrix.apply() while applying
bc. Every such call triggers a finalization of the matrix atm.

3. https://bugs.launchpad.net/dolfin/+bug/891452
(which for some reason does not show up in the bug overview, neither
has been sent to the list, strange.)

See the original post
https://bugs.launchpad.net/dolfin/+bug/864741/comments/2

for a more detailed explanation.
The point is that it seems like a more fundamental problem where I
haven't dared to dirty-hack around without opinions from the devs
/maintainers :)

- --
Andre

-----BEGIN PGP SIGNATURE-----
Version: GnuPG v2.0.17 (GNU/Linux)
Comment: Using GnuPG with Mozilla - http://enigmail.mozdev.org/

iQEcBAEBAgAGBQJOxXuoAAoJEA79ggnbq9dm0ToIAJYp09Rkq9midizQnt8zw4gr
197meLuTEuBn7fN4FR9EV3yBqhAwg7Ay/4UW96Da9mt09JoY+CRKEWyzBVhD0TC9
jcDw7aroXZVJMZ9KgHWKiwyBA/7l+4SRXen6DXAneEKIfBJGnNFiul3HTgjBxFdR
ZG7sMELsvynDM7ryt2FpZl89qI34YnOwCNj9vZAXNe8TWA4qsHmc7BTmBWojLYBS
eStuEd73652epXBOmmQHZ47fZg83SxMUJgPkPbDKxgHBaM0sp9Q2X+f1lu5uY/wF
4pIOPIfb8bSA/90Ud3JBIinin43gW4GfL661DzldNlki8kbU5Zzob97wOnCHOpI=
=l8m5
-----END PGP SIGNATURE-----