版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

北京石油化工学院2026年研究生招生接收调剂公告

返回列表

wdphysics

金虫 (小有名气)

应助: 0 (幼儿园)
金币: 1110.6
帖子: 114
在线: 90.9小时
虫号: 1000002
注册: 2010-04-18
专业: 凝聚态物性 II ：电子结构

[交流] 【求助】如何调用此函数分解具体矩阵？已有1人参与

对称矩阵的cholesky分解，a为矩阵，n为其维数。主要是不知道怎么弄调用接口！
void choldc(float **a, int n, float p[])
{
void nrerror(char error_text[]);
int i,j,k;
float sum;
for (i=1;i<=n;i++) {
for (j=i;j<=n;j++) {
for (sum=a[j],k=i-1;k>=1;k--) sum -= a[k]*a[j][k];
if (i == j) {
if (sum <= 0.0)
nrerror("choldc failed";
p=sqrt(sum);
} else a[j]=sum/p;
}
}
}

回复此楼

» 猜你喜欢

求老师收留已经有4人回复
一志愿郑州大学材料与化工085600，求调剂已经有11人回复
专硕085601求调剂已经有3人回复
重庆大学材料与化工085600，初试370+，求求调剂建议已经有7人回复
286求调剂已经有3人回复
260求调剂已经有3人回复
326求调剂已经有4人回复
调剂已经有4人回复
285求调剂已经有6人回复
材料与化工306分找调剂已经有3人回复

1楼 2010-08-20 10:38:10

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

yalefield

金虫 (文坛精英)

老汉一枚

程序强帖: 3
应助: 129 (高中生)
贵宾: 0.17
金币: 21238.9
散金: 3440
红花: 66
帖子: 12101
在线: 759.1小时
虫号: 96063
注册: 2005-10-07
专业: 高等教育学
管辖: 计算模拟

★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
余泽成(金币+3):辛苦了！ 2010-08-20 12:18:33

（1）
a是指针的指针，在这里，它就是指向一个2维的数组
C/C++里面，数组都是一维的，遇到二维数组，就把它“展开”成一维的。
因此，你要指定二维数组中，那个第一维有多少个元素，也就是参数n
   第0到第n-1个，是二维数组中的前n个，也就是第一行（或列）
p则是一维数组

（2）
因此，要调用这个函数，需要准备好一个二维数组和一个一维数组
例如：
int f ()
{
      float X[100][200];
      float Y[100];
      return choldc(X, 100, Y);
}

（3）
按照目前的写法，应该是：int choldc(float **a, int n, float *p)
然后，返回1表示成功，返回0表示出错（或者相反）。
在这样的函数里，直接调用错误输出，实在算不上的德艺双馨。

CODE:

int choldc(float **a, int n, float p[])

{

        int i, j, k;

        float sum = 0.0F;

        for (i=1; i<=n; i++) 

        {

                for (j=i; j<=n; j++) 

                {

                        for (sum=a[i][j], k=i-1; k>=1; k--)

                        {

                                sum -= a[i][k]*a[j][k];

                        }

                        if (i == j) 

                        {

                                if (sum <= 0.0)

                                {

                                        return 0;

                                }

                                p[i]=sqrt(sum);

                        } 

                        else

                        {

                                a[j][i]=sum/p[i];

                        }

                }

        }

        return 1;

}

-------------------
C语言九星大师，偶发三俗代码

赞一下(2人)