stephanmg

bool test_tetgen() {
	/// grid management
	Grid g;
	ug::SubsetHandler sh(g);
        g.attach_to_vertices(ug::aPosition);
        ug::Grid::VertexAttachmentAccessor<ug::APosition3> aaPos(g, ug::aPosition);
        
        /// unit cube
        ug::RegularVertex* temp1 = *g.create<ug::RegularVertex>();
        ug::RegularVertex* temp2 = *g.create<ug::RegularVertex>();

stephanmg

#include <cstdlib>
#include <iostream>
#include <math.h>

using namespace std;

/// matrix dimensions n >= m for my use case
const int n = 3;
const int m = 2;