[petsc-users] Galerkin multigrid coarsening

[Jed Brown]

Torquil Macdonald Sørensen <torquil at gmail.com> writes:

> On 30 January 2014 22:05, Jed Brown <jed at jedbrown.org> wrote:
>>
>> > I did not have an opportunity to test what happens if I don't set any
>> > restriction, but I'm hoping it will then use P^T, as before.
>>
>> Yes, it will.
>
> Only setting the interpolation, and not the restriction, results in a
> crash for me (with output given at the bottom of this message).
> Looking at the changed code, is it not the case that if "restrct" is
> unset, the code will still try to use it in the MatMatMatMult() inside
> the else-part of the new code? If restrct is not set, it will surely
> differ from "interpolate" in the tests on lines 628 and 641, and
> therefore trigger the else-part of the if-else constructs, thereby
> causing this crash, right?

Sorry, this is my fault.  (I thought restrct was configured before
getting to this bit of code, but clearly not.)  Fixed now.
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: application/pgp-signature
Size: 835 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140131/9cd0b8b9/attachment.pgp>