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)