R имеетqr() функция, которая выполняет декомпозицию QR с использованием LINPACK или LAPACK (по моему опыту, последняя на 5% быстрее). Основным возвращаемым объектом является матрица "qr", которая содержится в верхней треугольной матрице R (т.е.R=qr[upper.tri(qr)]). Все идет нормально. Нижняя треугольная часть qr содержит Q "в компактной форме". Можно извлечь Q из разложения qr, используяqr.Q(), Я хотел бы найти обратнуюqr.Q(), Другими словами, у меня есть Q и R, и я хотел бы поместить их в объект "qr". R тривиально, а Q нет. Цель состоит в том, чтобы обратиться к немуqr.solve(), что намного быстрее, чемsolve() на больших системах.