Biomechanical Joint Model  Author: Anderson Maciel

---

sample_linalg.cpp

Go to the documentation of this file.

// sample code for using amaury's linear algebra library

#include <linalg.h>

using namespace LinAlg;
using namespace std;


#define DOUBLE_PRECISION   // let's use doubles instead of floats
// bounds checking can be activated by using compilation flag BOUNDS_CHECK
#include "linalg.h"

void main()
{
   // vector stuff
   
   VEC A;                           // define an empty vector
   VEC B(dimspace);                 // define a vector with usual size 3 (default val = 0);
   VEC C = B;
   VEC D(dimspace, "1.0 4.5 -10");
   VEC E(dimspace, 0.1234);         // assign a default val
   float array_floats[] = {1., 10., 100.};
   VEC F(array_floats, dimspace);  // note the reverse order when using floats

   cout << "A = " << A << " B = " << B <<" C = " << C <<" D = " << D <<" E = " << E
        << " F = " << F << endl;

   // usual operations
   A = D;                // A is resized as its default size is 0, contents is lost when resizing
   A += A;
   A *= 0.5;
   B = -C + D*10 - E/3;
   

   // array of vectors stuff
   VECArray VA(3, 1, 1, 55);   // define a one-dimensional array of 3 vectors with val -5
   VECArray VB = -VA;          // copy constructor
   VECArray VC(4,4,4,100,2);   // define a 4x4x4 array of vectors of dim 2 with val 100
   VECArray VD(512);           // define an array of 512 vectors of size 3 with val 0
      
   // VECArray is a container so it can be used with std classes
   // copy(VA.begin(), VA.end(), ostream_iterator<REAL>(cout, " "));

   // usual operations are defined
   cout << endl << -VA + 2.5*abs(VB)/6 << endl;
   
   // one-dim array access
   cout << " dot product of " << VA[0] << " and " << VB[0] << " = " 
        << dot_prod(VA[0], VB[0]) << endl;

   // 3-dim array access
   //cout << sqrt(VC(idx3(2,2,2)));
}

---