Adding a new integrator type

For developers writing C++, Fortran, Java, code who have questions or comments to make.

Moderators: silvia, selimgunay, Moderators

Post Reply
Posts: 3
Joined: Fri Sep 08, 2017 7:18 am

Adding a new integrator type

Post by thomasJacob » Fri Sep 08, 2017 7:40 am

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?

Site Admin
Posts: 5866
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley

Re: Adding a new integrator type

Post by fmk » Thu Oct 19, 2017 8:00 pm

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.

Post Reply