* solveTriangularInPlace(): take a const ref and const_cast it, to allow passing temporary xprs. * improvements, simplifications in LU::solve() * remove remnant of old norm2()