def avgOplList(opllist, ebds, weight="constant", bound="error"):
getEnergies = lambda n: ebds
def getOpac(n):
rho_n, trad_n = opllist.getDensTemp(n)
return avgopac(opllist.getEnergies(n),
opllist.getOpac(n),
trad_n,
ebds,
weight=weight,
bound=bound)
return OplList(opllist.nopacs, opllist.getDensTemp, getEnergies, getOpac)
def oplListScat(self):
return OplList(self.nopacs, lambda n:
(self.opacs[n].rho, self.opacs[n].tele),
lambda n: self.opacs[n].energies,
lambda n: self.opacs[n].opac3)
def oplRosseland(self):
return OplList(len(self.dens),
lambda n: (self.dens[n], self.temps[n]),
lambda n: self.energies,
lambda n: self.rosseland[n])
def oplPlanck(self):
return OplList(len(self.dens),
lambda n: (self.dens[n], self.temps[n]),
lambda n: self.energies,
lambda n: self.planck[n])
def oplAbsorb(self):
return OplList(
len(self.opacs), lambda n: (self.opacs[n].rho, self.opacs[n].tele),
lambda n: self.opacs[n].energies, lambda n: self.opacs[
n].opac_bb + self.opacs[n].opac_bf + self.opacs[n].opac_ff)
def oplScat(self):
return OplList(len(self.opacs), lambda n:
(self.opacs[n].rho, self.opacs[n].tele),
lambda n: self.opacs[n].energies,
lambda n: self.opacs[n].opac_scat)
def oplTotal(self):
return OplList(len(self.opacs), lambda n:
(self.opacs[n].rho, self.opacs[n].tele),
lambda n: self.opacs[n].energies,
lambda n: self.opacs[n].opac_tot)
