可変次元を持つドットプロダクトのCマクロ

Question

数学的アルゴリズムの実装のヘッダファイルにドットプロダクトを定義する次のマクロコンストラクトがあります。

#define DIM 3 
#define MULTIPLY(a, b)  ((a)*(b)) 
#define DOT_PRODUCT(a, b) ( MULTIPLY((a)[0],(b)[0]) + \ 
           MULTIPLY((a)[1],(b)[1]) + \ 
           MULTIPLY((a)[2],(b)[2])) 

おそらく、ディメンションから実装を独立させるために使用されます。しかし明らかに、もし誰かが​​を変更したとしても、それを反映するためにDOT_PRODUCTを修正しなければなりません。

を定義して、設定された寸法​​を自動的に尊重する方法はありますか？おそらく、いくつかの（疑似）再帰的なマクロマジックで？私は何かを考え出そうとしましたが、私のC-fooは特に強くはありません。

Answer 1

私は@Olafに同意します：機能が使用できるマクロは使用しないでください。 しかし時々状況や欲求がその考えを覆してしまいます。

プリプロセッサは時には正しい場所に適切な数の迂回を望んでいるので、これを得るまで実験しました。私は、私が予測することができなかったことを認めて、少なくとも、面白い、自明ではない部分について、詳しく説明することはできません。

#define DIM 2 
#define CONCAT(x,y,a,b) x ## y(a,b) 
#define DOT_PRODUCT_DIM(y,a,b) CONCAT(DOT_PRODUCT_,y,a,b) 
#define MULTIPLY(a, b)  ((a)*(b)) 
#define DOT_PRODUCT_2(a, b) ( MULTIPLY((a)[0],(b)[0]) + \ 
           MULTIPLY((a)[1],(b)[1])) 
#define DOT_PRODUCT_3(a, b) ( MULTIPLY((a)[0],(b)[0]) + \ 
           MULTIPLY((a)[1],(b)[1]) + \ 
           MULTIPLY((a)[2],(b)[2])) 
#define DOT_PRODUCT_4(a, b) ( MULTIPLY((a)[0],(b)[0]) + \ 
           MULTIPLY((a)[1],(b)[1]) + \ 
           MULTIPLY((a)[2],(b)[2]) + \ 
           MULTIPLY((a)[3],(b)[3])) 
#define DOT_PRODUCT(a,b) DOT_PRODUCT_DIM(DIM,a,b) 

DOT_PRODUCT_DIM(2,a,b) 
DOT_PRODUCT_DIM(3,a,b) 
DOT_PRODUCT_DIM(4,a,b) 
DOT_PRODUCT_DIM(DIM,a,b) 
DOT_PRODUCT(a,b) 

OUTPUT（gcc -E Toy.c）のみ#define DIM 4を変更した後

((((a)[0])*((b)[0])) + (((a)[1])*((b)[1])) + (((a)[2])*((b)[2]))) 
((((a)[0])*((b)[0])) + (((a)[1])*((b)[1])) + (((a)[2])*((b)[2]))) 

OUTPUT（（出力の最後の2行を引用のみ#define DIM 3を変更した後、）

((((a)[0])*((b)[0])) + (((a)[1])*((b)[1]))) 
((((a)[0])*((b)[0])) + (((a)[1])*((b)[1])) + (((a)[2])*((b)[2]))) 
((((a)[0])*((b)[0])) + (((a)[1])*((b)[1])) + (((a)[2])*((b)[2])) + (((a)[3])*((b)[3]))) 
((((a)[0])*((b)[0])) + (((a)[1])*((b)[1]))) 
((((a)[0])*((b)[0])) + (((a)[1])*((b)[1]))) 

OUTPUT出力の最後の2行のみを引用する）：

((((a)[0])*((b)[0])) + (((a)[1])*((b)[1])) + (((a)[2])*((b)[2])) + (((a)[3])*((b)[3]))) 
((((a)[0])*((b)[0])) + (((a)[1])*((b)[1])) + (((a)[2])*((b)[2])) + (((a)[3])*((b)[3]))) 

（出力のみの最後の2行引用のみ#define DIM 5を変更した後、）OUTPUT：もちろん

DOT_PRODUCT_5(a,b) 
DOT_PRODUCT_5(a,b) 

、定義されなければならないDOT_PRODUCT_number()中間マクロを。私が提供した例は仕事であり、体系的な問題はなく、より高い次元を妨げることになります。
ユースケースに一致することを願っています。あなたは不確定な、静的に未知/予測不可能な多数の次元を考えていませんか？