Exemplo n.º 1
0
        public static CompressedColumnStorage <double> ConvertSparsityJacobian(AlgebraicSystem problem)
        {
            var sparseJacobian = new CoordinateStorage <double>(problem.NumberOfEquations, problem.NumberOfVariables, problem.Jacobian.Count);

            foreach (var entry in problem.Jacobian)
            {
                sparseJacobian.At(entry.EquationIndex, entry.VariableIndex, 1);
            }

            var compJacobian = CSparse.Converter.ToCompressedColumnStorage(sparseJacobian);

            return(compJacobian);
        }
Exemplo n.º 2
0
        public static CompressedColumnStorage <double> FillResidualAndJacobian(AlgebraicSystem system, Vector b)
        {
            var sparseJacobian = new CoordinateStorage <double>(system.NumberOfEquations, system.NumberOfVariables, system.Jacobian.Count);

            int    j    = -1;
            double grad = -1;

            for (int i = 0; i < system.Equations.Count; i++)
            {
                system.Equations[i].Reset();
                b[i] = system.Equations[i].Residual();

                foreach (var variable in system.Equations[i].Variables)
                {
                    if (system.VariableIndex.TryGetValue(variable, out j))
                    {
                        grad = system.Equations[i].Expression.Diff(variable);
                        sparseJacobian.At(i, j, grad);
                    }
                }
            }

            /*  int j = -1;
             * double grad = -1;
             *
             * Parallel.For(0, system.Equations.Count, i =>
             * {
             *    system.Equations[i].Reset();
             *    b[i] = system.Equations[i].Residual();
             *
             *    foreach (var variable in system.Equations[i].Variables)
             *    {
             *        if (system.VariableIndex.TryGetValue(variable, out j))
             *        {
             *            grad = system.Equations[i].Expression.Diff(variable);
             *            sparseJacobian.At(i, j, grad);
             *        }
             *    }
             * });*/


            return(CSparse.Converter.ToCompressedColumnStorage(sparseJacobian));
        }
Exemplo n.º 3
0
        public void Initialize(ILogger logger)
        {
            _logger = logger;

            StateValues   = new double[DifferentialVariables.Count];
            SystemToSolve = new AlgebraicSystem(Name);


            foreach (var vari in AlgebraicVariables)
            {
                SystemToSolve.AddVariables(vari);
            }
            foreach (var vari in DifferentialVariables)
            {
                SystemToSolve.AddVariable(vari);
            }
            foreach (var eq in Equations)
            {
                SystemToSolve.AddEquation(eq);
            }

            SystemToSolve.CreateIndex();
            SystemToSolve.GenerateJacobian();
        }