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)