right-look LU分解

最新推荐文章于 2024-10-06 21:57:13 发布

JMAVI

最新推荐文章于 2024-10-06 21:57:13 发布

阅读量613

点赞数

CC 4.0 BY-SA版权

分类专栏：计算电磁学文章标签： ini

本文链接：https://blog.youkuaiyun.com/wangjiannuaa/article/details/7025164

计算电磁学专栏收录该内容

41 篇文章

订阅专栏

本文详细阐述了矩阵右向看算法的实现过程，包括更新L矩阵和U矩阵的步骤，涉及数学运算和迭代操作。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

void Right_Looking(double**& A,int N)
{
   for (int i=0;i<N;i++)
   {
       for (int j=i+1;j<N;j++){
           A[j][i] = A[j][i]/A[i][i];//update L

           for (int k=i+1;k<N;k++)
               A[j][k] -= A[i][k]*A[j][i]; //update U
       }
   }
}

void Right_Looking2(double**& A,int N)
{
   for (int k=0;k<N-1;k++)
   {
       //A(k+1:N,k) = A(k+1:N,k) / A(k,k)
       for (int i=k+1;i<N;i++){
           A[i][k] = A[i][k]/A[k][k];//update L
       }
       //A(k+1:N,k+1:N) -= A(k+1:N,k)*A(k,k+1:N)
       for(int i=k+1;i<N;i++){
           for (int j=k+1;j<N;j++)
               A[i][j] -= A[i][k]*A[k][j]; //update U
       }
   }
}