Changeset 20936 for trunk/Ohana/src/opihi/cmd.data/test/vgauss.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/vgauss.sh (modified) (4 diffs)

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

@@ -20 +20 @@
 
  vgauss -q x y dy yfit
+ # lim x y; clear; box; plot -x 1 -c black x y; plot -c red x yfit
+ # cursor
 
  if (abs($C0 - 0.0) > 0.01)
@@ -51 +53 @@
 
  vgauss -q x y dy yfit
+ # lim x y; clear; box; plot -x 1 -c black x y; plot -c red x yfit
+ # cursor
 
  if (abs($C0 - 0.0) > 0.01)
@@ -71 +75 @@
  break -auto off
 
+ $C0o = 0.5
+ $C1o = 2
+ $C2o = 900
+ $C3o = 1
+
  create x -10 10 0.1
- set y = 1000 * exp(-0.5*x^2/3^2)
+ set y = $C3o + $C2o * exp(-0.5*(x - $C0o)^2/$C1o^2)
 
  gaussdev dY y[] 0.0 1.0
@@ -78 +87 @@
  set y = y + dy
 
- $C0 = 0.5
- $C1 = 2
- $C2 = 900
- $C3 = 1
+ $C0 = $C0o + 2
+ $C1 = $C1o + 2
+ $C2 = $C2o + 50
+ $C3 = $C3o + 10
 
  vgauss -q x y dy yfit
+ # lim x y; clear; box; plot -x 1 -c black x y; plot -c red x yfit
+ # cursor
 
- if (abs($C0 - 0.0) > 0.01)
+ $dS = 3.0 / sqrt(y[])
+
+ 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) > 1)
+ if (abs($C2 - $C2o)/$C2o > $dS)
    $PASS = 0
  end
- if (abs($C3 - 0.0) > 0.2)
+ if (abs($C3 - $C3o) > $dS)
    $PASS = 0
  end