def test_open_air_diffusivity():
    pn = OpenPNM.Network.Cubic([5, 5, 5], spacing=1)
    pn.add_boundaries()
    Ps = pn.pores('boundary', mode='not')
    Ts = pn.find_neighbor_throats(pores=Ps, mode='intersection', flatten=True)
    geom = OpenPNM.Geometry.Cube_and_Cuboid(network=pn, pores=Ps, throats=Ts)
    geom['pore.diameter'] = 0.999999
    geom['throat.diameter'] = 0.999999
    geom.regenerate(['pore.diameter', 'throat.diameter'], mode='exclude')
    Ps = pn.pores('boundary')
    Ts = pn.find_neighbor_throats(pores=Ps, mode='not_intersection')
    boun = OpenPNM.Geometry.Boundary(network=pn,
                                     pores=Ps,
                                     throats=Ts,
                                     shape='cubes')
    air = OpenPNM.Phases.Air(network=pn)
    Ps = pn.pores()
    Ts = pn.throats()
    phys_air = OpenPNM.Physics.Standard(network=pn,
                                        phase=air,
                                        pores=Ps,
                                        throats=Ts)
    BC1_pores = pn.pores(labels=['top_boundary'])
    BC2_pores = pn.pores(labels=['bottom_boundary'])
    Diff = OpenPNM.Algorithms.FickianDiffusion(network=pn,
                                               phase=air)
    # Assign Dirichlet boundary conditions to top and bottom surface pores
    Diff.set_boundary_conditions(bctype='Dirichlet',
                                 bcvalue=0.6,
                                 pores=BC1_pores)
    Diff.set_boundary_conditions(bctype='Dirichlet',
                                 bcvalue=0.4,
                                 pores=BC2_pores)
    Diff.run()
    Diff.return_results()
    Diff_deff = np.mean(air['pore.diffusivity'])/Diff.calc_eff_diffusivity()
    assert np.round(Diff_deff, 3) == 1
dataset/ETHPy150Open PMEAL/OpenPNM/test/integration/script_test.py/test_open_air_diffusivity/swapped operands of `/`