Adding a new integrator type

[thomasJacob]

I've been thinking about adding a transient integrator in which the constants c1, c2 and c3 are not double scalars anymore. In fact, they are matrices. Is it possible to overload the addKtToTang, addCtoTang and addMtoTang methods with something like addKtToTang(Matrix *fact), addCtoTang(Matrix *fact) and addMtoTang(Matrix *fact) owing to the fact that their parameter types are different? These matrices also affect P* vector. How can I handle that?

[fmk]

you might look at overriding the formTangent and formResidual methods in Transient integrator .. if you looped over the FE_Elements directly you could get the tangent, do whatever matrix operation you want and add to the SOE.