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 = []