Thank you for your help with this! After tinkering with it for a few days, I have another question if it isn’t too much of a bother.
I am trying to figure out what these sections are doing.
varf vL (uh, vh) = int2d(Th)(f*vh);
varf vb (uh, vh) = int2d(Th)(1.*vh);
I am under the impression here that f is the non homogeneous term in the original equation. If my equation is homogeneous, then can I set f = 0, or can I just delete this line with vL?
(if so, what do I do with b = vL(0, Vh); and bb = [b, b1] later on?)
Also, what is the vb term representing? Is that what the solution integrates to? For example, if my constraint is that int3d(u) = Constant, would I put that constant in the vb term above or would I put it something below like b1 = Constant?
matrix A = va(Vh, Vh);
real[int] b = vL(0, Vh);
real[int] B = vb(0, Vh);
// Block matrix
matrix AA = [ [ A, B ], [ B', 0 ] ];
set(AA, solver=sparsesolver);
real[int] bb(n+1), xx(n+1), b1(1), l(1);
b1 = 0;
I’m thinking I won’t have to change the block matrix stuff above, except maybe the b1 thing?
// Builds the right hand side block
bb = [b, b1];
// Solve
xx = AA^-1 * bb;
// Set values
[uh[],l] = xx;