#ifndef _PHYSICS_TXX
#define _PHYSICS_TXX

#include "Physics.h"

namespace mial
{
	// constructor:
	template<class Type, int nDims, class MType, class VType> 
		Physics<Type, nDims, MType, VType >
		::Physics()
	{
		numDeformations = 0;
		time =0;
	}

	// destructor:
	template<class Type, int nDims, class MType, class VType>
		Physics<Type,nDims,MType, VType>
		::~Physics()
	{
		
	}
	template<class Type, int nDims, class MType, class VType>
		bool 
		Physics<Type,nDims,MType, VType>
		::addDeformation(Deformation<Type,nDims>* b)
	{
		if(b!=NULL)
		{
			typename Deformation<Type,nDims>::Pointer dp = b;
			deformationsList.push_back(dp); //Copy by value
			numDeformations++;
			return true;
		}
		else
		{
			Error e;
			e.deformationNumber = numDeformations;
			e.msg = "Cannot add null deformation.";
			return false;
		}
	}

} // end namespace mial

#endif