main.cpp

/*
 * ### WIKI 1 ###
 *
 * ## Security enhancement example
 *
 * This example show how to see where an allocation or corruption happen offline and online.
 * Every time an error occur, the library output where the detection happen filename and line,
 *  in order to debug, there is an online option and an offline option
 *
 *  * online: put a breakpoint on the indicated line with your preferred debugger
 *  * offline: set ulimit -c unlimited to activate the core dump file and open the core dump with your debugger
 *
 */

#define SE_CLASS1
#define SE_CLASS2
#define SE_CLASS3
#define THROW_ON_ERROR
#include "Memleak_check.hpp"
#include "data_type/scalar.hpp"
#include "Grid/grid_dist_id.hpp"
#include "Decomposition/CartDecomposition.hpp"
#include "Point_test.hpp"

/*
 * ### WIKI END ###
 */


int main(int argc, char* argv[])
{
	//
	// ### WIKI 2 ###
	//
	// Here we Initialize the library,
	// * message_on_allocation set a message to print when one allocation is reached, the filename and line number can be used to set a breakpoint and analyze the stacktrace.
	// * throw_on_allocation throw when one allocation is reached, producing the termination of the program and a core dump (if no try catch is set-up)
	//
	openfpm_init(&argc,&argv);
	Vcluster & v_cl = create_vcluster();

	throw_on_alloc(10);
	// message_on_alloc(10);

	//
	// ### WIKI 3 ###
	//
	// Create several object needed later, in particular
	// * A 3D box that define the domain
	// * an array of 3 unsigned integer that define the size of the grid on each dimension
	// * A Ghost object that will define the extension of the ghost part for each sub-domain in physical units
	Box<3,float> domain({0.0,0.0,0.0},{1.0,1.0,1.0});
	size_t sz[3];
	sz[0] = 100;
	sz[1] = 100;
	sz[2] = 100;
	
	// Ghost
	Ghost<3,float> g(0.01);
	
	//
	// ### WIKI 4 ###
	//
	// Create a distributed grid in 3D (1° template parameter) defined in R^3 with float precision (2° template parameter)
	// using a CartesianDecomposition strategy (3° parameter) (the parameter 1° and 2° inside CartDecomposition must match 1° and 2°
	// of grid_dist_id)
	//
	// Constructor parameters:
	//
	// * sz: size of the grid on each dimension
	// * domain: where the grid is defined
	// * g: ghost extension
	//
	//
	grid_dist_id<3, float, scalar<float[3]>, CartDecomposition<3,float>> * g_dist = new grid_dist_id<3, float, scalar<float[3]>, CartDecomposition<3,float>>(sz,domain,g);

	//
	// ### WIKI 6 ###
	//
	// print allocated structures
	//

	if (v_cl.getProcessUnitID() == 0)
		print_alloc();

	//
	// ### WIKI 5 ###
	//
	// delete g_dist
	//

	delete g_dist;

	//
	// ### WIKI 6 ###
	//
	// On purpose we try to access a deleted object
	//

	g_dist->getGridInfo();

	//
	// ### WIKI 13 ###
	//
	// Deinitialize the library
	//
	openfpm_finalize();
}