In [1]:

    

import numpy as np
from IPython.display import Image
import h5py
import arrayfire as af

from physics_tests.nonrelativistic_boltzmann.run_cases import run_cases
from physics_tests.nonrelativistic_boltzmann.check_convergence import check_convergence


    

In [2]:

    

@af.broadcast
def tau_collisionless(q1, q2, p1, p2, p3):
    return(np.inf * q1**0 * p1**0)


    

In [3]:

    

run_cases(1, 1, -10, tau_collisionless)
check_convergence()


    

    




---------------------------------------------------------------------------
AssertionError                            Traceback (most recent call last)
<ipython-input-3-b13863e98836> in <module>()
      1 run_cases(1, 1, -10, tau_collisionless)
----> 2 check_convergence()

~/git_repos/Bolt/physics_tests/nonrelativistic_boltzmann/check_convergence.py in check_convergence()
     89     pl.savefig('convergence_plot.png')
     90 
---> 91     assert(abs(poly[0] + 2)<0.25)
     92 

AssertionError: 

In [4]:

    

Image(filename='comparison_N_032.png')


    

    
Out[4]:

In [5]:

    

Image(filename='comparison_N_064.png')


    

    
Out[5]:

In [6]:

    

Image(filename='comparison_N_128.png')


    

    
Out[6]:

In [7]:

    

Image(filename='comparison_N_256.png')


    

    
Out[7]:

In [8]:

    

Image(filename='comparison_N_512.png')


    

    
Out[8]:

In [9]:

    

Image(filename='convergence_plot.png')


    

    
Out[9]:

In [ ]: