Python solve_sn Examples

Example #1

File: homogSlab2GMC2011.py Project: ndarmage/RM1D

        'chi': chi,
        'D': D
    }
}
media = [['HM', L]]  # i.e. homogeneously filled

# set b.c.
LBC, RBC = 0, 0

# Initiate inpur data object
Homog2GSlab_data = input_data(xs_media, media, xi, geometry_type, LBC, RBC)
# Initiate solver options object
slvr_opts = solver_options()

# Initiate quadratures object
N = 16
qdata = quad_data(N, L=0)

if __name__ == "__main__":

    import logging as lg
    lg.info("*** Solve the MC 2011 problem ***")
    from snmg1dslab import solve_sn

    flxm, k = solve_sn(Homog2GSlab_data, slvr_opts, qdata)
    basefilen = "../output/kflx_MC2011_LBC%dRBC%d_I%d_L%d_N%d" % (
        Homog2GSlab_data.LBC, Homog2GSlab_data.RBC, Homog2GSlab_data.I,
        Homog2GSlab_data.xi[-1], qdata.N)
    np.save(basefilen + ".npy", np.array([k, flxm]), allow_pickle=True)
    #np.savez(basefilen + ".npz", k=k, flxm=flxm)

Example #2

File: Sood2003_heterSlab1G_P30.py Project: ndarmage/RM1D

sys.path.append('..')  # do not use this at the beginning!
from snmg1dslab import input_data, solver_options, quad_data
import numpy as np

# set b.c.
# LBC, RBC = 0, 0

# or initiate input data object here...
# Heter2GSlab_data = input_data(xs_media, media, xi, geometry_type, LBC, RBC)

# Initiate solver options object
slvr_opts = solver_options()

# Initiate quadratures object
N, L = 16, 0
qdata = quad_data(N, L)

if __name__ == "__main__":

    import logging as lg
    lg.info("*** Solve Sood 2003 problem 30 ***")
    from snmg1dslab import solve_sn

    flxm, k = solve_sn(Heter1GSlab_data, slvr_opts, qdata)

    basefilen = "../output/Sood2003_kflx_SN_LBC%dRBC%d_I%d_N%d_P30" % (
        Heter1GSlab_data.LBC, Heter1GSlab_data.RBC, Heter1GSlab_data.I,
        qdata.N)
    np.save(basefilen + ".npy", np.array([k, flxm]), allow_pickle=True)
    #np.savez(basefilen + ".npz", k=k, flxm=flxm)