Revision 1108 trunk/src/nh99/xmonte2.cpp

xmonte2.cpp (revision 1108)
2 2
 * $Id$
3 3
 *
4 4
 * Author: David Fournier
5
 * Copyright (c) 2008-2012 Regents of the University of California 
5
 * Copyright (c) 2008-2012 Regents of the University of California
6 6
 */
7 7
#include <admodel.h>
8 8

  
......
24 24
  a=a1;
25 25
  b=b1;
26 26
  wght=0;
27
  double ah; 
28
  double bl; 
27
  double ah;
28
  double bl;
29 29
  double upper;
30
  double lower; 
30
  double lower;
31 31
  double diff;
32 32
  int expflag;
33 33
  int in=0;
34 34
  int ie=0;
35 35
  for (int i=1;i<=nvar;i++)
36 36
  {
37
    ah=a(i)/ch(i,i); 
38
    bl=b(i)/ch(i,i); 
37
    ah=a(i)/ch(i,i);
38
    bl=b(i)/ch(i,i);
39 39
    upper=cumd_norm(bl);
40
    lower=cumd_norm(ah); 
40
    lower=cumd_norm(ah);
41 41
    diff=upper-lower;
42 42
    if (diff>1.e-5)
43 43
    {
......
89 89
  a=a1;
90 90
  b=b1;
91 91
  wght=0;
92
  double ah; 
93
  double bl; 
92
  double ah;
93
  double bl;
94 94
  double upper;
95
  double lower; 
95
  double lower;
96 96
  double diff;
97 97
  double diff1;
98 98
  //int in=0;
99 99
  //int ie=0;
100 100
  for (int i=1;i<=nvar;i++)
101 101
  {
102
    ah=a(i)/ch(i,i); 
103
    bl=b(i)/ch(i,i); 
102
    ah=a(i)/ch(i,i);
103
    bl=b(i)/ch(i,i);
104 104
    upper=cumd_norm(bl);
105
    lower=cumd_norm(ah); 
105
    lower=cumd_norm(ah);
106 106
    diff=upper-lower;
107 107
    upper=cumd_cauchy(bl);
108 108
    lower=cumd_cauchy(ah);
......
136 136
  a=a1;
137 137
  b=b1;
138 138
  wght=0;
139
  double ah; 
140
  double bl; 
139
  double ah;
140
  double bl;
141 141
  double upper;
142
  double lower; 
142
  double lower;
143 143
  double diff;
144 144
  for (int i=1;i<=nvar;i++)
145 145
  {
146
    ah=a(i)/ch(i,i); 
147
    bl=b(i)/ch(i,i); 
146
    ah=a(i)/ch(i,i);
147
    bl=b(i)/ch(i,i);
148 148
    lower=ffmax(-1.0,ah);
149 149
    upper=ffmin(1.0,bl);
150 150
    diff=upper-lower;

Also available in: Unified diff