29 Mar
2016
29 Mar
'16
5:49 p.m.
On 29 Mar 2016, at 16:09, George Ovchinnikov <lives9@gmail.com> wrote:
a_tilde = (w_tilde * dot(grad(u_tilde), grad(v))) * dx frm_tilde = assemble( a_tilde )
L = f * v * dx rhs = assemble(L)
w_sol = Function(V) solve(frm_tilde, rhs, w_sol, solver_parameters={"ksp_type": "preonly", "pc_type": "lu"})
Without boundary conditions this problem is ill posed. The lay factorisation is failing, presumably due to zero pivot. Lawrence