DOLFIN

Bug #891716
Comment #4

Comment 4 for bug 891716

Revision history for this message

Garth Wells (garth-wells) wrote on 2011-11-17: Re: [Bug 891716] Re: Unit tests fail to run if Trilinos is installed but not PETSc

On 17 November 2011 21:37, Anders Logg <email address hidden> wrote:
> On Thu, Nov 17, 2011 at 10:24:56PM +0100, Andre Massing wrote:
>> -----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 :)
>
> I know and I've read it carefully!
>
> I think I responded and said that this needs to be addressed by
> incorporating boundary conditions into the FunctionSpace by
> implementing constrained spaces. I even linked to a blueprint.
>

A blueprint on this has been about for coming up two years:

https://blueprints.launchpad.net/dolfin/+spec/bcs-in-sparsitypattern

Garth

> Bottom line: we won't fix this now so I'll add a check in that demo
> that it doesn't work with Epetra.
>
> --
> Anders
>
>
>> 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-----
>
> --
> You received this bug notification because you are a member of DOLFIN
> Core Team, which is subscribed to DOLFIN.
> https://bugs.launchpad.net/bugs/891716
>
> Title:
> Unit tests fail to run if Trilinos is installed but not PETSc
>
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/dolfin/+bug/891716/+subscriptions
>

On 17 November 2011 21:37, Anders Logg <logg@simula.no> wrote:
> On Thu, Nov 17, 2011 at 10:24:56PM +0100, Andre Massing wrote:
>> -----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 :)
>
> I know and I've read it carefully!
>
> I think I responded and said that this needs to be addressed by
> incorporating boundary conditions into the FunctionSpace by
> implementing constrained spaces. I even linked to a blueprint.
>

A blueprint on this has been about for coming up two years:

https://blueprints.launchpad.net/dolfin/+spec/bcs-in-sparsitypattern

Garth

> Bottom line: we won't fix this now so I'll add a check in that demo
> that it doesn't work with Epetra.
>
> --
> Anders
>
>
>> 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-----
>
> --
> You received this bug notification because you are a member of DOLFIN
> Core Team, which is subscribed to DOLFIN.
> https://bugs.launchpad.net/bugs/891716
>
> Title:
>  Unit tests fail to run if Trilinos is installed but not PETSc
>
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/dolfin/+bug/891716/+subscriptions
>