anneal=True,
init_mode="random",
alpha=0.0,
)
cmd_kwargs = parse_command_line()
parameters.update(**cmd_kwargs)
if parameters["restart_folder"]:
load_parameters(
parameters, os.path.join(parameters["restart_folder"],
"parameters.dat"))
parameters.update(**cmd_kwargs)
R = parameters["R"]
r = parameters["r"]
res = parameters["res"]
dt = TimeStepSelector(parameters["dt"])
tau = df.Constant(parameters["tau"])
h = df.Constant(parameters["h"])
M = df.Constant(parameters["M"])
geo_map = TorusMap(R, r)
geo_map.initialize(res, restart_folder=parameters["restart_folder"])
W = geo_map.mixed_space((geo_map.ref_el, ) * 4)
# Define trial and test functions
du = df.TrialFunction(W)
chi, xi, eta, etahat = df.TestFunctions(W)
# Define functions
u = df.TrialFunction(W)
verbose=True,
anneal=True,
init_mode="random",
alpha=0.0,
)
cmd_kwargs = parse_command_line()
parameters.update(**cmd_kwargs)
if parameters["restart_folder"]:
load_parameters(parameters, os.path.join(
parameters["restart_folder"], "parameters.dat"))
parameters.update(**cmd_kwargs)
R = parameters["R"]
Rz = parameters["Rz"]
res = parameters["res"]
dt = TimeStepSelector(parameters["dt"])
tau = df.Constant(parameters["tau"])
h = df.Constant(parameters["h"])
M = df.Constant(parameters["M"])
geo_map = EllipsoidMap(R, R, Rz)
geo_map.initialize(res, restart_folder=parameters["restart_folder"])
W = geo_map.mixed_space(4)
# Define trial and test functions
du = df.TrialFunction(W)
chi, xi, eta, etahat = df.TestFunctions(W)
# Define functions
u = df.TrialFunction(W)
verbose=True,
# init_mode="random",
init_mode="nothing",
alpha=0.0,
)
cmd_kwargs = parse_command_line()
parameters.update(**cmd_kwargs)
if parameters["restart_folder"]:
load_parameters(parameters, os.path.join(
parameters["restart_folder"], "parameters.dat"))
parameters.update(**cmd_kwargs)
R = parameters["R"]
r = parameters["r"]
res = parameters["res"]
dt = TimeStepSelector(parameters["dt"])
rho = df.Constant(parameters["rho"])
mu = df.Constant(parameters["mu"])
# geo_map = TorusMap(R, r)
geo_map = SphereMap(R, eps=1e-1)
geo_map.initialize(res, restart_folder=parameters["restart_folder"])
R_el = df.FiniteElement("R", geo_map.ref_mesh.ufl_cell(), 0)
W = geo_map.mixed_space((geo_map.ref_vel, geo_map.ref_el, R_el))
# Define trial and test functions
du = df.TrialFunction(W)
v, q, rt = df.TestFunctions(W)