Ejemplo n.º 1
0
 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))
Ejemplo n.º 2
0
 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))
Ejemplo n.º 3
0
 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)
Ejemplo n.º 4
0
 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)