Beispiel #1
0
 def __init__(self, mu, n_depth_pts, radial_grid):
     self.mu = mu
     self.I_lam = np.zeros(n_depth_pts)
     # outward pointing rays should have diffusion as their initial condition
     if self.mu > 0:
         self.I_lam[0] = planck_fn(1) # TODO: make this the diffusion condition
     # our object is not illuminated from some other source so inward pointing rays have zero intensity to start
     else:
         self.I_lam[0] = 0
     self.tau_grid = np.zeros(n_depth_pts)
Beispiel #2
0
def calc_source_fn(source_fn, epsilon, J_n):
    source_fn[:] = epsilon * planck_fn(1) + (1 - epsilon) * J_n[:]
    return (source_fn)
Beispiel #3
0
def calc_source_fn(source_fn, epsilon, J_n):
    source_fn[:] = epsilon * planck_fn(1) + (1 - epsilon) * J_n[:]
    return (source_fn)