Math operations in Chrono (demo_CH_math.cpp)

Tutorial on main mathematical tools:

• vectors
• quaternions
• matrices

No GUI: only text output.

// =============================================================================
// PROJECT CHRONO - http://projectchrono.org
//
//
// Use of this source code is governed by a BSD-style license that can be found
// in the LICENSE file at the top level of the distribution and at
//
// =============================================================================
// Authors: Alessandro Tasora, Radu Serban
// =============================================================================
//
// Demo on how to use Chrono mathematical functions
//
// =============================================================================
#include "chrono/core/ChMatrix.h"
#include "chrono/core/ChLog.h"
#include "chrono/core/ChVector.h"
#include "chrono/core/ChException.h"
#include "chrono/core/ChMathematics.h"
using namespace chrono;
int main(int argc, char* argv[]) {
GetLog() << "Copyright (c) 2017 projectchrono.org\nChrono version: " << CHRONO_VERSION << "\n\n";
// To write something to the console, use the chrono::GetLog()
// statement, which returns a global output stream to the console (just
// like the std::out stream).
GetLog() << "\n=== Computing integrals of functions in 1D/2D/3D ===\n\n";
// Define a y=f(x) function by inheriting ChIntegrable1D:
class MySine1d : public ChIntegrable1D<double> {
public:
void Evaluate(double& result, const double x) { result = sin(x); }
};
// Create an object from the function class
MySine1d mfx;
// Invoke 6th order Gauss-Legendre quadrature on 0..PI interval:
double qresult = 0;
GetLog() << "Quadrature 1d result: " << qresult << " (analytic solution: 2.0) \n";
// Other quadrature tests, this time in 2D
class MySine2d : public ChIntegrable2D<double> {
public:
void Evaluate(double& result, const double x, const double y) { result = sin(x); }
};
MySine2d mfx2d;
qresult = 0;
ChQuadrature::Integrate2D<double>(qresult, mfx2d, 0, CH_C_PI, -1, 1, 6);
GetLog() << "Quadrature 2d result: " << qresult << " (analytic solution: 4.0) \n";
// Other quadrature tests, this time with vector function (that is, integrates 2x1 matrix)
class MySine2dM : public ChIntegrable2D<ChMatrixNM<double, 1, 2>> {
public:
void Evaluate(ChMatrixNM<double, 1, 2>& result, const double x, const double y) {
result(0, 0) = x * y;
result(0, 1) = 0.5 * y * y;
}
};
MySine2dM mfx2dM;
resultM.setZero();
ChQuadrature::Integrate2D<ChMatrixNM<double, 1, 2>>(resultM, mfx2dM, 0, 1, 0, 3, 6);
GetLog() << "Quadrature 2d matrix result: " << resultM << " (analytic solution: 2.25, 4.5) \n";
return 0;
}
Base class for 1D integrand T=f(x) to be used in ChQuadrature.