Предыдущая страница реферата | 17  18  19  20  21  22  23  24  25  26  27 | Следующая страница реферата

DWORD RRow;

for (DWORD i = 0L; i < FE.Size(); i++) for (DWORD l = 0L; l < Dim; l++) for (DWORD j = 0L; j < FE.Size(); j++) for (DWORD k = 0L; k < Dim; k++)

{

Res = FEMatr[i * Dim + l][j * Dim + k]; if (Res) Add(FE[i],l,FE[j],k,Res);

} for (DWORD i = 0L; i < FE.Size(); i++) for (DWORD l = 0L; l < Dim; l++)

{

RRow = FE[UINT(i % (FE.Size()))] * Dim + l;

Res = FEMatr[i * Dim + l][FEMatr.Size1()]; if (Res) Add(RRow,Res);

}

}

DWORD TSMatrix::Find(DWORD I,DWORD J)

{

DWORD i;

for (i = 0; i < Links[I].Size(); i++) if (Links[I][i] == J) return i; return DWORD(-1);

}

void TSMatrix::Restore(Matrix& Matr)

{

DWORD i, j,

NRow,

NPoint,

NLink,

Pos;

Matr.ReSize(Size * Dim,Size * Dim + 1); for (i = 0; i < Size; i++) for (j = 0; j

{

NRow = j / (Array[i].Size() / Dim); // Number of row

NPoint = (j - NRow * (Array[i].Size() / Dim)) / Dim; // Number of points

NLink = j % Dim; // Number of link

Pos = Links[i][NPoint];

Matr[i * Dim + NRow][Pos * Dim + NLink] = Array[i][j];

} for (i = 0; i < Right.Size(); i++) Matr[i][Matr.Size1()] = Right[i];

}

void TSMatrix::Set(DWORD Index,DWORD Position,double Value,bool Case)

{

DWORD Row = Index,

Col = Position * Links[Index].Size() * Dim +
Find(Index,Index) * Dim + Position, i; double koff = Array[Row][Col], val;

if (!Case)

Right[Dim * Index + Position] = Value; else

{

Right[Index * Dim + Position] = Value * koff; for (i = 0L; i < Size * Dim; i++) if (i != Index * Dim + Position)

{

Set(Index * Dim + Position,i,0);

Set(i,Index * Dim + Position,0); if (Get(i,Index * Dim + Position,val))

Right[i] -= val * Value;

}

}

}

void TSMatrix::Mul(RVector& Arr,RVector& Res)

{

DWORD i, j,

NRow,

Рекомендуем скачать другие рефераты по теме: бесплатно рассказы, бесплатно реферат на тему.

Предыдущая страница реферата | 17  18  19  20  21  22  23  24  25  26  27 | Следующая страница реферата

Поделитесь этой записью или добавьте в закладки