#!/usr/bin/env python

import csymm
import cfrac

relpath = "reference-data"

nerr = 0

for fn in [ "population-red", "population-blue", "colour-red-blue" ]:
    nerr += csymm.checkMirrorSymmetryY2d("fluid-init-eqdistcylinderfrac-y-2d", fn, relpath)
    nerr += cfrac.checkEqualFrac2d("fluid-init-eqdistcylinderfrac-y-2d", "fluid-init-eqdistcylinder-y-2d", fn, relpath)

if ( nerr > 0 ):
    exit(1)

#!/usr/bin/env python

import cfrac
import csymm

relpath = "reference-data"

nerr = 0

for fn in [ "population-red", "population-blue", "colour-red-blue" ]:
    # Check mirror symmetry
    nerr += csymm.checkMirrorSymmetryY2d("fluid-init-eqdisttwospheresfrac-2d", fn, relpath)
    # Check if fractional is the same as absolute
    nerr += cfrac.checkEqualFrac2d("fluid-init-eqdisttwospheresfrac-2d", "fluid-init-eqdisttwospheres-2d", fn, relpath)

if ( nerr > 0 ):
    exit(1)

#!/usr/bin/env python

import cfrac
import csymm

relpath = "reference-data"

nerr = 0

for fn in ["population-red", "population-blue", "colour-red-blue"]:
    nerr += csymm.checkQuasi1dX2d("fluid-init-eqdistlamellaefrac-y-2d", fn, relpath)
    nerr += cfrac.checkEqualFrac2d("fluid-init-eqdistlamellaefrac-y-2d", "fluid-init-eqdistlamellae-y-2d", fn, relpath)

if nerr > 0:
    exit(1)
#!/usr/bin/env python

import csymm
import cfrac

relpath = "reference-data"

nerr = 0

for fn in [ "population-red", "population-blue", "colour-red-blue" ]:
    nerr += csymm.checkMirrorSymmetryY2d("fluid-init-eqdistcylinder-z-2d", fn, relpath)
    # Cylinder in z-direction should be the same as a 2d sphere
    nerr += cfrac.checkEqualFrac2d("fluid-init-eqdistcylinderfrac-z-2d", "fluid-init-eqdistcylinder-z-2d", fn, relpath)
    nerr += cfrac.checkEqualFrac2d("fluid-init-eqdistcylinderfrac-z-2d", "fluid-init-eqdistspherefrac-2d", fn, relpath)

if ( nerr > 0 ):
    exit(1)