PrintGrid.h

Go to the documentation of this file.

/***************************************************************************
                          PrintGrid.h  -  description
                             -------------------
    begin                : Tue Mar 19 2002
    copyright            : (C) 2002 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.                                   *
 *                                                                         *
 ***************************************************************************/

class PrintGrid
{
/*
______________________________________________________________________________________________

DESCRIPTION

    This is a grid of level L used to print data into file.

______________________________________________________________________________________________

PUBLIC METHODS

    Constructor and distructor
______________________________________________________________________________________________

*/

public:
    PrintGrid(int L); 
    ~PrintGrid();   
/*
______________________________________________________________________________________________

    Set and get procedures
______________________________________________________________________________________________

*/

void   setValue(const int i, const int j, const int k, const Vector& UserAverage);

Vector value(const int i, const int j, const int k) const;

real value(const int i, const int j, const int k, const int QuantityNo) const;

real cellValue(const int i, const int j, const int k, const int QuantityNo) const;

inline real density(const int i, const int j, const int k) const;

inline real psi(const int i, const int j, const int k) const;

real pressure(const int i, const int j, const int k) const;

real temperature(const int i, const int j, const int k) const;

real concentration(const int i, const int j, const int k) const;

inline real energy(const int i, const int j, const int k) const;

inline real velocity(const int i, const int j, const int k, const int AxisNo) const;

Vector  velocity(const int i, const int j, const int k) const;

real divergenceB(const int i, const int j, const int k) const;

real vorticity(const int i, const int j, const int k) const;

inline real magField(const int i, const int j, const int k, const int AxisNo) const;

Vector  magField(const int i, const int j, const int k) const;

/*
______________________________________________________________________________________________

    Refresh and predict procedure
______________________________________________________________________________________________

*/

void refresh();

void predict(const int l, const int i, const int j, const int k);

void computePointValue();
/*
______________________________________________________________________________________________

PRIVATE VARIABLES AND METHODS
______________________________________________________________________________________________

*/

private:
  Vector tempValue(const int l, const int i, const int j, const int k) const;

    int localScaleNb;   
    int elementNb;      
    Vector *Q;      
    Vector *Qt;     
};


/*
______________________________________________________________________________________________

INLINE FUNCTIONS
______________________________________________________________________________________________

*/

inline real PrintGrid::density(const int i, const int j, const int k) const
{
    return value(i,j,k,1);
}

/*
______________________________________________________________________________________________

*/

inline real PrintGrid::psi(const int i, const int j, const int k) const
{
    return value(i,j,k,6);
}

/*
______________________________________________________________________________________________

*/

inline real PrintGrid::energy(const int i, const int j, const int k) const
{
    return value(i,j,k,5);
}


/*
______________________________________________________________________________________________

*/

inline real PrintGrid::velocity(const int i, const int j, const int k, const int AxisNo) const
{
    return cellValue(i,j,k,AxisNo+1)/cellValue(i,j,k,1);
}

/*
______________________________________________________________________________________________

*/

inline real PrintGrid::magField(const int i, const int j, const int k, const int AxisNo) const
{
    return cellValue(i,j,k,AxisNo+6);
}