aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorIan Hinder <ian.hinder@aei.mpg.de>2013-09-12 15:20:58 +0200
committerIan Hinder <ian.hinder@aei.mpg.de>2013-09-12 15:20:58 +0200
commit7daca94f282a8b09a6ed0f431311b5c60a71150f (patch)
tree6ac6129fe1bef6bdcf9fbd55de50cd00ad131546
parentd7af35bceb859a45abc4d96805ceda90ca1874ac (diff)
KrancThorn.m: Reformat conservation calculation section
-rw-r--r--Tools/CodeGen/KrancThorn.m44
1 files changed, 21 insertions, 23 deletions
diff --git a/Tools/CodeGen/KrancThorn.m b/Tools/CodeGen/KrancThorn.m
index e2a0c98..588db40 100644
--- a/Tools/CodeGen/KrancThorn.m
+++ b/Tools/CodeGen/KrancThorn.m
@@ -178,29 +178,27 @@ CreateKrancThorn[groupsOrig_, parentDirectory_, thornName_, opts:OptionsPattern[
------------------------------------------------------------------------ *)
Module[
- {consCalcsIn, consCalcs, consGroups},
-
- consCalcsIn = Append[#,Groups -> groups]& /@
- OptionValue[ConservationCalculations];
-
- (* Add in calculations to solve any conservation laws that have
- been provided *)
-
- consCalcs = Flatten[Map[ProcessConservationCalculation[#,thornName] &,
- consCalcsIn],1];
-
- consCalcs = Map[Join[#, {PartialDerivatives -> partialDerivs,
- Implementation -> implementation}] &, consCalcs];
-
- calcs = Join[calcs,consCalcs];
- (* Print["consCalcs = ", consCalcs]; *)
-
- consGroups = Union@Flatten[
- Map[ConservationCalculationDeclaredGroups, consCalcsIn],1];
-
- groups = Join[groups, consGroups];
-
- declaredGroups = Join[declaredGroups, Map[groupName, consGroups]]];
+ {inputConsCalcs, outputConsCalcs, consGroups},
+
+ inputConsCalcs = Map[Append[#, Groups -> groups] &, OptionValue[ConservationCalculations]];
+
+ outputConsCalcs =
+ Flatten[
+ Map[
+ ProcessConservationCalculation[#, thornName] &,
+ inputConsCalcs],
+ 1];
+
+ outputConsCalcs = Map[Join[#, {PartialDerivatives -> partialDerivs,
+ Implementation -> implementation}] &,
+ outputConsCalcs];
+
+ consGroups = Union@Flatten[
+ Map[ConservationCalculationDeclaredGroups, inputConsCalcs],1];
+
+ calcs = Join[calcs,outputConsCalcs];
+ groups = Join[groups, consGroups];
+ declaredGroups = Join[declaredGroups, Map[groupName, consGroups]]];
(* ------------------------------------------------------------------------
ODEs