@Test public void updateHessian() { DMatrixRMaj J = RandomMatrices_DDRM.rectangle(30,6,rand); DMatrixRMaj H = new DMatrixRMaj(1,1); CommonOps_DDRM.multTransA(J,J,H); ((HessianLeastSquares)alg).updateHessian(convert(J)); DMatrixRMaj expected = new DMatrixRMaj(6,1); DMatrixRMaj found = expected.createLike(); alg.extractDiagonals(found); CommonOps_DDRM.extractDiag(H,expected); assertTrue(MatrixFeatures_DDRM.isIdentical(expected,found,UtilEjml.TEST_F64)); } }