CLAPACKの動作について(121)

1 名無し＠タコグラマー sage 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;
}

総レス数 121

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

---

IE等普通のブラウザで見る場合 http://tubo.80.kg/tubo_and_maru.html
専用のブラウザで見る場合 http://www.monazilla.org/

２ちゃんねる　Viewerを使うと、すぐに読めます。 http://2ch.tora3.net/
この Viewer(通称●) の売上で、２ちゃんねるは設備を増強しています。
●が売れたら、新しいサーバを投入できるという事です。

よくわからない場合はソフトウェア板へGo http://hibari.2ch.net/software/

モリタポを持っていれば、50モリタポで表示できます。

---

read.cgi@hanako ver 2011/01/22
FOX ★ DSO(Dynamic Shared Object)