Changeset 20936 for trunk/Ohana/src/opihi/cmd.data/test/vmaxwell.sh

Timestamp:

Dec 7, 2008, 3:31:01 PM (18 years ago)

Author:

eugene

Message:

big update from eam_branch_20081124 with updates to Opihi math

File:

• trunk/Ohana/src/opihi/cmd.data/test/vmaxwell.sh (modified) (1 diff)

trunk/Ohana/src/opihi/cmd.data/test/vmaxwell.sh

@@ -40 +40 @@
  break -auto off
 
+ $C0o = 250
+ $C1o = 25
+ $C2o = 100
+ $C3o = 1
+ $C4o = 10
+
  create x 0 1000 1
- set y = 1000 * (x-200)^2 * exp(-0.5*(x-300)^2/30^2)
- set dy = (rnd(y) - 0.5)/0.5
+ set y = $C2o * (x-$C4o)^2 * exp(-0.5*(x-$C0o)^2/$C1o^2) + $C3o
+ gaussdev dY y[] 0.0 1.0
+ set dy = dY * sqrt(y)
  set y = y + dy
 
- lim x y; clear; box; plot -x 1 -c black x y
+ $C0 = $C0o + 20
+ $C1 = $C1o + 2
+ $C2 = $C2o + 100
+ $C3 = $C3o + 10
+ $C4 = $C4o + 20
 
- $C0 = 250
- $C1 = 25
- $C2 = 900
- $C3 = 1
- $C4 = 190
+ vmaxwell -q x y dy yfit
+ lim x y; clear; box; plot -x 1 -c black x y; plot x yfit -c red
 
- vmaxwell x y dy yfit
- plot x yfit -c red
+ $dS = 6.0 / sqrt(y[])
 
- if (abs($C0 - 0.0) > 0.01)
+ if (abs($C0 - $C0o) > $dS)
    $PASS = 0
  end
- if (abs($C1 - 3.0) > 0.01)
+ if (abs($C1 - $C1o)/$C1o > $dS)
    $PASS = 0
  end
- if (abs($C2 - 1000.0) > 0.1)
+ if (abs($C2 - $C2o)/$C2o > $dS)
    $PASS = 0
  end
- if (abs($C3 - 0.0) > 0.1)
+ if (abs($C3 - $C3o) > $dS)
+   $PASS = 0
+ end
+ if (abs($C4 - $C4o)/$C4o > $dS)
    $PASS = 0
  end