ru.algorithms

 
 - RU.ALGORITHMS ----------------------------------------------------------------
 From : Eugene Belogurov                     2:461/1148     04 Apr 2003  02:28:02
 To : All
 Subject : СЛАУ...
 -------------------------------------------------------------------------------- 
 

 
 Решаю сабж, только немного расширенную: не матрица*столбец=столбец, а
 4-х мерная матрица*матрица = матрица.
 
  I   Q
 -+- ---                        q=1..Q
 \   \    Yijpq * Vip = Fqj     j=1..I
 /__ /__
 i=1  p=1
 
 Пытаюсь это делать так:
 
 //  Y2(Q*I,Q*I)
 //  V1(I*Q),F1(Q*I)
 // 4->2
 
    for(i=0;i<I;++i)
     for(j=0;j<I;++j)
      for(p=0;p<Q;++p)
       for(q=0;q<Q;++q) Y2(Q*i+p,Q*j+q) = Y(i,j,p,q);
 
 //2-> 1
 
    for(i=0;i<I;++i)
     for(p=0;p<Q;++p) F1[Q*i+p] = F(i,p);
 
 Далее решаю это дело и получаю матрицу решений из столбца:
 
    for(i=0;i<I;++i)
     for(p=0;p<Q;++p) V(i,p) = V1[Q*i+p];
 
 Так вот вопросы у меня: правильное решение получается таким образом или нет?
 Есть ли более эффективный способ решения?
 Eugene.
 
 ... Hе грусти,медитиламинотрифосфодиметилфенилпироизонафтола отглотни!
 --- Весна: [ЫЫЫЫЫЫЫЫЫЫЫЫ°°°°°°°°°°°°°°°°°°°°°] 37.46% complete
  * Origin: астрой букву "" !!! (2:461/1148)
 
 

Вернуться к списку тем, сортированных по: возрастание даты  уменьшение даты  тема  автор