Хорошо в реальности код следующий:
_mY m_mY; //исходная матрица проводимостей.

_Builder<_mY> builder;
typedef _LU< _mY > _LU;
_LU::_X B;

Build_Jac(&builder , &Jacobi , &Data);
Build_B(&builder , &Data);

//генерим матрицу
nRes = builder.Build( m_mY ); ATLASSERT( nRes > 0 );
nRes = builder.BuildB( B ); ATLASSERT( nRes > 0 );

std::auto_ptr< _mY > apL ( new _mY );
std::auto_ptr< _mY > apU ( new _mY );
std::auto_ptr< _mY > apmY( &m_mY );

//треангуляция
std::auto_ptr< _LU > apLU( new _LU( *apmY, *apL, *apU ) );

LU разложение реализовано не мной и на входе принимает apmY , вот собственно и вопрос как мне передать этот параметр не создавая в функции переменной m_mY, и если я создам в другой функции m_mY (разделю код) сдается мне что от проблемы двойного удаления это не избавит.