def plot_scalar_field( field, dof_name, filename, plot_order=1, plot_edges=True, scaling=[-6.0, 7.0, 40.0, 40.0], color_scheme=[(-2.0, [0.2, 0.2, 1.0]), (0.0, [0.8, 0.8, 0.8]), (2.0, [1.0, 0.2, 0.2])], ): ocaml.plot_scalar_field(field, dof_name, filename, color_scheme, plot_order, plot_edges, scaling) return None
def plot_scalar_field(field, dof_name, filename, plot_order=1, plot_edges=True, scaling=[-6.0, 7.0, 40.0, 40.0], color_scheme=[(-2.0, [0.2, 0.2, 1.0]), (0.0, [0.8, 0.8, 0.8]), (2.0, [1.0, 0.2, 0.2])]): ocaml.plot_scalar_field(field, dof_name, filename, color_scheme, plot_order, plot_edges, scaling) return None
[the_field_sigma, field_J], []) return the_field_sigma update_sigma(field_sigma) ocaml.plot_scalar_field( last_field_phi[0], "phi", "/tmp/mesh-phi-initial.ps", # Simple color scheme: # [(-1.2,[0.0,0.0,0.0]),(1.2,[1.0,1.0,1.0])], # More elaborate color scheme - the trick is to include surface lines... [ (-0.5, [0.1, 0.1, 0.4]), (-0.333333, [0.2, 1.0, 0.2]), (-0.166666, [0.1, 0.1, 0.4]), (0.0, [0.2, 1.0, 0.2]), (0.166666, [0.1, 0.1, 0.4]), (0.333333, [0.2, 1.0, 0.2]), (0.5, [0.1, 0.1, 0.4]), ], 2, # plot order 1, # plot edges? [-6.0, 7.0, 40.0, 40.0]) # [(-0.6,[0.1,0.1,0.4]), # (-0.5,[0.2,1.0,0.2]), # (-0.4,[0.1,0.1,0.4]), # (-0.3,[0.2,1.0,0.2]), # (-0.2,[0.1,0.1,0.4]),
recompute_conductivity([h_x, h_y, sigma0, alpha], [the_field_sigma, field_J], []) return the_field_sigma update_sigma(field_sigma) ocaml.plot_scalar_field( last_field_phi[0], "phi", "/tmp/mesh-phi-initial.ps", # Simple color scheme: # [(-1.2,[0.0,0.0,0.0]),(1.2,[1.0,1.0,1.0])], # More elaborate color scheme - the trick is to include surface lines... [ (-0.5, [0.1, 0.1, 0.4]), (-0.333333, [0.2, 1.0, 0.2]), (-0.166666, [0.1, 0.1, 0.4]), (0.0, [0.2, 1.0, 0.2]), (0.166666, [0.1, 0.1, 0.4]), (0.333333, [0.2, 1.0, 0.2]), (0.5, [0.1, 0.1, 0.4]), ], 2, # plot order 1, # plot edges? [-6.0, 7.0, 40.0, 40.0], ) # [(-0.6,[0.1,0.1,0.4]), # (-0.5,[0.2,1.0,0.2]), # (-0.4,[0.1,0.1,0.4]), # (-0.3,[0.2,1.0,0.2]),