[완료] 동적할당 부분에서 멈추는데...왜 이런지 모르겟습니다.

using namespace techsoft;
typedef techsoft::matrix<double>	dMatrix;

template <class T>
class matrix
{
public:
   typedef T value_type;
   typedef valarray<T> array_type;
 
   // Constructors
   matrix ();
   matrix (size_t nRow, size_t nCol);
   matrix (size_t nRow, size_t nCol, const T& e);
   matrix (size_t nRow, size_t nCol, const T* Array, MatArrayType aType = C_ARRAY);
   matrix (size_t nRow, size_t nCol, const valarray<T>& vec);
    ~matrix ();
}

 struct base_mat
   {
      valarray<T> val;
      size_t nrow, ncol;
      int	 refcnt;
 
      base_mat (size_t row, size_t col)
            : nrow( row), ncol( col), val( row * col) { refcnt = 1; }
 
      base_mat (const T& v, size_t row, size_t col)
            : nrow( row), ncol( col), val( v, row * col) { refcnt = 1; }
 
      base_mat (const T* v, size_t row, size_t col)
            : nrow( row), ncol( col), val( v, row * col) { refcnt = 1; }
 
      base_mat (size_t row, size_t col, const valarray<T>& v)
            : nrow( row), ncol( col), val( v) { refcnt = 1; }
 
      ~base_mat () {;}
   }

template <class T> inline
matrix<T>::matrix (size_t nRow, size_t nCol)
{
    mPtr = allocator( Allocate, nRow, nCol);
}
 
template <class T>
matrix<T>::matrix (size_t nRow, size_t nCol, const T& e)
{
    mPtr = new base_mat( e, nRow, nCol);
    allocator( Record);
}
 
template <class T>
matrix<T>::matrix (size_t row, size_t col, const T* val, MatArrayType aType)
{
    if (aType == C_ARRAY)
    {
        mPtr = new base_mat( val, row, col); <= 이부분에서 멈춥니다.( base_mat 함수는 위에 있습니다)
        allocator( Record);
    }
    else
    {
        mPtr = allocator( Allocate, row, col);
        T *pv = &mPtr->val[0];
 
        for (size_t k=0, j=0; j < col; j++)
            for (size_t i=0; i < row; i++)
                pv[i*col+j] = val[k++];
    }
}
 
 
template <class T>
matrix<T>::matrix (size_t nRow, size_t nCol, const valarray<T>& vec)
{
   mPtr = new base_mat( nRow, nCol, vec);
   allocator( Record);
}

        double _b[4][3] = {
		{ x.s_,    -x.v_[2],  x.v_[1] },
		{ x.v_[2],  x.s_,    -x.v_[0] },
		{-x.v_[1],  x.v_[0],  x.s_    },
		{-x.v_[0], -x.v_[1], -x.v_[2] }};
 
        dMatrix B(4,3,&_b[0][0]); <=이부분에서 멈춥니다.(하위 소스는 위에 있습니다.)
 
        B = B*0.5*dt;
 
	const double _q[3][3] = { 
		{1, 0, 0 },
		{0, 1, 0 },
		{0, 0, 1},
	};
 
	// Covariance matrix of porcess noises
	const dMatrix Q(3, 3, &_q[0][0]); 
 
	// Predicted state estimate
	x = x * Quaternion (w[0]*dt, w[1]*dt, w[2]*dt);