From 6c574f5f886e0614a415593e2e12494e21469c6c Mon Sep 17 00:00:00 2001
From: Eugeniy Mikhailov <evgmik@gmail.com>
Date: Tue, 27 Nov 2012 14:15:55 -0500
Subject: move Maxwell distribution from computed vector to just a vector
MIME-Version: 1.0
Content-Type: text/plain; charset=UTF-8
Content-Transfer-Encoding: 8bit

Thanks to suggestion of Graham Dennis, this skip unnecessary calculations
of the Maxwell probability distributions at every step since. Thus
avoids recalculation.

Below the original comment:
Your computed_vector's "Maxwell_distribution_probabilities" and
"Maxwell_distribution_probabilities_norm" are constant.  If you just
change them to be <vector> elements changing the <evaluation> tags to
<initialisation> tags, this will ensure that they are not evaluated
every time they are needed.
---
 xmds2/realistic_Rb/realistic_Rb.xmds | 16 ++++++++--------
 1 file changed, 8 insertions(+), 8 deletions(-)

diff --git a/xmds2/realistic_Rb/realistic_Rb.xmds b/xmds2/realistic_Rb/realistic_Rb.xmds
index 94b398c..57e7425 100644
--- a/xmds2/realistic_Rb/realistic_Rb.xmds
+++ b/xmds2/realistic_Rb/realistic_Rb.xmds
@@ -185,29 +185,29 @@ const double rt2 = 1.4142135623730951;
 	</vector>
 
 	<!--Maxwell distribution probability p(v)-->
-	<computed_vector name="Maxwell_distribution_probabilities" dimensions="v" type="real">
+	<vector name="Maxwell_distribution_probabilities" dimensions="v" type="real">
 		<components>probability_v</components>
-		<evaluation>
+		<initialisation>
 			<![CDATA[
 			// TODO: move to the global space/function. This reevaluated many times since it called from dependency requests but it never changes during  the script lifetime since 'v' is fixed.
 			probability_v=1.0/(v_thermal_averaged*sqrt(2*M_PI)) * exp( - mod2(v/v_thermal_averaged)/2.0 ); 
 			]]>
-		</evaluation>
-	</computed_vector>
+		</initialisation>
+	</vector>
 
 	<!--Maxwell distribution norm sum(p(v))
 			 Needed since we sum over the grid instead of true integral,
 			 we also have finite cut off velocities-->
-	<computed_vector name="Maxwell_distribution_probabilities_norm" dimensions="" type="real">
+	<vector name="Maxwell_distribution_probabilities_norm" dimensions="" type="real">
 		<components>probability_v_norm</components>
-		<evaluation>
+		<initialisation>
 			<dependencies basis="v">Maxwell_distribution_probabilities</dependencies>
 			<![CDATA[
 			// TODO: move to the global space/function. This reevaluated many times since it called from dependency requests but it never changes during  the script lifetime since 'v' is fixed.
 			probability_v_norm=probability_v;
 			]]>
-		</evaluation>
-	</computed_vector>
+		</initialisation>
+	</vector>
 
 
 	<!-- Averaged across Maxwell distribution fields amplitudes -->
-- 
cgit v1.2.3