inline TridiagonalOperator operator+(const QL_ADD_NAMESPACE(PDE,Identity)<Array<double> >& I, const TridiagonalOperator& D) {
inline TridiagonalOperator operator+(const TridiagonalOperator& D, const QL_ADD_NAMESPACE(PDE,Identity)<Array<double> >& I) {
inline TridiagonalOperator operator+(const QL_ADD_NAMESPACE(PDE,Identity)<Array<double> >& I, const TridiagonalOperator& D) {	Array<double> temp(D.size(),1.0);	TridiagonalOperator result(D.belowDiagonal,D.diagonal+temp,D.aboveDiagonal);	result.setLowerBC(D.theLowerBC);	result.setHigherBC(D.theHigherBC);	return result;}