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)