遅い

Question

Aで解決の任意の選択肢は、（）国の各産業が他の国の他の産業に関連しているどのくらいを表します。それは約1.00から値を取ることができます。したがって、行列Aの対角要素は1に近いが、行列Aのオフ対角要素は1よりもはるかに小さい。なぜなら、各業界は主にそれ自体に関連しているからである。

I、Jはそれぞれ1:44の国を表します。
注意！ （ "I" で（I-A_omitは（J、J）））DIAG（2464で解決する。
i、jは1：56それぞれの産業を表す。

Xが出力されます。長さのベクトル44の国X 56の業種= F（特定の国における産業の最終需要を表す2464 * 1行列である行列を乗じ
それは（DIAG（2464）-A）の逆数によって算出される2464年）

ここで、行列Xのj + 56 *（J-1）番目の列は、
ところで、解く（I-A_omit（J、j））の "I"は診断（2464）。 （solveを使用して）X_omitの

1つの計算は、約20秒を要します。 Aは2464 * 2464マトリックスで、Xは2464 * 1マトリックスです。 計算すると仮定しますb_linkは約20秒* 2464の計算時間を要します。これは1つの単一値に対して14時間です。
しかし、この場合には、私は問題ではありませんIとiを分離し、なぜならこの場合にはsumを使用して簡単な式に変更するために管理しています。 しかし、次の計算では、I,i,J,jのサイズが1:44,1:56,1:44,1:56であり、合計を使用するだけでなく、各文字を区切る必要があります。

私はより速くそれを作るための唯一の方法は、より高速な20秒かかりsolve()を作ることによってだと思います。 solve(I-A_omit(J,j))の[A]行列は対称ではなく、通常の行列には数字が入っています。

これはあなたがループの外X_omit(J,j)を計算することができます私のコード

X_omit <- function (J,j) {solve((I-A_omit(J,j)),F)} 
B_linkages_inter <- function (I,J) { 
    for (j in 1:56) { 
    for (i in 1:56) { 
     sum((X[i+56*(I-1)]-X_omit(J,j)[i+56*(I-1)])/X[j+56*(J-1)]) 
    } 
    } 
} 

Answer 1

です。また、ループ外でX[j+56*(J-1)])を計算することもできます。