Adding a new integrator type

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

Moderators: silvia, Moderators

Adding a new integrator type

Postby 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?
thomasJacob
 
Posts: 2
Joined: Fri Sep 08, 2017 7:18 am

Re: Adding a new integrator type

Postby 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.
fmk
Site Admin
 
Posts: 5755
Joined: Fri Jun 11, 2004 2:33 pm
Location: UC Berkeley


Return to Framework

Who is online

Users browsing this forum: No registered users and 1 guest

cron