Алгоритм компактного хранения и решения СЛАУ высокого порядка - (реферат)
p>void Convert(Vector& X, Vector& Y, Vector& Z, Matrix& FE, DWORD NumTr, Matrix& Bounds, DWORD& BnCount) {int cData8[6][5] = {{0, 4, 5, 1, 7},
{6, 2, 3, 7, 0},
{4, 6, 7, 5, 0},
{2, 0, 1, 3, 5},
{1, 5, 7, 3, 4},
{6, 4, 0, 2, 1}},
cData4[4][4] = {{0, 1, 2, 3},
{1, 3, 2, 0},
{3, 0, 2, 1},
{0, 3, 1, 2}},
cData10[4][7] = {{0, 1, 2, 4, 5, 6, 3},
{0, 1, 3, 4, 8, 7, 2},
{1, 3, 2, 8, 9, 5, 0},
{0, 2, 3, 6, 9, 7, 1}},
cData[6][7],
Data[6],
l,
Num1,
Num2,
m;
DWORD i,
j,
p[6],
pp[6],
Index;
Matrix BoundList(4 * NumTr, 6);
double cx,
cy,
cz,
x1,
y1,
z1,
x2,
y2,
z2,
x3,
y3,
z3;
Bounds. ReSize(4 * NumTr, 6);
switch (CurrentType)
{
case BASE3D_4:
Num1 = 4;
Num2 = 3;
for (l = 0; l < Num1; l++)
for (m = 0; m < Num2+1; m++)
cData[l][m] = cData4[l][m];
break;
case BASE3D_8:
Num1 = 6;
Num2 = 4;
for (l = 0; l < Num1; l++)
for (m = 0; m < Num2 + 1; m++)
cData[l][m] = cData8[l][m];
break;
case BASE3D_10:
Num1 = 4;
Num2 = 6;
for (l = 0; l < Num1; l++)
for (m = 0; m < Num2+1; m++)
cData[l][m] = cData10[l][m];
}
printf("Create bounds.... \r");
for (i = 0; i < NumTr - 1; i++)
for (int j = 0; j < Num1; j++)
if (! BoundList[i][j])
{
for (l = 0; l < Num2; l++)
p[l] = FE[i][cData[j][l]];
for (DWORD k = i + 1; k < NumTr; k++)
for (int m = 0; m < Num1; m++)
if (! BoundList[k][m])
{
for (int l = 0; l < Num2; l++)
pp[l] = FE[k][cData[m][l]];
if (Test(p, pp))
BoundList[i][j] = BoundList[k][m] = 1;
}
}
for (i = 0; i < NumTr; i++)
for (j = 0; j < (DWORD)Num1; j++)
if (BoundList[i][j] == 0)
{
if (CurrentType == BASE3D_4)
{
cx = X[FE[i][cData[j][3]]];
cy = Y[FE[i][cData[j][3]]];
cz = Z[FE[i][cData[j][3]]];
}
else
if (CurrentType == BASE3D_10)
{
cx = X[FE[i][cData[j][6]]];
cy = Y[FE[i][cData[j][6]]];
cz = Z[FE[i][cData[j][6]]];
}
else
{
cx = X[FE[i][cData[j][4]]];
cy = Y[FE[i][cData[j][4]]];
cz = Z[FE[i][cData[j][4]]];
}
x1 = X[FE[i][cData[j][0]]];
y1 = Y[FE[i][cData[j][0]]];
z1 = Z[FE[i][cData[j][0]]];
x2 = X[FE[i][cData[j][1]]];
y2 = Y[FE[i][cData[j][1]]];
z2 = Z[FE[i][cData[j][1]]];
x3 = X[FE[i][cData[j][2]]];
y3 = Y[FE[i][cData[j][2]]];
z3 = Z[FE[i][cData[j][2]]];
for (l = 0; l < Num2; l++)
Data[l] = cData[j][l];
if ( ((cx-x1)*(y2-y1)*(z3-z1) + (cy-y1)*(z2-z1)*(x3-x1) + (y3-y1)*(cz-z1)*(x2-x1)
(x3-x1)*(y2-y1)*(cz-z1) - (y3-y1)*(z2-z1)*(cx-x1) - (cy-y1)*(z3-z1)*(x2-x1)) > 0)
{
if (CurrentType == BASE3D_4)
{
Data[0] = cData[j][0];
Data[1] = cData[j][2];
Data[2] = cData[j][1];
}
else
if (CurrentType == BASE3D_10)
{
Data[0] = cData[j][0];
Data[1] = cData[j][2];
Data[2] = cData[j][1];
Data[3] = cData[j][5];
Data[5] = cData[j][3];
}
else
{
Data[0] = cData[j][0];
Data[1] = cData[j][3];
Data[2] = cData[j][2];
Data[3] = cData[j][1];
}
}
for (l = 0; l < Num2; l++)
Bounds[BnCount][l] = FE[i][Data[l]];
BnCount++;
}
}
void main(int argc, char** argv)
{
char *input1,
*input2,
*input3,
*op = "",
*sw;
bool CreateFile(char*, char*, char*, char*);
printf("ANSYS->FORL file convertor. ZSU(c) 1998. \n\n");
if (argc < 5 || argc > 6)
{
PrintHeader();
return;
}
sw = argv[1];
input1 = argv[2];
input2 = argv[3];
input3 = argv[4];
if (! strcmp(sw, "-t10"))
CurrentType = BASE3D_10;
else
if (! strcmp(sw, "-c8"))
CurrentType = BASE3D_8;
else
{
printf("Unknown switch %s\n\n", sw);
PrintHeader();
return;
}
if (argc == 6)
{
op = argv[5];
if (strcmp(op, "/8") && strcmp(op, "/6"))
{
printf("Unknown options %s\n\n", op);
PrintHeader();
return;
}
}
if (CreateFile(input1, input2, input3, op))
printf("OK\n");
}
bool CreateFile(char* fn1, char* fn2, char* fn3, char* Op)
{
FILE *in1,
*in2,
*in3;
Vector X(1000),
Y(1000),
Z(1000);
DWORD NumPoints,
NumFE,
NumBounds = 0,
tmp;
Matrix FE(1000, 10),
Bounds;
bool
ReadTetraedrData(char*, char*, FILE*, FILE*, Vector&, Vector&, Vector&, Matrix&, DWORD&, DWORD&),
ReadCubeData(char*, char*, FILE*, FILE*, Vector&, Vector&, Vector&, Matrix&, DWORD&, DWORD&);
void Convert824(Matrix&, DWORD&),
Convert1024(Matrix&, DWORD&);
if ((in1 = fopen(fn1, "r")) == NULL)
{
printf("Unable open file %s", fn1);
return false;
}
if ((in2 = fopen(fn2, "r")) == NULL)
{
printf("Unable open file %s", fn2);
return false;
}
if (CurrentType == BASE3D_10)
{
if (! ReadTetraedrData(fn1, fn2, in1, in2, X, Y, Z, FE, NumPoints, NumFE)) return false; if (! strcmp(Op, "/8"))
{
// Create 8*Tetraedr(4)
Convert1024(FE, NumFE);
}
Convert(X, Y, Z, FE, NumFE, Bounds, NumBounds);
