def generate_fast_HTI_PDE_solution(self, domain):
pde = WavePDE(domain, [("c11", self.c11),
("c23", self.c23), ("c13", self.c13), ("c33", self.c33),
("c44", self.c44), ("c66", self.c66)])
pde.getSolverOptions().setSolverMethod(SolverOptions.HRZ_LUMPING)
pde.setSymmetryOn()
dim = pde.getDim()
X = domain.getX()
y = Vector([2.,3.,4.][:dim], DiracDeltaFunctions(domain))
du = grad(X*X)
D = 2500.*kronecker(dim)
pde.setValue(D=D, y_dirac=y, du=du)
return pde.getSolution()
def generate_fast_HTI_PDE_solution(self, domain):
pde = WavePDE(domain, [("c11", self.c11),
("c23", self.c23), ("c13", self.c13), ("c33", self.c33),
("c44", self.c44), ("c66", self.c66)])
pde.getSolverOptions().setSolverMethod(SolverOptions.HRZ_LUMPING)
pde.setSymmetryOn()
dim = pde.getDim()
X = domain.getX()
y = Vector([2.,3.,4.][:dim], DiracDeltaFunctions(domain))
du = grad(X*X)
D = 2500.*kronecker(dim)
pde.setValue(D=D, y_dirac=y, du=du)
return pde.getSolution()