public static double[,,] Generate(LabelPoint[,,] finiteElement)
{
double[,,] result = new double[3, 3, 27];
_3DPoint[,,] matrix = new SquareGenerator(-1 * Math.Sqrt(0.6), 0, Math.Sqrt(0.6)).getMatrix();
Dictionary <int, _3DPoint> magicDictionary = Globals.magicDictionary;
//// to delete
//double[] oksana = new double[] { -1 * Math.Sqrt(0.6), 0, Math.Sqrt(0.6) };
//// to delete
//for (int i = 0; i < 3; i++)
//{
// for (int j = 0; j < 3; j++)
// {
// for (int k = 0; k < 3; k++)
// {
// lst.Add(new _3DPoint(oksana[i], oksana[j], oksana[k]));
// }
// }
//}
List <_3DPoint> lst = new List <_3DPoint>();
for (int z = 0; z < 3; z++)
{
for (int y = 0; y < 3; y++)
{
for (int x = 0; x < 3; x++)
{
lst.Add(matrix[x, y, z]);
}
}
}
_3DPoint p = null;
double globalCoordinate = 0;
double diFi = 0;
double sum = 0;
for (int i = 0; i < 3; i++)
{
for (int j = 0; j < 3; j++)
{
for (int k = 0; k < 27; k++)
{
sum = 0;
for (int l = 1; l <= 20; l++)
{
p = magicDictionary[l];
globalCoordinate = finiteElement[(int)p.X, (int)p.Y, (int)p.Z].value[i];
diFi = FI.getDiFi(j + 1, l, lst[k].X, lst[k].Y, lst[k].Z);
sum += globalCoordinate * diFi;
}
result[i, j, k] = sum;
}
}
}
return(result);
}
public static double[,,] Generate()
{
double[,,] result = new double[27, 3, 20];
_3DPoint[,,] matrix = new SquareGenerator(-1 * Math.Sqrt(0.6), 0, Math.Sqrt(0.6)).getMatrix();
List <_3DPoint> lst = new List <_3DPoint>();
//double[] otherCountWay = new double[3] { -1 * Math.Sqrt(0.6), 0, Math.Sqrt(0.6) };
Dictionary <int, _3DPoint> magicDictionary = Globals.magicDictionary;
for (int z = 0; z < 3; z++)
{
for (int y = 0; y < 3; y++)
{
for (int x = 0; x < 3; x++)
{
lst.Add(matrix[x, y, z]);
//lst.Add(new _3DPoint(otherCountWay[z], otherCountWay[y], otherCountWay[x]));
}
}
}
for (int i = 0; i < 27; i++)
{
for (int j = 0; j < 3; j++)
{
for (int k = 0; k < 20; k++)
{
result[i, j, k] = FI.getDiFi(j + 1, k + 1, lst[i].X, lst[i].Y, lst[i].Z);
//_3DPoint point = magicDictionary[i+1];
//_3DPoint coord = matrix[(int)point.X, (int)point.Y, (int)point.Z];
//result[i, j, k] = FI.getDiFi(j + 1, k + 1, coord.X, coord.Y, coord.Z);
}
}
}
return(result);
}