n = FacetNormal(mesh)
r0 = Expression(("0.0"))
b0 = Expression(("1.0", "0.0", "0.0"))
u0 = Expression(("1.0", "0", "0.0"))
Hiptmairtol = 1e-6
HiptmairMatrices = PrecondSetup.MagneticSetup(Magnetic, Lagrange, b0, r0,
Hiptmairtol, params)
B0 = 1.
delta = 0.1
x_on = 4.
x_off = 6.
u0, p0, b0, r0, Laplacian, Advection, gradPres, NScouple, CurlCurl, gradLagr, Mcouple = Generator.ExactSolution(
params, B0, delta, x_on, x_off)
F_NS = -MU * Laplacian + Advection + gradPres - kappa * NScouple
if kappa == 0.0:
F_M = Mu_m * CurlCurl + gradLagr - kappa * Mcouple
else:
F_M = Mu_m * kappa * CurlCurl + gradLagr - kappa * Mcouple
MO.PrintStr("Seting up initial guess matricies", 2, "=", "\n\n", "\n")
u_k, p_k = Generator.Stokes(Velocity, Pressure, F_S, u0, params,
boundaries, domains)
b_k, r_k = Generator.Maxwell(Magnetic, Lagrange, F_M, b0, params,
HiptmairMatrices, Hiptmairtol)
x = Iter.u_prev(u_k, p_k, b_k, r_k)
