5ちゃんねる ★スマホ版★ ■掲示板に戻る■ 全部 1- 最新50  

■ このスレッドは過去ログ倉庫に格納されています

CLAPACKの動作について

1 :名無し@タコグラマー:01/11/01 15:09
CLAPACK利用して以下のプログラムを走らせると
正しい答えだしてくれませんん。
何がわるいかわかりますか?

ちなみに、一般帯行列Aを利用した[A]{x}={b}
を解いています。
サブルーチンにいれる行列形式は帯格納形式です。

#include <stdio.h>
#define N 4
int main(void)
{
double A[N*N];
double x[N];
int i;
long int n=N,kl=1,ku=1,nrhs=1,ldab,ldb,info,piv[N];
ldab=2*kl+ku+1;
ldb=N;
A[0]=0.;A[1]=0.;A[2]=0.;A[3]=0.;
A[4]=0.;A[5]=-1.;A[6]=-1.;A[7]=-1.;
A[8]=2.;A[9]=2.;A[10]=2.;A[11]=2.;
A[12]=-1.;A[13]=-1.;A[14]=-1.;A[15]=0.;
x[0]=5.;x[1]=5.;x[2]=5.;x[3]=5.;
printf("N = %d\n",N);

dgbsv_(&n,&kl,&ku,&nrhs,A,&ldab,piv,x,&ldb,&info);

printf("info= %ld\n",info);
for(i=0;i<N;++i) printf("%f\n",x[i]);
return 0;
}

126 :名無しさん@5周年:2006/06/26(月) 23:16:51

【数学】フェルマーの最終定理に反例が・・・ワイルズ会見
http://news18.2ch.net/test/read.cgi/news7/1149592020/



36 KB
■ このスレッドは過去ログ倉庫に格納されています

★スマホ版★ 掲示板に戻る 全部 前100 次100 最新50

read.cgi ver 05.04.00 2017/10/04 Walang Kapalit ★
FOX ★ DSO(Dynamic Shared Object)