Ejemplo n.º 1
0
        public static AlphaBeta CreateFrom(ISection1D <IsotropyLayer> section1D)
        {
            int numberOf1DLayers = section1D.NumberOfLayers;
            var zeta             = section1D.GetAllSection1DZeta();

            var alpha1 = Enumerable.Repeat(Complex.One, numberOf1DLayers).ToArray();
            var beta1  = Enumerable.Repeat(Complex.One, numberOf1DLayers).ToArray();

            var alphaSigma = new Complex[numberOf1DLayers];
            var betaSigma  = new Complex[numberOf1DLayers];

            for (int k = 0; k < numberOf1DLayers - 1; k++)
            {
                alphaSigma[k] = zeta[k + 1] / zeta[k];
            }

            for (int k = 1; k < numberOf1DLayers; k++)
            {
                betaSigma[k] = zeta[k - 1] / zeta[k];
            }

            return(new AlphaBeta(alpha1, beta1, alphaSigma, betaSigma));
        }