def set_mesh(self, mesh):
""" Sets the computational mesh
to be used.
"""
self.mesh = mesh
self.mfd = mfd.MFD()
self.mfd.set_mesh(mesh)
self.mfd.set_m_e_construction_method(6)
import sys
# Include relative path for mimpy library.
import mimpy.mesh.hexmeshwmsfracs as mesh
import mimpy.mfd.mfd as mfd
import numpy as np
res_mfd = mfd.MFD()
res_mfd.set_compute_diagonality(True)
res_mfd.set_m_e_construction_method(0)
#Define the permeability function
def K(p, i, j, k):
return np.eye(3)
#set the mesh and an instance of the HexMesh class
res_mesh = mesh.HexMeshWMSFracs()
#The modification function is applied to the points of the mesh.
#In this case no change is applied.
def mod_function(p, i, j, k):
return p
frac_file = open("fracs.dat")
frac_file.readline()
frac_list = []