aboutsummaryrefslogtreecommitdiffstats
path: root/astro/xmm
diff options
context:
space:
mode:
authorAaron LI <aly@aaronly.me>2017-11-18 13:03:29 +0800
committerAaron LI <aly@aaronly.me>2017-11-18 13:03:29 +0800
commitb7beb1a6d7abe192c7b97a18ec389b0270af0bc6 (patch)
tree29b47fa15c5752930a72d32e6508045dcd3b9393 /astro/xmm
parentb82305558b8c574092aaf6bfda515190918f85bc (diff)
downloadatoolbox-b7beb1a6d7abe192c7b97a18ec389b0270af0bc6.tar.bz2
xspec_instlines_mc.tcl: Freeze line norms before Monte Carlo
Diffstat (limited to 'astro/xmm')
-rw-r--r--astro/xmm/xspec_instlines_mc.tcl4
1 files changed, 2 insertions, 2 deletions
diff --git a/astro/xmm/xspec_instlines_mc.tcl b/astro/xmm/xspec_instlines_mc.tcl
index 4a3d579..dc2d059 100644
--- a/astro/xmm/xspec_instlines_mc.tcl
+++ b/astro/xmm/xspec_instlines_mc.tcl
@@ -249,6 +249,7 @@ puts "Number of free parameters: [ llength $freeparameters ]"
print_header $freeparameters $outfd
puts "-----------------------------------------------------------------------"
+freeze_norms $gaussians
set tstart [ clock seconds ]
for {set i 0} {$i < $mc_times} {incr i} {
puts -nonewline "... [ expr {$i + 1} ] / ${mc_times} ..."
@@ -261,14 +262,13 @@ for {set i 0} {$i < $mc_times} {incr i} {
puts [ format " Elapsed %.1f min / ETA %.1f min ..." $elapsed $eta ]
}
randomize_norms $gaussians
- freeze_norms $gaussians
fit
print_data $freeparameters $outfd
}
puts "-----------------------------------------------------------------------"
set tnow [ clock seconds ]
set elapsed [ expr {($tnow - $tstart) / 60.0} ]
-puts [ format "Total elapsed time: %.1f min" $elapsed ]
+puts [ format "Total Monte Carlo time: %.1f min" $elapsed ]
close $outfd
# Recover query and chatter level