diff options
| author | Aaron LI <aly@aaronly.me> | 2017-11-18 13:03:29 +0800 |
|---|---|---|
| committer | Aaron LI <aly@aaronly.me> | 2017-11-18 13:03:29 +0800 |
| commit | b7beb1a6d7abe192c7b97a18ec389b0270af0bc6 (patch) | |
| tree | 29b47fa15c5752930a72d32e6508045dcd3b9393 | |
| parent | b82305558b8c574092aaf6bfda515190918f85bc (diff) | |
| download | atoolbox-b7beb1a6d7abe192c7b97a18ec389b0270af0bc6.tar.bz2 | |
xspec_instlines_mc.tcl: Freeze line norms before Monte Carlo
| -rw-r--r-- | astro/xmm/xspec_instlines_mc.tcl | 4 |
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 |