private void ApplyBCs() { KK = new double[model.LQ * model.ln, model.LQ *model.ln]; KK = MatrixTransformation.Multiply <int, double>(model.Identity.Ip, K); KK = MatrixTransformation.Multiply <double, int>(KK, model.Identity.Ip); KK = MatrixTransformation.Sum(KK, MatrixTransformation.ToDouble <int>(model.Identity.Id)); F = new double[model.LQ * model.ln, 1]; F = MatrixTransformation.Multiply <int, double>(model.Identity.Ip, model.Load.Matrix); }
private void ApplyBCs() { DoubleMatrixBoundaryGlobal = new double[SolverModel.IntDOFPerNode * SolverModel.IntModelNodes, SolverModel.IntDOFPerNode *SolverModel.IntModelNodes]; DoubleMatrixBoundaryGlobal = MatrixTransformation.Multiply(SolverModel.Identity.Ip, DoubleMatrixStiffnessGlobal); DoubleMatrixBoundaryGlobal = MatrixTransformation.Multiply(DoubleMatrixBoundaryGlobal, SolverModel.Identity.Ip); DoubleMatrixBoundaryGlobal = MatrixTransformation.Sum(DoubleMatrixBoundaryGlobal, MatrixTransformation.ToDouble(SolverModel.Identity.Id)); DoubleMatrixLoads = new double[SolverModel.IntDOFPerNode * SolverModel.IntModelNodes, 1]; DoubleMatrixLoads = MatrixTransformation.Multiply(SolverModel.Identity.Ip, SolverModel.Load.Matrix); }