Amoiensis / Matrix_hub

A lib of Matrix operation for C language. (矩阵运算库--C语言)

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

求特征值时的QR分解M_QR以及转上hessenberg矩阵函数M_householder有误

645770225 opened this issue · comments

您好,感谢作者开源
在使用时,我发现特征值特征向量计算时,出现了和[schumipush]一样的问题,检查发现,是作者做householder变换时,没有把原矩阵乘以当前循环得到的变换矩阵再传入下一次循环,而是直接把原矩阵直接传入下一次循环,因此我修改如下:
1613行

Matrix* x = M_Cut(Ri,i+1,_END_,i,i);//此处"_mat"修改为Ri

1657行

x = M_Cut(Ri,i+1,_END_,i+1,i+1);//此处"_mat"修改为Ri

你好,感谢你对本项目的关注和使用!
确实如您发现的,在当时实现的时候误写,导致了这个QR分解、特征计算部分出现了计算问题。这个问题,将在下一版本中修正,并将参考您提出的这个修改方案,在代码中进行说明和感谢。

您好,感谢作者开源 在使用时,我发现特征值特征向量计算时,出现了和[schumipush]一样的问题,检查发现,是作者做householder变换时,没有把原矩阵乘以当前循环得到的变换矩阵再传入下一次循环,而是直接把原矩阵直接传入下一次循环,因此我修改如下: 1613行

Matrix* x = M_Cut(Ri,i+1,_END_,i,i);//此处"_mat"修改为Ri

1657行

x = M_Cut(Ri,i+1,_END_,i+1,i+1);//此处"_mat"修改为Ri

用户"645770225"您好:

感谢您对本项目的关注和使用,这个问题已经在 最新发布的Matrix_Hub_v1.52 版本中得到修复,欢迎您去更新使用。