33 static double get(
const double R)
51 static double get(
const double R)
69 static double get(
const double R)
83 int main(
int argc,
char **argv)
93 JParser<> zap(
"Example program to test integration of sphere in any number of dimensions.");
101 catch(
const exception &error) {
102 FATAL(error.what() << endl);
110 const double R = 1.0;
112 typedef JGridPolint1Function1D_t JFunction1D_t;
113 typedef JMultipleMap<
N-2, JPolint1FunctionalGridMap>::typelist JMaplist_t;
114 typedef JMultiFunction<JFunction1D_t, JMaplist_t> JMultiFunction_t;
121 gs.configure(make_multigrid<N-2>(grid));
123 for (JMultiFunction_t::super_iterator i = gs.super_begin(); i != gs.super_end(); ++i) {
125 const double x = (*i).getKey().getLength();
127 for (
int __i = 0; __i != grid.getSize(); ++__i) {
129 const double y = grid.getX(__i);
130 const double z =
R*
R - x*x - y*y;
133 (*i).getValue()[y] = 2.0 * sqrt(z);
135 (*i).getValue()[y] = 0.0;
143 JTimer timer(
"integrator");
151 NOTICE(
"Sphere " <<
N <<
"D" << endl);
157 ASSERT(fabs(U - W) < precision * U);
Utility class to parse command line options.
int main(int argc, char *argv[])
T get(const JHead &header)
Get object from header.
then JShowerPostfit f $INPUT_FILE o $OUTPUT_FILE N
Various implementations of functional maps.
#define ASSERT(A,...)
Assert macro.
I/O formatting auxiliaries.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
static const double PI
Mathematical constants.
then usage $script[distance] fi case set_variable R
General purpose messaging.
Utility class to parse command line options.
Auxiliary data structure for floating point format specification.
int numberOfBins
number of bins for average CDF integral of optical module