diff options
author | Ian Hinder <ian.hinder@aei.mpg.de> | 2013-09-12 15:20:58 +0200 |
---|---|---|
committer | Ian Hinder <ian.hinder@aei.mpg.de> | 2013-09-12 15:20:58 +0200 |
commit | 7daca94f282a8b09a6ed0f431311b5c60a71150f (patch) | |
tree | 6ac6129fe1bef6bdcf9fbd55de50cd00ad131546 | |
parent | d7af35bceb859a45abc4d96805ceda90ca1874ac (diff) |
KrancThorn.m: Reformat conservation calculation section
-rw-r--r-- | Tools/CodeGen/KrancThorn.m | 44 |
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 |