---

mesh.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_NUM_SEQUENCER_MESH_H__
#define __CSFL_NUM_SEQUENCER_MESH_H__

// Include =====================================================================
#include <csfl/num/sequencer/sequencer.h>
#include <csfl/math/equation/mesh_x.h>
#include <csfl/math/equation/mesh_y.h>
#include <csfl/sys/geom/dline.h>
#include <csfl/sys/geom/dpolygon.h>


// Namespace
namespace csfl {

//==============================================================================
//              enum MeshGenerationType: Mesh Generation Type
//==============================================================================
typedef
        enum _MeshGenerationType {
                mgtKsiInterpolation,
                mgtEtaInterpolation,
                mgtElliptical
        }
        MeshGenerationType;
//==============================================================================
//              Class ISequencerMesh
//              Description: Class for the sequence objects
//==============================================================================

class ISequencerMesh : 
        public ISequencer
{

public:
        ISequencerMesh(
                const IString &_id = IString( "MESH Sequencer") ,
                IContainerEquation *_contEquation =NULL ,
                ISolver *_solver = NULL,
                MeshGenerationType _mot = mgtElliptical
                );
        ISequencerMesh( ISequencerMesh *_s );
        ~ISequencerMesh();


public:

        virtual IEquationMeshX * EquationX() { return eqX; }
        
        virtual IEquationMeshY * EquationY() { return eqY; }
        
        const IString ID() { return id; }

        virtual void Run();

        virtual void SetTolerance( const double & _d ){ tolerance =_d; }

        virtual void EstimateMesh();

        virtual void SetEstimateMesh( bool _b ){ estimatemesh = _b; }

        virtual IMesh * Mesh(){ return mesh ;}

        virtual void InitVariables();

        virtual void ActualizeVariables();

        virtual ILinearSystem * LSys(){ return eqX->LinearSystem(); }
        
        /* O Argumento “_v” define o número máximo de iterações do Solver.*/
        void SetMaximumIteratorSolver( int _v ){ solver->SetIterMax( _v ); }

        /* O Argumento “_mot” define o método de obtenção da malha.*/
        virtual void SetMeshGenerationType( MeshGenerationType _mot ){ mot = _mot; }

        IAction         ActionCalculateX, ActionCalculateY;

        void OutPutMesh( ostream & );

        double MeshError(){ return mesherror; }

        void SetAtraction( const bool _b )
                { isatracted = _b; }

        void SetFixedPoints();
                
        void SetYRelaxation( const double & _d ){ alphay = _d; }
        
        void SetXRelaxation( const double & _d ){ alphax = _d; }



protected:      

        void SetX();
        void SetY();

        IEquationMeshX *eqX;
        IEquationMeshY *eqY;

        MeshGenerationType mot;

        bool estimatemesh, isatracted;
        double tolerance, mesherror;
        double alphax, alphay;

        IMesh *mesh, *meshAnt;

        virtual void Begin();
    virtual void Solution();
        virtual void End(){};
        virtual void Find();
        
        virtual bool SolutionError();

        double VariableError( IMesh *, IMesh * );
};
} // csfl namespace

#endif // __CSFL_NUM_SEQUENCER_MESH_H__

 

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