FEBio/FEMaterialPoint.h

#pragma once

#include "mat3d.h"
#include "DumpFile.h"
#include <vector>
using namespace std;

//-----------------------------------------------------------------------------


class FEMaterialPoint
{
public:
        FEMaterialPoint(FEMaterialPoint* ppt = 0) : m_pt(ppt){}
        virtual ~FEMaterialPoint() { if (m_pt) { delete m_pt; m_pt = 0; } }

        virtual void Init(bool bflag) = 0;

        virtual FEMaterialPoint* Copy() = 0;

        virtual void Serialize(DumpFile& ar) = 0;

        template <class T>
        T* ExtractData()
        {
                T* p = dynamic_cast<T*>(this);
                if (p) return p; 
                else
                {
                        if (m_pt) return m_pt->ExtractData<T>();
                        else return 0;
                }
        }

protected:
        FEMaterialPoint*        m_pt;   //<! nested point data

public:
        static double time;     // time value
        static double dt; // time increment
};