Re: Help with Differentiation

John Fletcher <J.P.Fletcher <at> aston.ac.uk> writes:


> Just playing at the moment. I can see a lot of uses for this.

The reason why I have written Daixtrose is: I had to implement a time-implicit
FEM-FCT method (as described by Kuzmin,
http://www.mathematik.uni-dortmund.de/mitarbeiter_Dmitri_Kuzmin_SS05.html
e.g. in http://www3.interscience.wiley.com/cgi-bin/abstract/104530112/ABSTRACT
or http://portal.acm.org/citation.cfm?id=585284) for the solution of systems of
partial differential equations like the Euler equations or the equations for a
thermoelastic linear or geometrically nonlinear solid in 2D with up to 23
partial differential equations. Since the equations I had to investigate
changed a lot I finally got fed up and created an Interface that allows C++
code of the Form

Solver.AddPDE<1>(Diff<Time>(u) + Diff<SpaceX>(Pow<2,1>(u)) + ... == 0);
etc...

which automagically generates a Newton-Raphson Time-Intergration scheme and a
solver which outputs a latex report about the equation solved. I still hope to
use this for my doctoral thesis, so please forgive for not publishing the code
- yet.

Since you are working on Clifford Algebra: I once had correspondence with Jaap
Suter who wrote a geometric algebra package based on Expression Templates. I
could not find it anymore. But this code might be inspiring, too. So probably
contact him directly.

best regards,

Markus





-------------------------------------------------------
This SF.net email is sponsored by: Splunk Inc. Do you grep through log files
for problems?  Stop!  Download the new AJAX search engine that makes
searching your log files as easy as surfing the  web.  DOWNLOAD SPLUNK!
http://ads.osdn.com/?ad_id=7637&alloc_id=16865&op=click