def test_get_LHS_matrix(self):
self.DUT.set_frequency(self.frequency)
self.DUT.init_problem()
# test data generated using a suspected-working version on 31 May 2011
actual_LHSmat = self.DUT.get_LHS_matrix().todense()
desired_file = Paths.get_module_path_filename('LHS_matrix.npy', __file__)
desired_LHSmat = N.load(desired_file)
self.assertTrue(N.allclose(
actual_LHSmat, desired_LHSmat, rtol=1e-10, atol=3e-15))
def test_get_LHS_matrix(self):
self.DUT.set_frequency(self.frequency)
self.DUT.init_problem()
# test data generated using a suspected-working version on 31 May 2011
actual_LHSmat = self.DUT.get_LHS_matrix().todense()
desired_file = Paths.get_module_path_filename('LHS_matrix.npy',
__file__)
desired_LHSmat = N.load(desired_file)
self.assertTrue(
N.allclose(actual_LHSmat, desired_LHSmat, rtol=1e-10, atol=3e-15))
def setUp(self):
data_dir = Paths.get_module_path_filename(self.data_dir, __file__)
data_file = open(os.path.join(data_dir, self.test_data_file))
data = self.data = pickle.load(data_file)
self.mesh = dolfin.Mesh(os.path.join(data_dir, data['meshfile']))
self.discretisation_order = data['order']
self.function_space = dolfin.FunctionSpace(
self.mesh, "Nedelec 1st kind H(curl)", self.discretisation_order)
self.k0 = 2*np.pi*data['freq']/c0
self.DUT = power_flux.VariationalSurfaceFlux(self.function_space)
def setUp(self):
data_dir = Paths.get_module_path_filename(self.data_dir, __file__)
data_file = open(os.path.join(data_dir, self.test_data_file))
data = self.data = pickle.load(data_file)
self.mesh = dolfin.Mesh(os.path.join(data_dir, data['meshfile']))
self.discretisation_order = data['order']
self.function_space = dolfin.FunctionSpace(self.mesh,
"Nedelec 1st kind H(curl)",
self.discretisation_order)
self.k0 = 2 * np.pi * data['freq'] / c0
self.DUT = power_flux.VariationalSurfaceFlux(self.function_space)
