testGrid.cpp

This example demonstrates the Grid class.

1. Construct two grids in a rectangular region defined by the lower and upper limits that describe the region's body diagonal endpoints. 'grid1' is a 3D Reference grid, while 'grid2' is a 1D Current grid

       Vector3d lowerLimit(1,2,3);
       Vector3d upperLimit(2,4,6);
       int ngrid= 10;
       Grid<Reference> grid1(lowerLimit,upperLimit,ngrid,ngrid,ngrid);
       Grid<Current>   grid2(lowerLimit,upperLimit,ngrid);

2. Write the coordinates of grids to files grid1.grid and grid2.grid

       grid1.write("grid1");
       grid2.write("grid2");

3. An alternate way to construct a grid is to read its coordinates from a file. The following snippet initializes a 3D Current grid and reads its coordinates from a file

       Grid<Current> grid3(ngrid*ngrid*ngrid);
       grid3.read("grid1.grid");

4. We can keep track of the total number of Reference and Current grids using static variables.

       std::cout << "Number of Current grids = " << Grid<Current>::numberOfCurrentGrids << std::endl;
       std::cout << "Number of Reference grids = " << Grid<Reference>::numberOfReferenceGrids << std::endl;

Full code:

/*
 * grid_test.cpp
 *
 *  Created on: Nov 29, 2019
 *      Author: Nikhil
 */
#include "Grid.h"
#include <iostream>
#include "typedef.h"
#include <fstream>
 
int main()
{
    Vector3d lowerLimit(1,2,3);
    Vector3d upperLimit(2,4,6);
    int ngrid= 10;
    Grid<Reference> grid1(lowerLimit,upperLimit,ngrid,ngrid,ngrid);
    Grid<Current>   grid2(lowerLimit,upperLimit,ngrid);
    grid1.write("grid1");
    grid2.write("grid2");
    Grid<Current> grid3(ngrid*ngrid*ngrid);
    grid3.read("grid1.grid");
    std::cout << "Number of Current grids = " << Grid<Current>::numberOfCurrentGrids << std::endl;
    std::cout << "Number of Reference grids = " << Grid<Reference>::numberOfReferenceGrids << std::endl;
}