Torquil Macdonald Sørensen <[email protected]> writes: > On 30 January 2014 22:05, Jed Brown <[email protected]> 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.
pgp1mS0E0xnYn.pgp
Description: PGP signature
