Blender Git Loki

Cleanup, removed unused old code.

Optimized matrix filling using the Eigen triplets method.

Otherwise the construction of matrices becomes very slow for larger
vertex counts because adding a new element is O(n), making it O(n^2) in
total.

Added back external effector forces for the cloth solver.

Reenabled goal springs and removed the arbitrary structural scaling
factor.

Fixed bending forces for cloth springs.

Correctly implemented and verified gravity, drag, structural springs.

Minor fix in debug print, so solver prints can be diff'ed.

A couple of defines to disable various influences on the cloth sim for
debugging.

Debug printing code for large sparse vectors and matrices, to compare
solver input/output of the old and new methods.

Added a new internal header for implicit solver defines.

Should become its own intern wrapper library eventually.

Reenabled air drag force for cloth sim.

Arbitrary scaling factor for the structural spring forces (stretching).

This is not nice at all, but for some reason (possibly time scale) the
old force values are much too high and cause the solver to become
unstable. These will be revisited later anyway, so for now such scaling
should be fine.

Switched to the modified CG method that supports constraints, and
added back structural stretch springs.

Implemented a modifier conjugate gradient algorithm to support stable
constraining.

The algorithm is described in the paper "Large Steps in Cloth Simulation"
(Baraff/Witkin 1998). The same method was (incorrectly) implemented in
the old cloth solver.

It is based on restricting the degrees of freedom (ndof) of vertices
using a block matrix and a vector of target velocity deltas.
See chapter 5 of the paper for details.

Reimplemented Goal springs for the Eigen CG solver method.

Note that goal springs currently are really bad ... They have a factor
on hairs that "fades" goal influence from the root to the tip. The last
point on the hair is completely free, which makes the goal springs
pretty much useless on their own without supporting bend stiffness.
Can only assume this was added to compensate unphysical behavior of
goal springs when using uniform weight, but it's a poor replacement for
true localized bending forces ...

To make gravity a true force vector, multiply by the mass.

Alternative new solver for cloth using the Eigen CG solver instead of
a custom built solver.

The old cloth solver is broken unfortunately. Eigen is a designated
linear algebra library and very likely their implementation is a lot
better (can't compare until it's implemented though).

Only basic gravity is active atm, spring forces, external force fields,
damping and volumetric friction have to be added back by converting
the data into the Eigen format.

Don't try to enforce a velocity change for pinned vertices; their
velocity is defined externally by hair root motion.

Use squareroot of the current epsilon error target for CG, since we
square this value anyway.

Disabled collision culling on the inside of the collider faces for now,
this seems to remove too many contact points somehow ...