---

scalar_face.h

Go to the documentation of this file.

/*
 * Copyright (c) 2001 The AUTHORS 
 *   Romeu Andre' Pieritz, Ph.D.        - romeu_pieritz@hotmail.com
 *   Rafael Mendes, Eng.                – mendes_rafael@yahoo.com
 *   Rodrigo Ferraz de Andrade, Eng.    – rferraz@iname.com 
 * All rights reserved.
 *
 * Permission to use, copy and distribute this software and its
 * documentation for educational and personal use, without fee is hereby granted, 
 * provided that the above copyright notice and the following
 * two paragraphs appear in all copies of this software.
 *
 * IN NO EVENT SHALL THE AUTHORS BE LIABLE TO ANY PARTY FOR
 * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, OR CONSEQUENTIAL DAMAGES ARISING OUT
 * OF THE USE OF THIS SOFTWARE AND ITS DOCUMENTATION, EVEN IF THE AUTHORS
 * HAS BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
 *
 * THE AUTHORS SPECIFICALLY DISCLAIMS ANY WARRANTIES,
 * INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY
 * AND FITNESS FOR A PARTICULAR PURPOSE.  THE SOFTWARE PROVIDED HERE UNDER IS
 * ON AN "AS IS" BASIS, AND THE AUTHORS HAVE NO OBLIGATION TO
 * PROVIDE MAINTENANCE, SUPPORT, UPDATES, ENHANCEMENTS, OR MODIFICATIONS.
 *
 * SINMEC Lab. - CFD Sinflow Project - http://www.sinmec.ufsc.br/cfd
 */


// File Define
#ifndef __CSFL_FIELD_VALUE_SCALAR_FACE_H__
#define __CSFL_FIELD_VALUE_SCALAR_FACE_H__


// Include
#include <csfl/field/core/field.h>
#include <csfl/sys/math/matrix.h>


// Namespace
namespace csfl {

//==============================================================================
//              Class IFaceScalar
//              Description: Represents a scalar in the cell faces
//==============================================================================

class IFaceScalar 


{
public:

    IFaceScalar();
    ~IFaceScalar();

        inline IFaceScalar * operator = ( const IFaceScalar & );

    double *N, *S, *E, *W;
};

// =============================================================================
// Inline functions
// =============================================================================

inline IFaceScalar * IFaceScalar::operator = ( const IFaceScalar &_fv )
{
    N = _fv.N;
    S = _fv.S;
    E = _fv.E;
    W = _fv.W;
    
    return this;
}

//==============================================================================
//              Class IFieldScalarFace
//              Description: Represents a discretized physical face scalar field 
//==============================================================================

class IFieldScalarFace : 
        public IField


{
public:
        
        IFieldScalarFace( const IString &_name, int _ni, int _nj );

    IFieldScalarFace( const IString &_name, const ISize &_s );

    IFieldScalarFace( IFieldScalarFace *_f );

        ~IFieldScalarFace();

    inline double Value( const IPoint &_p, TDirection _dir );   

  inline double Value( int _i, int _j, TDirection _dir );

        inline void SetValue( int _i, int _j, TDirection _dir, double _v );
  
protected:
        IArray2D<IFaceScalar> *data;
};



// =============================================================================
// Inline functions
// =============================================================================
        
inline double IFieldScalarFace::Value( int _i, int _j, TDirection _dir )
{
    return Value( IPoint( _i, _j ), _dir );
}

inline double IFieldScalarFace::Value( const IPoint & _p, TDirection _dir )
{

    int i, j;

    i = _p.x;
    j = _p.y;

    if ( i>=0.0 && i<fsize.width && j>=0.0 && j<fsize.height) {
        
        switch (_dir) {
        case (diNorth) :
            return *(*data)(i,j).N;
        case (diSouth) :
            return *(*data)(i,j).S;
        case (diEast) :
            return *(*data)(i,j).E;
        case (diWest) :
            return *(*data)(i,j).W;
        case (diNULL):
            return 0.0;
        }
    }
    return 0.0;
}
        
inline void IFieldScalarFace::SetValue( int _i, int _j, TDirection _dir, double _v )
{  

    if ( _i>=0.0 && _i<fsize.width && _j>=0.0 && _j<fsize.height) {        
        switch (_dir) {
        case (diNorth) :
            *(*data)(_i,_j).N = _v;
            break;
        case (diSouth) :
            *(*data)(_i,_j).S = _v;
            break;
        case (diEast) :
            *(*data)(_i,_j).E = _v;
            break;
        case (diWest) :
            *(*data)(_i,_j).W = _v;
            break;
        }
    }
}    } // csfl namespace

#endif// __CSFL_FIELD_VALUE_SCALAR_FACE_H__

 

CFD Project | Documents | Downloads | Contact us | Use Terms
© SINMEC/EMC/UFSC, 2001.	All rights reserved.	Last Update: Jan. 18, 2002	Webmaster