On Fri, Jun 20, 2014 at 10:53 AM, Sahai, Amal (ARC-TN)[Universities Space Research Association (USRA)] <amal.sa...@nasa.gov> wrote: > Hi what would be the best way of solving a second order linear partial > differential equation in time of the type in LibMesh: > > [M] d^2u/dt^2 = [K] u + [F] > > I could obviously break it down into sets of equations - one solving for the > first derivative and other for u itself but I'd like to know if there is an > easier way to implement this. Also, how does one account for [M]?
I think this is what the NewmarkSystem is good for? It allows you to specify mass, damping, and stiffness matrices... see transient_ex2 for an example. -- John ------------------------------------------------------------------------------ HPCC Systems Open Source Big Data Platform from LexisNexis Risk Solutions Find What Matters Most in Your Big Data with HPCC Systems Open Source. Fast. Scalable. Simple. Ideal for Dirty Data. Leverages Graph Analysis for Fast Processing & Easy Data Exploration http://p.sf.net/sfu/hpccsystems _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
