def get_template(self, beta, DRbar):
if (beta, DRbar) in self.cached:
return self.cached[(beta, DRbar)]
if self.interpolate and DRbar <= max(self.DRbars) and DRbar >= min(self.DRbars):
values = map(lambda interp: interp([beta, DRbar])[0], self.interpolators)
# Logarithm of the flux
return numpy.array(values)
else:
_, eps, f_time_int = fast_flux_computation.go(0, 0, beta, DRbar, 1.0, 1.0, False)
# Logarithm of the flux
return numpy.array(numpy.log10(f_time_int))
help="",
required=False,type=float,default=0)
parser.add_argument("-b",help="",
required=False, default=0, type=float)
parser.add_argument("-a",help="Photon index",
required=True, default=2.0, type=float)
parser.add_argument("-d",help="\Delta R / R0", required=True, type=float)
parser.add_argument("-r",help="R0", required=False, type=float, default=1.0)
parser.add_argument("-t",help="tau_star", required=False, type=float, default=1.0)
parser.add_argument("--plot",help="If yes, make a plot with the results",
required=False,choices=['yes','no'],default='no')
if __name__=="__main__":
args = parser.parse_args()
if args.plot.lower()=='yes':
plot = True
else:
plot = False
fast_flux_computation.go(args.m, args.b, args.a, args.d, args.r, args.t, plot)
parser.add_argument("-d", help="\Delta R / R0", required=True, type=float)
parser.add_argument("-r", help="R0", required=False, type=float, default=1.0)
parser.add_argument("-t",
help="tau_star",
required=False,
type=float,
default=1.0)
parser.add_argument("--plot",
help="If yes, make a plot with the results",
required=False,
choices=['yes', 'no'],
default='no')
if __name__ == "__main__":
args = parser.parse_args()
if args.plot.lower() == 'yes':
plot = True
else:
plot = False
fast_flux_computation.go(args.m, args.b, args.a, args.d, args.r, args.t,
plot)
