19int main(
int argc,
char**argv)
29 JParser<> zap(
"Example program to test matrix operations.");
37 catch(
const exception &error) {
38 FATAL(error.what() << endl);
48 for (
int i = 0; i != N; ++i) {
52 for (
int j = 0; j != i; ++j) {
53 A(j,i) = A(i,j) = 0.0;
60 for (
int i = 0; i != numberOfEvents; ++i) {
85 ERROR(
"Matrix product error" << endl);
96 for (
int i = 0; i != numberOfEvents; ++i) {
100 B = A * A * A * A * A;
106 if (!B.isIdentity()) {
107 ERROR(
"Matrix product error" << endl);
int main(int argc, char **argv)
General purpose messaging.
Utility class to parse command line options.
#define make_field(A,...)
macro to convert parameter to JParserTemplateElement object
Auxiliary class for CPU timing and usage.
void print(std::ostream &out, const JScale_t scale=milli_t) const
Print timer data.
Utility class to parse command line options.
This name space includes all other name spaces (except KM3NETDAQ, KM3NET and ANTARES).
JMatrixND & div(const double factor)
Scale matrix.
bool isIdentity(const double eps=std::numeric_limits< double >::min()) const
Test identity.
JMatrixND & mul(const double factor)
Scale matrix.