Hello Anna, under the Inequality constraint folder, I'm running "python buoy-swe.py" and, besides the TSFC/COFFEE prints, I'm obtaining: mu_bar ------ [ 0.34693095 0.28326793 0.20030068 0. -0. ] Matrix b_mu ----------- [[ 1.37331778 0.6536432 0. 0. 0. ] [ 0.6536432 2.33042494 0.50155807 0. 0.] [ 0. 0.50155807 1.54031885 0.20775216 0.] [ 0. 0. 0.20775216 0.20775216 0. ] [ 0. 0. 0. 0. 0. ]] What values do you expect for b_mu ? Am I running the right thing? I tried switching off some optimisations but I still get the same output. Thanks! -- Fabio 2016-06-08 15:16 GMT+02:00 Anna Kalogirou <a.kalogirou@leeds.ac.uk>:
Hi Fabio,
Sorry for that. It should be public now.
Anna.
On 08/06/16 13:56, Fabio Luporini wrote:
Hello Anna,
the code is not accessible as it appears to be in a private repository. Could you fix this? Thanks
-- Fabio
2016-06-08 14:10 GMT+02:00 Anna Kalogirou <a.kalogirou@leeds.ac.uk>:
Dear all,
I have a simple question regarding a matrix assembly in Firedrake which I believe is not computed correctly.
The relevant code can be found here <https://bitbucket.org/annakalog/buoy2d/src/14334d3c20b9f10ed7c1246cde9e3cb60b1c75e4/Inequality%20constraint/?at=master> and the issue is with the B_mu matrix in solvers.py, line 51. I compare the result with my own Matlab code and the two don't agree. Also, the rest of the matrices M, A, C_lam agree in both sets of code.
Can someone help with what the problem might be?
Thanks,
Anna.
--
Dr Anna Kalogirou Research Fellow School of Mathematics University of Leeds
_______________________________________________ firedrake mailing listfiredrake@imperial.ac.ukhttps://mailman.ic.ac.uk/mailman/listinfo/firedrake