Example #1
0
 def test_gpu(self):
     
     # just makes sure that CPU and GPU match
     if not GPU: raise SkipTest
     
     gpu_I = gpuscatter.simulate(self.num_molecules, self.q_grid, self.xyzlist, 
                                 self.atomic_numbers, rfloats=self.rfloats, 
                                 poisson_parameter=1.0)
                                 
     cpu_I = cpuscatter.simulate(self.num_molecules, self.q_grid, self.xyzlist, 
                                 self.atomic_numbers, rfloats=self.rfloats, 
                                 poisson_parameter=1.0)
Example #2
0
    def test_gpu_scatter(self):

        if not GPU: raise SkipTest
            
        gpu_I = gpuscatter.simulate(self.num_molecules, self.q_grid, self.xyzlist,
                                    self.atomic_numbers, rfloats=self.rfloats)

        print "GPU", gpu_I
        print "REF", self.ref_I
        
        assert_allclose(gpu_I, self.ref_I, rtol=1e-03,
                        err_msg='scatter: gpu/cpu reference mismatch')
        assert not np.all( gpu_I == 0.0 )
        assert not np.sum( gpu_I == np.nan )