// This is a symetric matrix: m=n (= global m) /* * J = m*n * W = n*n * JT = n*m * WJT = n*m * JWJT = m*m */ public Eq11LHS(BlockSparseMatrix a_J, double[] a_W) { m_J = a_J; m_W = a_W; m_Size = a_J.GetM(); //=m m_TempVector = new double[m_W.Length]; //=n }
public ParticleSystem(LinearSolver a_Solver, float a_ConstraintSpringConstant, float a_ConstraintDampingConstant) { m_Solver = a_Solver; if (m_Solver == null) { throw new Exception("Please provide a valid solver."); } m_Particles = new List <Particle>(); m_Forces = new List <Force>(); m_Constraints = new List <Constraint>(); m_J = new BlockSparseMatrix(); m_JDot = new BlockSparseMatrix(); m_ConstraintSpringConstant = a_ConstraintSpringConstant; m_ConstraintDampingConstant = a_ConstraintDampingConstant; }