Холеский openMP в C - PullRequest
       4

Холеский openMP в C

0 голосов
/ 30 января 2020

В следующем коде реализован факторинг Холески, где для таблицы 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));
}
}
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...