---

linearsys.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_SYS_MATH_LINEARSYSTEM_H__
#define __CSFL_SYS_MATH_LINEARSYSTEM_H__


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


// Namespace
namespace csfl {


class ILinearSystem : 
        virtual public ISFLObject 


{

public:
        
        ILinearSystem( int _variables );
        ILinearSystem( int _row, int _col );
        ILinearSystem( const ISize &_size );
        ILinearSystem( ILinearSystem *_lsys );    
        ~ILinearSystem();

        Vector * VectorB() { return vecB; }

        Vector * VectorX() { return vecX; }   
    
        inline void SetVectorB( Vector *_vecB );

        inline void SetVectorX( Vector *_vecX );

    inline virtual double MatrixA( int _i, int _j );

    inline double VectorX( int _i );

    inline double VectorB( int _i );
     
    inline virtual void MatrixA( int _i, int _j, double _val );

    inline void VectorX( int _i, double _val );

    inline void VectorB( int _i, double _val );

        ISize GridSize() { return gsize; }

        ISize Size() { return lsize; }
        
        int Variables() { return lsize.width; }
  
protected:
        ISize gsize, lsize;
        Matrix *matA;
    Vector *vecB, *vecX;
};

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

inline void ILinearSystem::SetVectorB( Vector *_vecB )
{ 
    if (vecB)
        delete vecB;
    vecB = new Vector( _vecB );
}
 
inline void ILinearSystem::SetVectorX( Vector *_vecX )
{
    if (vecX)
        delete vecX;
    vecX = new Vector( _vecX );
}

inline double ILinearSystem::MatrixA( int _i, int _j )
{
    int lni = lsize.width;
    int lnj = lsize.height;
        
    if ( _i>=0 && _i<lni && _j>=0 && _j<lnj ) return (*matA)( _i, _j );
        
        else return 0.0;
}

inline double ILinearSystem::VectorX( int _i )
{
    int ni = vecX->Size();
   
    if ( _i>=0 && _i<ni )
            return (*vecX)( _i );
    else
        return  0.0;
}

inline double ILinearSystem::VectorB( int _i )
{
    int ni = vecB->Size();

    if ( _i>=0 && _i<ni )
            return (*vecB)( _i );
    else
        return 0.0;
}

inline void ILinearSystem::MatrixA( int _i, int _j, double _val )
{    
        int lni = lsize.width;
    int lnj = lsize.height;

    if ( _i>=0 && _i<lni && _j>=0 && _j<lnj ) 
                (*matA)( _i, _j )= _val;
}

inline void ILinearSystem::VectorX( int _i, double _val )
{
    int ni = vecX->Size();

    if ( _i>=0 && _i<ni )
           (*vecX)( _i ) = _val;
}

inline void ILinearSystem::VectorB( int _i, double _val )
{
    int ni = vecB->Size();

    if ( _i>=0 && _i<ni )
           (*vecB)( _i ) = _val;
}

} // csfl namespace

#endif// __CSFL_SYS_MATH_LINEARSYSTEM_H__

 

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