def __init__(self, gb, mp, mix, el, kw): self.mp = mp self.mix = mix EllipticModel.__init__(self, gb, **kw) self.el = False if len(el) > 1: self.el = True self.full_grid = gb_full self.el_data = el_data
def __init__(self, gb): EllipticModel.__init__(self, gb, **kw)
def __init__(self, gb, mp=False, kw={}): self.multi_point = mp EllipticModel.__init__(self, gb, **kw)
def __init__(self, gb, mp, mix, kw): self.mp = mp self.mix = mix EllipticModel.__init__(self, gb, **kw)
plt.savefig('figures/' + plot_type + title + '.png') def split_mrst_to_gb(gb, u): Both.el.flux_disc().split(gb, 'pressure', u) if __name__ == "__main__": # Set up grid nc = 4 gb = define_grid(nc, nc) assign_data(gb, FlowData, 'problem') # Copy it for the elimination gb_el, _ = gb.duplicate_without_dimension(0) # Initialize model and solver class problem = EllipticModel(gb) problem_el = EllipticModel(gb_el) # Merge for convenience Both = BothProblems(problem, problem_el) # Define permeability ranges horizontal_permeabilities = [10**i for i in range(-3, 4)] vertical_permeabilities = horizontal_permeabilities.copy() global_errors_porepy = np.zeros( (len(horizontal_permeabilities), len(vertical_permeabilities))) global_errors_mrst = np.zeros( (len(horizontal_permeabilities), len(vertical_permeabilities))) improvement_porepy = np.zeros( (len(horizontal_permeabilities), len(vertical_permeabilities))) improvement_mrst = np.zeros( (len(horizontal_permeabilities), len(vertical_permeabilities)))
def __init__(self, gb, el=''): # Initialize base class kw = {'folder_name': global_folder_name + el} EllipticModel.__init__(self, gb, **kw)