def update_sigma(the_field_sigma): compute_grad_phi = ocaml.prematrix_to_applicator( [], # interface_coeffs [the_field_sigma], "", prematrix_grad_phi) # # laplace_solver = ocaml.laplace_solver( [], # if coeffs [(-1, 1, laplace_dbc)], # dbcs [the_field_sigma], # mid field "", # Petsc name prematrix_laplace) # # result_field_phi = laplace_solver([], laplace_dbc_value, cofield_drho_by_dt) last_field_phi[0] = result_field_phi # # field_J = ocaml.cofield_to_field( [], compute_grad_phi( [], # no target result_field_phi)) # # # Next, let us compute the new condictivity by site-wise operation on # field_J. We just overwrite field_sigma: # recompute_conductivity = ocaml.site_wise_applicator( [mwe_sigma, mwe_J], # all the fields we use [], ["H_x", "H_y", "sigma0", "alpha"], # all the names of extra parameters "", code_recompute_conductivity, True) # # recompute_conductivity([h_x, h_y, sigma0, alpha], [the_field_sigma, field_J], []) return the_field_sigma
def update_sigma(the_field_sigma): compute_grad_phi = ocaml.prematrix_to_applicator([], [the_field_sigma], "", prematrix_grad_phi) # interface_coeffs # # laplace_solver = ocaml.laplace_solver( [], # if coeffs [(-1, 1, laplace_dbc)], # dbcs [the_field_sigma], # mid field "", # Petsc name prematrix_laplace, ) # # result_field_phi = laplace_solver([], laplace_dbc_value, cofield_drho_by_dt) last_field_phi[0] = result_field_phi # # field_J = ocaml.cofield_to_field([], compute_grad_phi([], result_field_phi)) # no target # # # Next, let us compute the new condictivity by site-wise operation on # field_J. We just overwrite field_sigma: # recompute_conductivity = ocaml.site_wise_applicator( [mwe_sigma, mwe_J], # all the fields we use [], ["H_x", "H_y", "sigma0", "alpha"], # all the names of extra parameters "", code_recompute_conductivity, True, ) # # recompute_conductivity([h_x, h_y, sigma0, alpha], [the_field_sigma, field_J], []) return the_field_sigma
def cofield_to_field(cofield, target=None, box_method=False): opt_target = [] if target: opt_target = [target] return ocaml.cofield_to_field(opt_target, cofield, box_method)
def cofield_to_field(cofield, target=None, box_method=False): opt_target = [] if target: opt_target = [target] return ocaml.cofield_to_field(opt_target, cofield, box_method)