FineMesh.h

Go to the documentation of this file.

/***************************************************************************
                          FineMesh.h  -  description
                             -------------------
    begin                : Wed Jun 13 2001
    copyright            : (C) 2001 by Olivier Roussel & Alexei Tsigulin
    email                : roussel@ict.uni-karlsruhe.de, lpsoft@mail.ru
 ***************************************************************************/

/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/

/***************************************************************************
 * An object FineMesh is a regular fine mesh, used for
 * finite volume computations.
 * It is not used for multiresolution computations.
 *
 * It contains an array of cells <i>*MeshCell</i>.
 *
 * NOTE: for reasons of simplicity, only periodic and Neuman boundary conditions have been
 * implemented.
 */

class FineMesh
{
/*
______________________________________________________________________________________________

PUBLIC METHODS

    Constructor and distructor
______________________________________________________________________________________________

*/
public:

FineMesh();

~FineMesh();
/*
______________________________________________________________________________________________

    Time evolution procedures
______________________________________________________________________________________________

*/

void store();

void storeGrad();

void computeDivergence(int);

void computeDivergence_cell(int);

void RungeKutta_cell(int);

    //each cells
void RungeKutta(int);

void computeIntegral();

void computeCorrection(int);

void computeCorrection_cell(int);

void computeGradient(int);

void computeGradient_cell(int);

void computeTimeAverage();

void checkStability() const;

/*
______________________________________________________________________________________________

    Output procedures
______________________________________________________________________________________________

*/

void writeHeader(const char* FileName) const;

void writeAverage(const char* FileName);

void writeTimeAverage(const char* FileName) const;
/*
______________________________________________________________________________________________

    Backup-restore procedure (to restart a computation)
______________________________________________________________________________________________

*/
void backup();

void restore();


//Parallel
    Cell ***Neighbour_iL,
         ***Neighbour_iU,
         ***Neighbour_jL,
         ***Neighbour_jU,
         ***Neighbour_kL,
         ***Neighbour_kU; 
/*
______________________________________________________________________________________________

PRIVATE METHODS AND VARIABLES
______________________________________________________________________________________________

*/

//parallel modification

/***************************************************************************
 * Returns pointer to the cell located at <i>i, j, k</i>.
 *
 */
public:
  Cell *MeshCell; 
private:
    Cell* cell(const int i, const int j, const int k) const; 
/***************************************************************************
 * Returns pointer to the <i>n</i>-th cell.
 *
 */
    inline Cell* cell(const int n) const; 
/***************************************************************************
 * Compute coarse grid data and store fine grid data in temporary
 *
 */
    void coarsen(); 
/***************************************************************************
 * Restore fine grid data from temporary. Inverse of : coarsen()
 *
 */
    void refine();  
/***************************************************************************
 * Array of cells
 *
 */
    TimeAverageGrid *MyTimeAverageGrid; 
};


/*
______________________________________________________________________________________________

INLINE FUNCTIONS
______________________________________________________________________________________________

*/

inline Cell* FineMesh::cell(const int n) const
{
    return (MeshCell +n);
}