primaryfield.h

Go to the documentation of this file.

/*
 *
 *                 #####    #####   ######  ######  ###   ###
 *               ##   ##  ##   ##  ##      ##      ## ### ##
 *              ##   ##  ##   ##  ####    ####    ##  #  ##
 *             ##   ##  ##   ##  ##      ##      ##     ##
 *            ##   ##  ##   ##  ##      ##      ##     ##
 *            #####    #####   ##      ######  ##     ##
 *
 *
 *             OOFEM : Object Oriented Finite Element Code
 *
 *               Copyright (C) 1993 - 2025   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 primaryfield_h
#define primaryfield_h
 
#include "field.h"
#include "interface.h"
#include "floatarray.h"
#include "valuemodetype.h"
#include "contextioresulttype.h"
#include "contextmode.h"
#include "timestep.h"
 
#include <vector>
 
namespace oofem {
class PrimaryField;
class Dof;
class BoundaryCondition;
class InitialCondition;
class UnknownNumberingScheme;

class OOFEM_EXPORT EIPrimaryFieldInterface : public Interface
{
public:
 
    EIPrimaryFieldInterface() : Interface() { }

    virtual int EIPrimaryFieldI_evaluateFieldVectorAt(FloatArray &answer, PrimaryField &pf,
                                                      const FloatArray &coords, IntArray &dofId, ValueModeType mode, TimeStep *tStep) = 0;
};
 

class OOFEM_EXPORT PrimaryField : public Field
{
protected:
    int actualStepNumber;
    int actualStepIndx;
    int nHistVectors;
    std :: vector< FloatArray >solutionVectors;
    std :: vector< FloatArray >prescribedVectors;
    std :: vector< TimeStep >solStepList;
    EngngModel *emodel;
    int domainIndx;
 
public:
    PrimaryField(EngngModel * a, int idomain, FieldType ft, int nHist);
    virtual ~PrimaryField();

    virtual void initialize(ValueModeType mode, TimeStep *tStep, FloatArray &answer, const UnknownNumberingScheme &s);
 
    // These functions are hardcoded to assume the default numbering scheme (as is the rest of primaryfield)
    void storeDofManager(TimeStep *tStep, DofManager &dman);
    void storeInDofDictionaries(TimeStep *tStep);
    void readDofManager(TimeStep *tStep, DofManager &dman);
    void readFromDofDictionaries(TimeStep *tStep);

    virtual void applyDefaultInitialCondition();
    void applyInitialCondition(InitialCondition &ic);
    virtual void applyBoundaryCondition(TimeStep *tStep);
    void applyBoundaryCondition(BoundaryCondition &bc, TimeStep *tStep);

    virtual double giveUnknownValue(Dof *dof, ValueModeType mode, TimeStep *tStep);
 
    int evaluateAt(FloatArray &answer, const FloatArray &coords, ValueModeType mode, TimeStep *tStep) override;
    int evaluateAt(FloatArray &answer, DofManager *dman, ValueModeType mode, TimeStep *tStep) override;

    virtual int __evaluateAt(FloatArray &answer, DofManager *dman,
                             ValueModeType mode, TimeStep *tStep, IntArray *dofId);
    virtual int __evaluateAt(FloatArray &answer, const FloatArray &coords,
                             ValueModeType mode, TimeStep *tStep, IntArray *dofId);
    virtual FloatArray *giveSolutionVector(TimeStep *tStep);

    virtual void update(ValueModeType mode, TimeStep *tStep, const FloatArray &vectorToStore, const UnknownNumberingScheme &s);

    virtual void advanceSolution(TimeStep *tStep);
 
    void saveContext(DataStream &stream) override;
    void restoreContext(DataStream &stream) override;
 
    const char *giveClassName() const override { return "PrimaryField"; }
 
    int giveActualStepNumber() { return actualStepNumber; }
protected:
    int resolveIndx(TimeStep *tStep, int shift);
    FloatArray *giveSolutionVector(int);
    FloatArray *givePrescribedVector(int);
};
} // end namespace oofem
#endif // primaryfield_h