fetisolver.h

Go to the documentation of this file.

/*
 *
 *                 #####    #####   ######  ######  ###   ###
 *               ##   ##  ##   ##  ##      ##      ## ### ##
 *              ##   ##  ##   ##  ####    ####    ##  #  ##
 *             ##   ##  ##   ##  ##      ##      ##     ##
 *            ##   ##  ##   ##  ##      ##      ##     ##
 *            #####    #####   ##      ######  ##     ##
 *
 *
 *             OOFEM : Object Oriented Finite Element Code
 *
 *               Copyright (C) 1993 - 2013   Borek Patzak
 *
 *
 *
 *       Czech Technical University, Faculty of Civil Engineering,
 *   Department of Structural Mechanics, 166 29 Prague, Czech Republic
 *
 *  This library is free software; you can redistribute it and/or
 *  modify it under the terms of the GNU Lesser General Public
 *  License as published by the Free Software Foundation; either
 *  version 2.1 of the License, or (at your option) any later version.
 *
 *  This program is distributed in the hope that it will be useful,
 *  but WITHOUT ANY WARRANTY; without even the implied warranty of
 *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 *  Lesser General Public License for more details.
 *
 *  You should have received a copy of the GNU Lesser General Public
 *  License along with this library; if not, write to the Free Software
 *  Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301  USA
 */

#ifndef fetisolver_h
#define fetisolver_h

#include "../sm/FETISolver/feticommunicator.h"
#include "sparselinsystemnm.h"
#include "sparsemtrx.h"
#include "floatarray.h"
#include "floatmatrix.h"
#include "processcomm.h"


#define _IFT_FETISolver_Name "feti"
#define _IFT_FETISolver_maxiter "maxiter"
#define _IFT_FETISolver_maxerr "maxerr"
#define _IFT_FETISolver_limit "limit"
#define _IFT_FETISolver_energynormflag "energynormflag"


namespace oofem {
class Domain;
class EngngModel;


 #define FETISOLVER_MAX_RBM 6
 #define FETISOLVER_ZERONUM 1.e-40

class FETISolver : public SparseLinearSystemNM
{
private:
    int nse;
    int ni;
    double err;
    double limit;
    FloatMatrix rbm;
    FloatMatrix l;
    IntArray se;
    IntArray rbmAddr;
    FloatArray w;
    //
    FloatArray qq, q, dd, g, d, p, pp, gamma, localGammas;
    //
    IntArray nsem;
    ProcessCommunicatorBuff pcbuff;
    ProcessCommunicator processCommunicator;
    CommunicatorBuff *commBuff;
    FETICommunicator *masterCommunicator;
    IntArray masterCommMap;
    int energyNorm_comput_flag;
public:
    FETISolver(Domain * d, EngngModel * m);
    virtual ~FETISolver();

    virtual NM_Status solve(SparseMtrx &A, FloatArray &b, FloatArray &x);

    int estimateMaxPackSize(IntArray &, DataStream &, int &);
    void setUpCommunicationMaps();

    virtual IRResultType initializeFrom(InputRecord *ir);

    // identification
    virtual const char *giveClassName() const { return "FETISolver"; }
    virtual LinSystSolverType giveLinSystSolverType() const { return ST_Feti; }

    void projection(FloatArray &v, FloatMatrix &l, FloatMatrix &l1);

    int packRBM(ProcessCommunicator &processComm);
    int masterUnpackRBM(ProcessCommunicator &processComm);
    int packQQProducts(ProcessCommunicator &processComm);
    int masterUnpackQQProduct(ProcessCommunicator &processComm);
    int packSolution(ProcessCommunicator &processComm);
    int unpackSolution(ProcessCommunicator &processComm);
    int packResiduals(ProcessCommunicator &processComm);
    int unpackResiduals(ProcessCommunicator &processComm);
    int packDirectionVector(ProcessCommunicator &processComm);
    int unpackDirectionVector(ProcessCommunicator &processComm);
    int packPPVector(ProcessCommunicator &processComm);
    int unpackPPVector(ProcessCommunicator &processComm);
    int packGammas(ProcessCommunicator &processComm);
    int unpackGammas(ProcessCommunicator &processComm);

    int masterMapRBM();
    int masterMapQQProduct();
    int masterMapSolution();
    int masterMapResiduals();
    int masterMapDirectionVector();
    int masterMapPPVector();
    int masterMapGammas();

    enum { FETISolverZeroTag, NumberOfRBMMsg, RBMMessage, QQMessage, SolutionMessage, ResidualMessage, DirectionVectorMessage, PPVectorMessage, GammasMessage, FETISolverIterationContinue, FETISolverIterationBreak };
    virtual SparseMtrxType giveRecommendedMatrix(bool symmetric) const { return symmetric ? SMT_Skyline : SMT_SkylineU; }
};
} // end namespace oofem

#endif // fetisolver_h