main.cpp
No OneTemporary
Actions

Subscribers

None

File Metadata

	#include <iostream>
	#include <sstream>

	#include <mpi.h>

	#include "simulation.h"


	// Define the constants
	#define G 127267.200000000 // Gravity, 9.82(3.6)^21000 in [km / hr^2]
	#define SIZE 500 // Size of map, Size*Size [km]
	#define NX 2001 // Number of cells in each direction on the grid
	#define TEND 0.20 // Simulation time in hours [hr]
	#define DIR "../MatlabSource/"
	#define OUTPUT "../MatlabSource/OUT.bin"


	int main (int argc, char *argv[]) {

	// initialize MPI
	MPI_Init(&argc, &argv);
	int prank, psize;

	// Usefull to debug
	MPI_Comm_rank(MPI_COMM_WORLD, &prank);
	MPI_Comm_size(MPI_COMM_WORLD, &psize);

	printf("I’m process %d out of %d\n", prank+1, psize);

	// compute the constant
	double DX = SIZE*1.0/NX;

	//std::cout << G << " " << SIZE << " " << NX << " " << TEND << " dx: "<< DX << std::endl;

	// Construct the simulation
	Simulation simu(NX,G,SIZE,TEND,DX, DIR,OUTPUT, MPI_COMM_WORLD);

	// Compute the simulation
	float totalTime;
	int k;
	std::tie(totalTime, k) = simu.compute();

	std::cout << "Total time : " << totalTime << " Number of steps : " << k << std::endl;


	// Close MPI
	MPI_Finalize();
	}