構造体内の配列

Question

構造体に埋め込まれていると、C言語の新機能で配列型に問題があります。以下は私の問題の一例です：

typedef struct { 
    double J[151][151]; 
} *UserData; 

static int PSolve(void *user_data, N_Vector solution) 
{ 
UserData data; 
data = (UserData) user_data; 

double J[151][151]; 
J = data->J; 

/* Solve a matrix equation that uses J, stored in 'solution' */ 

return(0); 
} 

私はこれをコンパイルしようとすると、私は エラーが出る：型に割り当てる際のタイプから「ダブル[151] [151]」互換性のない型「をダブル（*）[ 151] '

現在の回避策は、行列式を解くために' J [x] [y] 'をコードの' data-> J [x] [y] 'に置き換えることでしたが、効率が悪いことを示しています。

引数の型と順序を規定するsundials-cvodeソルバーを使用しているため、PSolveへの引数の変更はオプションではありません。

は、助けを求め アンドリューありがとう

Answer 1

typedef struct { 
    double J[151][151]; 
} UserData; // This is a new data structure and should not a pointer! 

static int PSolve(void *user_data, N_Vector solution) 
{ 
UserData* data; // This should be a pointer instead! 
data = (UserData*) user_data; 

double J[151][151]; 
memcpy(J, data->J, sizeof(double) * 151 * 151); // use memcpy() to copy the contents from one array to another 

/* Solve a matrix equation that uses J, stored in 'solution' */ 

return(0); 
} 

Answer 2

double **Jとしてローカル変数double J[151][151]を定義してみます。

Answer 3

配列に割り当てることはできません。

double (*J)[151] = data->J; 

長151の配列へのポインタです：あなたは、あなたのコード内で必要なもの

C.

における配列とポインタの違いについては、このようなだけで何かで調べる必要があります。あなたがすべてだのtypedef

typedef double line[151]; 
line *J = data->J; 

てみたい場合は、あなたはデータをコピーしますが、データへのポインタだけではないはずです。

編集：しかし、回答のスレッド全体を見て、私はそのすべてがあなたのボトルネックがあるかもしれない、純粋な憶測だと思います。たとえば、あなたが「間違った方法で丸ごと」、列方向にアクセスしているなど、どこにいてもかまいません。あるいは、メモリからデータをポンピングすることは、あなたの計算を支配しているだけです。

おそらく、あなたのコンパイラが生成するアセンブラ（オプション-S）を調べて、何か怪しいものがあるかどうかを調べることを検討してください。

Answer 4

問題の根本原因は、Cで配列タイプを直接割り当てることができないことです。memcpy()を明示的に使用する必要があります。karlphillip's answerが表示されます。

ただし、コピーを実行すると、残りの機能で行った最適化の効果がなくなる可能性があります。恐らくsolution引数はポインタであり、オプティマイザはuser_data/dataとsolutionの間の潜在的なエイリアシングを心配しています。あなたはC99のrestrictキーワードをサポートコンパイラを持っている場合は別の方法として、引数にその修飾子を使用することです：

static int PSolve(void * restrict user_data, N_Vector restrict solution) 

これは、これらのポインタがないエイリアスを行うコンパイラを約束し、あなたが直接できるようにする必要がありますコンパイラの最適化を犠牲にすることなくdata->Jを使用してください。

キーワードは、__restrictなどのスペルでC89モードで使用できます。詳細については、コンパイラのドキュメントを参照してください。