---

linearsys9d.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_LINEARSYSTEM9D_H__
#define __CSFL_SYS_MATH_LINEARSYSTEM9D_H__


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


// Namespace
namespace csfl {



class ILinearSystem9D : public ILinearSystem 


{

public:
        
    ILinearSystem9D( int _row, int _col );

        ILinearSystem9D( const ISize &_size );
        ILinearSystem9D( ILinearSystem9D *_lsys );  
        ~ILinearSystem9D();
        
        inline double MatrixA( int _i, int _j );
    
    inline void MatrixA( int _i, int _j, double _val );
     
};


// Inline functions

inline double ILinearSystem9D::MatrixA( int _i, int _j )
{
    //Mudando o formato de A
    int lni = lsize.width;
    int lnj = lsize.height;
    int ni = gsize.width;
    int nj = gsize.height;

    int x = _i;
    int y = _j;
     
    if ( _i>=0 && _i<lni && _j>=0 && _j<lnj ) {
     
    if (x==y)
        return (*matA)( 0, y );

    else if (x==(y+1))
        return (*matA)( 1, y );

    else if (x==(y-1))
        return (*matA)( 5, y );

    else if (x==(y+(ni-1)))
        return (*matA)( 2, y );

    else if (x==(y+ni))
        return (*matA)( 3, y );

    else if (x==(y+(ni+1)))
        return (*matA)( 4, y );

    else if (x==((y-ni)+1))
        return (*matA)( 6, y );

    else if (x==(y-ni))
        return (*matA)( 7, y );

    else if (x==((y-ni)-1))
        return (*matA)( 8, y );

    else 
                return 0.0;
   }
   else
       return 0.0;
}

inline void ILinearSystem9D::MatrixA( int _i, int _j, double _val )
{
    //Mudando o formato de A
    int lni = lsize.width;
    int lnj = lsize.height;
    int ni = gsize.width;
    int nj = gsize.height;

    int x = _i;
    int y = _j;

    if ( _i>=0 && _i<lni && _j>=0 && _j<lnj )
        {
        
        if (x==y)
            (*matA)( 0, y )= _val;
        else
        if (x==(y+1))
            (*matA)( 1, y )= _val;
        else
        if (x==(y-1))
            (*matA)( 5, y )= _val;
        else
        if (x==((y+ni)-1))
            (*matA)( 2, y )= _val;
        else   
        if (x==(y+ni))
            (*matA)( 3, y )= _val;
        else   
        if (x==(y+ni)+1)
            (*matA)( 4, y )= _val;
        else  
        if (x==((y-ni)+1))
            (*matA)( 6, y )= _val;
        else    
        if (x==(y-ni))
            (*matA)( 7, y )= _val;
        else   
        if (x==((y-ni)-1))
            (*matA)( 8, y )= _val;         
    } 
}

} // csfl namespace

#endif // __CSFL_SYS_MATH_LINEARSYSTEM9D_H__

 

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