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"}];