---

vector_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_VECTOR_FACE_H__
#define __CSFL_FIELD_VALUE_VECTOR_FACE_H__


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


// Namespace
namespace csfl {

//==============================================================================
//              Class IFaceVector
//              Description: Represents a vector in the cell faces
//==============================================================================

class IFaceVector 


{

public:

    static const IFaceVector zero;

    IFaceVector();
    ~IFaceVector();

        inline IFaceVector * operator = ( const IFaceVector & );

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

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

inline IFaceVector * IFaceVector::operator = ( const IFaceVector &_fv )
{

    S->x = _fv.S->x;
    S->y = _fv.S->y;

    W->x = _fv.W->x;
    W->y = _fv.W->y;

    return this;
}

//==============================================================================
//              Class IFieldVectorFace
//              Description: Represents a discretized physical face vector field 
//==============================================================================

class IFieldVectorFace : 
        public IField


{

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

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

    IFieldVectorFace( IFieldVectorFace *_f );

        ~IFieldVectorFace();
 
        inline IVector Vector( const IPoint &_p, TDirection _dir );     
 
        inline IVector Vector( int _i, int _j, TDirection _dir );

        void SetVector( int _i, int _j, TDirection _dir, const IVector &_v );
 
protected:
        IArray2D<IFaceVector> *data;

};

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

inline IVector IFieldVectorFace::Vector( int _i, int _j, TDirection _dir )
{ 
        return Vector( IPoint( _i, _j ), _dir );
}

inline IVector IFieldVectorFace::Vector( const IPoint & _p, TDirection _dir )
{


    if ( _p.x>=0.0 && _p.x<fsize.width && _p.y>=0.0 && _p.y<fsize.height) {
        
        switch (_dir) {
        case (diNorth) :
            return IVector( data->Value(_p).N->x, data->Value(_p).N->y );
        case (diSouth) :
            return IVector( data->Value(_p).S->x, data->Value(_p).S->y );
        case (diEast) :
            return IVector( data->Value(_p).E->x, data->Value(_p).E->y );
        case (diWest) :
            return IVector( data->Value(_p).W->x, data->Value(_p).W->y );
        case (diNULL) :
            return IVector( 0.0, 0.0 );
            
        }
    }
 
    return IVector( 0.0, 0.0 );
    
}

inline void IFieldVectorFace::SetVector( int _i, int _j, TDirection _dir, const IVector & _v )
{

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

     }
}} // csfl namespace

#endif// __CSFL_FIELD_VALUE_VECTOR_FACE_H__

 

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