[Commits] [svn:einsteintoolkit] WeylScal4/trunk/m/ (Rev. 78)

[ian.hinder at aei.mpg.de]

User: hinder
Date: 2011/09/25 02:20 PM

Modified:
 /trunk/m/
  WeylScal4.m

Log:
 Use run-time finite difference order selection
 
 Keep the old method for backwards compatibility for now

File Changes:

Directory: /trunk/m/
====================

File [modified]: WeylScal4.m
Delta lines: +19 -2
===================================================================
--- trunk/m/WeylScal4.m	2011-09-25 19:19:52 UTC (rev 77)
+++ trunk/m/WeylScal4.m	2011-09-25 19:20:01 UTC (rev 78)
@@ -25,6 +25,11 @@
   PDstandard4th[i_, j_] -> StandardCenteredDifferenceOperator[1,2,i] *
     StandardCenteredDifferenceOperator[1,2,j],
 
+  PDstandard[i_] -> StandardCenteredDifferenceOperator[1,fdOrder/2,i],
+  PDstandard[i_, i_] -> StandardCenteredDifferenceOperator[2,fdOrder/2,i],
+  PDstandard[i_, j_] -> StandardCenteredDifferenceOperator[1,fdOrder/2,i] *
+    StandardCenteredDifferenceOperator[1,fdOrder/2,j],
+
   PDplus[i_] -> DPlus[i],
   PDminus[i_] -> DMinus[i],
   PDplus[i_,j_] -> DPlus[i] DPlus[j],
@@ -277,8 +282,8 @@
 fdOrderParam = 
 {
   Name -> "fd_order",
-  Default -> "2nd",
-  AllowedValues -> {"2nd", "4th"}
+  Default -> "Nth",
+  AllowedValues -> {"Nth", "2nd", "4th"}
 };
 
 keywordParameters = 
@@ -286,8 +291,19 @@
   fdOrderParam
 };
 
+intParameters =
+{
+  {
+    Name -> fdOrder,
+    Default -> 2,
+    AllowedValues -> {2,4,6,8}
+  }
+};
+
+
 calculations = 
 {
+  PsisCalc["Nth", PDstandard],
   PsisCalc["2nd", PDstandard2nd],
   PsisCalc["4th", PDstandard4th]
 };
@@ -299,4 +315,5 @@
   UseLoopControl -> True,
   KeywordParameters -> keywordParameters,
   RealParameters -> realParameters,
+  IntParameters -> intParameters,
   InheritedImplementations -> {"admbase", "methodoflines"}];