The Method of Lines thorn (version 2 - see below).

MoL provides generic integration methods for multiple thorns simultaneously. By providing a layer between the driver and evolution thorns, this should mean that some technical issues to do with mesh refinement can be ignored. It also allows you to choose different evolution methods (in time, at least). But the primary purpose is to unambiguously evolve models in different thorns at the same time.

This is version 2 - the one that will work with mesh refinement. It's a straight copy of HawkeCVS/Public/CactusMoL2/MoL2 and I haven't checked that it will work "as is" with the new name (one of the reasons it goes into Alpha).

At the moment the only evolution method guaranteed to work with mesh refinement is ICN, although RK2 should. The generic RK methods will do something that will be subtly wrong...

Note that the "old" way of registering variables (through the functions declared in header files) is still there. As soon as function aliasing settles down, this will be removed.

Also to be done:

Better documentation
Tidy up the code (especially the debugging statements)
Optimize
Add various useful time evolution methods.


git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/MoL/trunk@2 578cdeb0-5ea1-4b81-8215-5a3b8777ee0b

1 files changed, 27 insertions, 0 deletions

diff --git a/src/MoLFunctions.h b/src/MoLFunctions.h
new file mode 100644
index 0000000..9f99e73
--- /dev/null
+++ b/src/MoLFunctions.h
@@ -0,0 +1,27 @@
+ /*@@
+   @file      MoLFunctions.h
+   @date      Thu Jun 13 22:17:32 2002
+   @author    Ian Hawke
+   @desc 
+   A header file with the definitions of the functions to be 
+   called by other thorns. This will exist until the function
+   aliasing is sorted out.
+   @enddesc 
+ @@*/
+
+#ifndef MOL_FUNCTIONS_H
+#define MOL_FUNCTIONS_H
+
+CCTK_INT MoL_RegisterEvolved(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex);
+CCTK_INT MoL_RegisterConstrained(CCTK_INT ConstrainedIndex);
+CCTK_INT MoL_RegisterSaveAndRestore(CCTK_INT SandRIndex);
+CCTK_INT MoL_ChangeToEvolved(CCTK_INT EvolvedIndex, CCTK_INT RHSIndex);
+CCTK_INT MoL_ChangeToConstrained(CCTK_INT ConstrainedIndex);
+CCTK_INT MoL_ChangeToSaveAndRestore(CCTK_INT SandRIndex);
+CCTK_INT MoL_ChangeToNone(CCTK_INT RemoveIndex);
+CCTK_INT MoL_RegisterEvolvedGroup(CCTK_INT EvolvedGroupIndex, 
+                                  CCTK_INT RHSGroupIndex);
+CCTK_INT MoL_RegisterConstrainedGroup(CCTK_INT ConstrainedGroupIndex);
+CCTK_INT MoL_RegisterSaveAndRestoreGroup(CCTK_INT SandRGroupIndex);
+
+#endif