Beispiel #1
0
        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);
        }
Beispiel #2
0
        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);
        }