В следующем коде реализован факторинг Холески, где для таблицы A применяется A = LL *, где L - нижнее три angular, а L * - его сопряженное обратное. Может ли кто-нибудь подтвердить, если я прав, потому что я не уверен, что понял, когда мне следует использовать omp. (Подготовка к экзамену)
for(j=0;j<n;j++){
s=0;
for(k=0;k<j;k++)
s+=L[j][k] * L[j][k];
L[j][j] = sqrt(A[j][j]-s);
#pragma omp parallel for
for(i=j+1;i<n;i++){
m=0;
for(k=0;k<j;k++)
m+=L[i][k]*L[j][k];
L[i][j]=(1.0/L[j][j]*(A[i][j]-m));
}
}