예제 #1
0
int main(void)
{
	float a=0.0,b=PIO2,s;

	printf("Integral of func computed with QROMB\n\n");
	printf("Actual value of integral is %12.6f\n",fint(b)-fint(a));
	s=qromb(func,a,b);
	printf("Result from routine QROMB is %11.6f\n",s);
	return 0;
}
예제 #2
0
void fint(int L,int R,POINT v)
{
    if(L>R)return;
    int M=(L+R)>>1;
    ll tmp=tree[M].dis(v);
    if(tmp!=0) ans=min(ans,tmp);
    if(flag[M]) {
        ll d=(v.x-tree[M].x)*(v.x-tree[M].x);
        if(v.x<=tree[M].x) {
            fint(L,M-1,v);
            if(d<ans) fint(M+1,R,v);
        } else {
            fint(M+1,R,v);
            if(d<ans) fint(L,M-1,v);
        }
    } else {
        ll d=(v.y-tree[M].y)*(v.y-tree[M].y);
        if(v.y<=tree[M].y) {
            fint(L,M-1,v);
            if(d<ans) fint(M+1,R,v);
        } else {
            fint(M+1,R,v);
            if(d<ans) fint(L,M-1,v);
        }
    }
}
예제 #3
0
int main(void)
{
	int i;
	float a=0.0,b=PIO2,s;

	printf("\nIntegral of func with 2^(n-1) points\n");
	printf("Actual value of integral is %10.6f\n",fint(b)-fint(a));
	printf("%6s %24s\n","n","approx. integral");
	for (i=1;i<=NMAX;i++) {
		s=trapzd(func,a,b,i);
		printf("%6d %20.6f\n",i,s);
	}
	return 0;
}
예제 #4
0
파일: xmidpnt.c 프로젝트: bamford/astrobamf
int main(void)
{
	float a=0.0,b=1.0,s;
	int i;

	printf("\nIntegral of func computed with MIDPNT\n");
	printf("Actual value of integral is %7.4f\n",(fint(b)-fint(a)));
	printf("%6s %29s \n","n","Approx. integral");
	for (i=1;i<=NMAX;i++) {
		s=midpnt(func,a,b,i);
		printf("%6d %24.6f\n",i,s);
	}
	return 0;
}
예제 #5
0
void MainWindow::on_pushButton_clicked()
{
    QString qstr=ui->comboBox->currentText();
    ui->lineEdit->setText(qstr);

    //qDebug()<<qstr;
    //const char *psearch=qstr.toStdString().c_str();//转化为C风格的字符串
    //qDebug()<<psearch;  // 字符串严格转换
    QByteArray ba=qstr.toLocal8Bit();
    const char *psearch=ba.data();

    ifstream fint("D:\\dangdangwang.txt");  //路径不能有中文
    if(!fint)
    {
        qDebug()<<"open fail";

    }
    else
    {
        qDebug()<<"open success";
    }
    int i=0;
    ;
     ;
    while(!fint.eof())
    {
        i++;
        if(i>4000)
        {
            break;
        }
        char strall[1024]={0};
        fint.getline(strall,1023);
       // qDebug()<<strall;  //QT读文件乱码

        QTextCodec *codec = QTextCodec::codecForName("gb2312");
       QString qstrall=codec->toUnicode(strall);  //中文

      //  qDebug()<<qstrall;
        const char *pstrall=qstrall.toStdString().c_str();//转化为C风格的字符串
        qDebug()<<pstrall;  // 字符串严格转换   //这种形式的转换进行比较会使得数据缺失


        char *pfind=strstr(strall,psearch);
        if(pfind!=NULL)
        {
            qDebug()<<pstrall;
            sl->append(qstrall);//添加到列表
            slm->setStringList(*sl);
            ui->listView->setModel(slm);
        }
        else
        {
            qDebug()<<"not find";
        }

    }


}
예제 #6
0
void Saturate(float *pb, int numSamples) {
	int		 i;
	double	 fr,s1,s2;

	//int oldctrlword = _control87(0, 0);
	//_control87(_RC_DOWN, _MCW_RC);			

	do {
		double s = *pb;
		s *= 127.0/(3.0*SatAmp);
		i = fint(s);
		if(i>127) 
			*pb++ = SatAmp;
		else if(i<-127)
			*pb++ = -SatAmp;
		else {
			fr = s - i; i+=127;
			s1 = Tanh[i++];
			s2 = Tanh[i];
			*pb++ = (float)((s1+(s2-s1)*fr)*SatAmp);
		}
	} while(--numSamples);

	//_control87(oldctrlword, _MCW_RC);
}
예제 #7
0
/*
 * remove rightest zeroes from the number
 */
var_num_t rmzeros(var_num_t n) {
  var_int_t i = fint(n);
  while (i > 0 && i % 10 == 0) {
    i /= 10;
  }
  return i;
}
예제 #8
0
 void add(int key1, int key2, long long val) {
   int tmp = fint(key1, key2);
   if (tmp == -1) {
     insert(key1, key2);
     ele[N - 1].val += val;
   } else
     ele[tmp].val += val;
 }
예제 #9
0
 long long get(int key1, int key2) {
   int tmp = fint(key1, key2);
   if (tmp == -1) {
     insert(key1, key2);
     return ele[N - 1].val;
   } else
     return ele[tmp].val;
 }
예제 #10
0
void MonitorSampledBar::draw_long_term(bool mustDraw) {
  if (!SpyShowLongTerm) return;

  // calculate long-term average (use floats to avoid overflow)
  if (elapsedT == 0) return;
  fint lta1 = fint((w-2) * ((float)sum1T / elapsedT) + 1);
  fint lta2 = fint((w-2) * ((float)sum2T / elapsedT) + 1);
  fint lta3 = fint((w-2) * ((float)sum3T / elapsedT) + 1);
  fint lta4 = fint((w-2) * ((float)sum4T / elapsedT) + 1);

  const fint markerH = 3;
  fint ys = y + (h - markerH) / 2;
  TheSpy->mw()->pw->set_xor();
  if (mustDraw || lta1 != old_lta1) {
    if (!mustDraw) TheSpy->mw()->pw->draw_line(x+old_lta1, ys, x+old_lta1, ys+markerH);
    TheSpy->mw()->pw->draw_line(x + lta1, ys, x + lta1, ys+markerH);
    old_lta1 = lta1;
  }
  // draw 2nd bar only if not at left end, and if not at same position as 1st
  if (mustDraw || lta2 != old_lta2) {
    if (!mustDraw && old_lta2 != 1)
      TheSpy->mw()->pw->draw_line(x+old_lta2, ys, x+old_lta2, ys+markerH);
    if (lta2 != 1 && lta2 != lta1)
      TheSpy->mw()->pw->draw_line(x + lta2, ys, x + lta2, ys+markerH);
    old_lta2 = lta2;
  }
  if (mustDraw || lta3 != old_lta3) {
    if (!mustDraw && old_lta3 != 1)
      TheSpy->mw()->pw->draw_line(x+old_lta3, ys, x+old_lta3, ys+markerH);
    if (lta3 != 1 && lta3 != lta2 && lta3 != lta1)
      TheSpy->mw()->pw->draw_line(x + lta3, ys, x + lta3, ys+markerH);
    old_lta3 = lta3;
  }
  if (mustDraw || lta4 != old_lta4) {
    if (!mustDraw && old_lta4 != 1)
      TheSpy->mw()->pw->draw_line(x+old_lta4, ys, x+old_lta4, ys+markerH);
    if (lta4 != 1 && lta4 != lta3 && lta4 != lta2 && lta4 != lta1)
      TheSpy->mw()->pw->draw_line(x + lta4, ys, x + lta4, ys+markerH);
    old_lta4 = lta4;
  }
  TheSpy->mw()->pw->set_copy();
}
예제 #11
0
파일: axlcomp.c 프로젝트: nilqed/aldor
void
compGLoopEval(FILE * fin, FILE * fout, EmitInfo finfo)
{
        Stab            stab = stabFile();
	AbSyn		ab;
	Foam		foam;
        static int      lineno = 0;
      
        comsgInit();
      
        while (!osFEof(fin)) {
			   Bool ok;
               intStepNo += 1;
               comsgFini();
               comsgInit();
               car(stab)->isChecked = false;

               if (SetJmp(compFintJmpBuf))
                        if (osFEof(fin))
                                break;

               comsgPromptPrint(fin, fout,
                                fintHistory ? "%%%d := " : "%%%d >> ",
                                intStepNo);
               osSetBreakHandler(compFintBreakHandler0);
               ab = compFileFront(finfo, stab, fin, &lineno);
               breakSetRoot(ab);

               if (!compIsMoreAfterFront(finfo) ||
                   abIsEmptySequence(ab))
                      continue;
                
               ab = (AbSyn) fintWrap(ab, intStepNo);
               foam = compFileMiddle(finfo, stab, ab);

               if (!foam)   continue;

               ok = fint(foam);
	       if (!ok)
		       fprintf(fout, "Unhandled Exception!\n");
               foamFree(foam);

               if (fintVerbose && ok)
                      fintPrintType(fout, ab);
               fintDisplayTimings();
        }
      
        comsgFini();                
}
예제 #12
0
int main() {
  int x[DIM1][DIM2];
  int i, j;
  for (i = 0; i < DIM1; i++)
    for(j = 0; j < DIM2; j++)
      x[i][j] = i+j*i;
  globalBoo.a = 42;
  globalBoo.b = 9876;
  globalBoo.c[0][0] = 12345;
  globalBoo.c[1][0] = 54321;
  globalBoo.c[2][0] = 42;
  globalBoo.c[9][49] = 949;
  fint(x);
  fcrashint(x, globalBoo);
  printf("sizeof(struct boo): %d\n", sizeof(struct boo));
  return 0;
}
예제 #13
0
int main()
{
    int T;
    scanf("%d",&T);
    while(T--) {
        int n;
        scanf("%d",&n);
        for(int i=0; i<n; i++) {
            scanf("%lld%lld",&p[i].x,&p[i].y);
            q[i]=p[i];
        }
        built(0,n-1);
        for(int i=0; i<n; i++) {
            ans=INF;
            fint(0,n-1,q[i]);
            printf("%I64d\n",ans);
        }
    }
    return 0;
}
예제 #14
0
enum judge_status check_user_output()
{
    struct stat st,ss;
    string strt,strs;
    stat(code_output_path,&st);
    stat(input_path,&ss);

    if(st.st_size < ss.st_size / 2 || st.st_size > ss.st_size * 2) return JUDGE_WA; //accelerate
    bool flag = false;
    ifstream fint(code_output_path,ifstream::in), fins(input_path,ifstream::in);


    while(fint.good() || fins.good())
    {
        if(fint.good())
            getline(fint,strt);
        else strt = "";

        if(fins.good())
            getline(fins,strs);
        else strs = "";

        if(strt != strs)
        {
            trim(strt);
            trim(strs);
            if(strt == strs)
                flag = true;
            else
                return JUDGE_WA;
        }
    }
    if(flag)
        return JUDGE_PE;
    return JUDGE_AC;
}
예제 #15
0
파일: floatOop.cpp 프로젝트: AdamSpitz/self
   Self_float = { int sign:1; int fract: 21; int exp: 8; int tag: 2}
*/
  
static const fint fractSize = 23;
static const fint expSize = 8;
static const fint expOffset = fractSize;

# ifdef FAST_FLOATS
  floatOop infinityOop = floatOop(nthMask(expSize) << expOffset
                                  |  Float_Tag);
# else

  static const fint signSize = 1;
  static const fint signOffset = fractSize + expSize;
  
  static const fint selfExpSize   = fint(expSize   - Tag_Size);
  static const fint selfExpOffset = fint(expOffset + Tag_Size);
  
  static const fint bias = nthBit(expSize) / 2  -  1;
  static const fint selfBias = nthBit(selfExpSize) / 2  -  1;
  
  floatOop infinityOop = floatOop(nthMask(selfExpSize) << selfExpOffset
                                  |  Float_Tag);
  
  floatOop as_floatOop(float value) {
    union { float f; uint32 i; } x;
    x.f = value;
    uint32 i = x.i;
    int32 exp     = i >> expOffset  &  nthMask(expSize);
    int32 selfExp = exp - bias + selfBias;
    int32 fract   = i & nthMask(fractSize);
예제 #16
0
/*
 * best float to string (lib)
 *
 * This is the real float-to-string routine.
 * It used by the routines:
 *   bestfta(double x, char *dest)
 *   expfta(double x, char *dest)
 */
void bestfta_p(var_num_t x, char *dest, var_num_t minx, var_num_t maxx) {
  var_num_t ipart, fpart, fdif;
  var_int_t power = 0;
  int sign, i;
  char *d = dest;
  char buf[64];

  memset(buf, 0, sizeof(buf));

  if (fabs(x) == 0.0) {
    strcpy(dest, "0");
    return;
  }

  // find sign
  sign = sgn(x);
  if (sign < 0) {
    *d++ = '-';
  }
  x = fabs(x);

  if (x >= 1E308) {
    *d = '\0';
    strcat(d, WORD_INF);
    return;
  } else if (x <= 1E-307) {
    *d = '\0';
    strcat(d, "0");
    return;
  }

  // find power
  if (x < minx) {
    for (i = 37; i >= 0; i--) {
      if (x < nfta_eminus[i]) {
        x *= nfta_eplus[i];
        power = -((i + 1) * 8);
      } else {
        break;
      }
    }

    while (x < 1.0 && power > -307) {
      x *= 10.0;
      power--;
    }
  } else if (x > maxx) {
    for (i = 37; i >= 0; i--) {
      if (x > nfta_eplus[i]) {
        x /= nfta_eplus[i];
        power = ((i + 1) * 8);
      } else {
        break;
      }
    }

    while (x >= 10.0 && power < 308) {
      x /= 10.0;
      power++;
    }
  }

  // format left part
  ipart = fabs(fint(x));
  fpart = fround(frac(x), FMT_RND) * FMT_xRND;
  if (fpart >= FMT_xRND) {      // rounding bug
    ipart = ipart + 1.0;
    if (ipart >= maxx) {
      ipart = ipart / 10.0;
      power++;
    }
    fpart = 0.0;
  }

  fptoa(ipart, buf);
  strcpy(d, buf);
  d += strlen(buf);

  if (fpart > 0.0) {
    // format right part
    *d++ = '.';

    fdif = frac(x) * FMT_xRND;
    if (fdif < fpart) {
      // rounded value has greater precision
      fdif = fpart;
    }

    while (fdif < FMT_xRND2) {
      fdif *= 10;
      *d++ = '0';
    }

    fptoa(rmzeros(fpart), buf);
    strcpy(d, buf);
    d += strlen(buf);
  }

  if (power) {
    // add the power
    *d++ = 'E';
    if (power > 0) {
      *d++ = '+';
    }
    fptoa(power, buf);
    strcpy(d, buf);
    d += strlen(buf);
  }

  // finish
  *d = '\0';
}
예제 #17
0
파일: axlcomp.c 프로젝트: nilqed/aldor
int
compInteractiveLoop(int argc, char **argv, FILE *fin, FILE *fout)
{
	int		iargc, totErrors, lineno;
	FileName	fn;
	EmitInfo	finfo;
	AbSyn		ab;
	Stab		stab;
	Foam		foam;
	Bool		readingInitFile = true, tmpHistory;
	Bool		endOfInput = false; 
	FILE		*fin0 = fin;

	compInit();

	iargc = cmdArguments(1, argc, argv);
 
	argc	 -= iargc;
	argv	 += iargc;
	emitDoneOptions(argc, argv);


	if (osFileIsThere(cmdInitFile)) {
		FileName	fname;

		fname = fnameParse(cmdInitFile);
		fin = fileMustOpen(fname, osIoRdMode);
		fnameFree(fname);
		intStepNo = 0;
		tmpHistory = fintHistory;
		fintHistory = false;
		fprintf(fout, "Reading %s...\n", cmdInitFile);
	}
	else {
		intStepNo = 1;
		readingInitFile = false;
		tmpHistory = fintHistory;
	}

	fn  = fnameStdin();

	finfo  = emitInfoNew(fn);
 
	lineno	  = 0;
 
	compFileInit(finfo);
	stab = stabFile();
	comsgFini();
 
	fintInit();

	for (; !endOfInput; intStepNo++) {
		comsgInit();
		breakSetRoot(NULL);
		car(stab)->isChecked = false;

		if (SetJmp(compFintJmpBuf)) {
			if (feof(fin)) break;
		}

		comsgPromptPrint(fin, fout, 
		       fintHistory ? "%%%d := " : "%%%d >> ",
		       intStepNo);
 
		osSetBreakHandler(compFintBreakHandler0);

		ab = compFileFront(finfo, stab, fin, &lineno);
		breakSetRoot(ab);

		if (compIsMoreAfterFront(finfo) &&
		    !abIsEmptySequence(ab)) {
			ab = (AbSyn) fintWrap(ab, intStepNo);
 
			foam = compFileMiddle(finfo, stab, ab);
			if (foam) {
 				Bool ok = fint(foam);
				if (ok && fintVerbose)
					fintPrintType(fout, ab);

				foamFree(foam);
				fintDisplayTimings();
			}
		}

		comsgFini();
		/* abFree(ab); !! ab is seeping into types. */

		if (feof(fin)) {
			if (readingInitFile) {
				fclose(fin);

				fin = fin0;
				lineno	  = 0;
				readingInitFile = false;
				fintHistory = tmpHistory;
				intStepNo = 0;
 
				comsgFini();
			}
			else
				endOfInput = true;
		}
	}			
 
	totErrors = comsgErrorCount();

	fintFini();
 
	compFileFini(finfo);
	emitAllDone();
	emitInfoFree(finfo);
	fnameFree(fn);
	compFini();

	return totErrors;
}
예제 #18
0
        GwtVantageFxBundle::GwtVantageFxBundle() {

            add("com.optionfair.client.model.ClientUpdatePack/839029662", {
                fstd("some_1"),
                fptr("assetUpdates"),
                fptr("money"),
                fptr("some_3"),
                fstd("instrumentType"),
                fstd("some_5"),
                fstd("some_6"),
                fstd("some_7"),
                fptr("options"),
                fstd("some_10"),
                fptr("some_11"),
                fptr("optionUpdates"),
                fstd("some_12"),
                fstd("some_13"),
                fptr("some_14"),
                fstd("some_15"),
                fptr("positionsSentimentDto"),
                fptr("positionUpdates"),
                fint("some_17"),
                fint("some_18"),
                fint("some_19"),
                fptr("some_20"),
                flng("date"),
                fptr("some_22"),
                fint("some_23"),
                fint("some_24")
            });

            add("com.optionfair.client.model.PositionsSentimentDTO/3338833266", {
                fptr("map"),
                fstd("some_1"),
                fstd("some_2")
            });

            add("com.optionfair.client.model.FeedTickUpdate/4036343884", {
                fdbl("some_1"),
                fptr("date"),
                fstd("some_2"),
                fstd("some_3")
            });

            add("com.optionfair.client.model.Option/1343902663", {
				fstr("serialuid"),
                fint("assetId"), // Для какой пары текущий опцион
                fstd("some_2"),
                fint("optionSeconds"),
                fstd("some_4"),
				fint("some_5"),
				fptr("closeDate"),
                fptr("some_7"),
                fint("expiryTypeId"),
                flng("id"),
                fdbl("some_10"),
                fint("instrumentId"),
                fstd("some_12"),
                fptr("some_13"),
                fdbl("some_14"),
                fdbl("some_15"),
                flng("some_16"),
                fint("some_17"),
                fint("some_18"),
                fint("some_19"),
                flng("some_20"),
				fstd("some_21"),
				fstd("some_22"),
				fint("some_23"),
				fint("return"),
				fstr("some_24"),
                fdbl("price"),
                fstd("some_26"),
                fptr("openFrom"),
				fint("some_int"),
				fptr("some_date"),
				fint("optionStatus"),
				fdbl("some_29"),
                fdbl("some_30"),
                fdbl("some_31"),
                fdbl("some_32"),
                flng("some_33"),
            });


/**
0
124
bVQI
7.679500e-01
311
0
6 [java.util.Date/3385151746]
Vo7NMsA
0.000000e+00
0.000000e+00
HUw
5 [com.optionfair.client.model.Option/1343902663]
1 [com.optionfair.client.model.ClientUpdatePack/839029662]
1 [com.optionfair.client.model.ClientUpdatePack/839029662]
ZN
0
0
0
170
7 [1]
7.679500e-01
0
0
0
6 [java.util.Date/3385151746]
Vo7NMsA
14
7.679500e-01
0.000000e+00
0.000000e+00
0.000000e+00
O
5 [com.optionfair.client.model.Option/1343902663]
5 [com.optionfair.client.model.Option/1343902663]
0
60
0
-1
6 [java.util.Date/3385151746]
Vo$o7MA
0
124
bVP2
1.252100e+00
311
0
6 [java.util.Date/3385151746]
Vo7NMsA
0.000000e+00
0.000000e+00
Opg
1 [com.optionfair.client.model.ClientUpdatePack/839029662]
1 [com.optionfair.client.model.ClientUpdatePack/839029662]
1 [com.optionfair.client.model.ClientUpdatePack/839029662]
P_
0
0
0
170
7 [1]
1.252100e+00
0
0
0
6 [java.util.Date/3385151746]
Vo7NMsA
14
1.252100e+00
0.000000e+00
0.000000e+00
0.000000e+00
O
5 [com.optionfair.client.model.Option/1343902663]
14
0
60
0
-1
6 [java.util.Date/3385151746]
Vo9WhwA
0
124
bVTl
5.340500e+01
311
0
6 [java.util.Date/3385151746]
Vo7owXA
0.000000e+00
0.000000e+00
HUw
1 [com.optionfair.client.model.ClientUpdatePack/839029662]
1 [com.optionfair.client.model.ClientUpdatePack/839029662]
1 [com.optionfair.client.model.ClientUpdatePack/839029662]
Tr
0
0
0
165
7 [1]
5.340500e+01
0
0
0
6 [java.util.Date/3385151746]
Vo7oqgA
14
5.340500e+01
0.000000e+00
0.000000e+00
0.000000e+00
O
5 [com.optionfair.client.model.Option/1343902663]
4 [java.lang.Long/4227064769]
0
300
0
-1
6 [java.util.Date/3385151746] */

			add("com.optionfair.client.model.AssetCurrencyRel/2204523642", {
				fstr("serialuid"), //6[rO0ABXcEAAAAAA == ]
				fstd("some1"), //43[Wagered]
                fptr("pk"),
				fstd("std1"), //1[com.optionfair.client.model.LutResolverInitPackage / 2546929785]
				fstd("std2"), //1[com.optionfair.client.model.LutResolverInitPackage / 2546929785]
			});

            add("com.optionfair.client.model.AssetCurrencyRelPK/4186144263", {
                fint("some_id"), //43[Wagered]
            });

            add("com.optionfair.client.model.PositionUpdate/3887973789", {
                fdbl("closedPrice"),//0.000000e+00
                fptr("positionCloseDate"),//0
                flng("transactionId"),//KlvU
                fint("instrumentId"),//311
                flng("date"),//VRdvf30
                flng("returned"),//A
                flng("some_7"),//A
                fint("positionStatus"),//18
                fdbl("some_9"),//0.000000e+00
            });

            add("com.optionfair.client.model.AssetUpdate/3263927368", {
				fstr("serialuid"), //6[rO0ABXcEAAAAAA == ]
				fint("assetStatusId"), // статус лота (нужен Production)
                fint("assetId"), // Идентификатор лота
                fdbl("targetPrice"), // Текущая цена
            });

            add("com.optionfair.client.model.OptionUpdateChangeStrike/3365438592", {
                fdbl("targetPrice"),
                fdbl("some_2"),
                flng("optionId"),
                fdbl("some_4"),
                fdbl("some_5"),
                fint("optionStatus"),
            });

			add("com.optionfair.client.resolver.BrandCurrencyCountryResolver/4149870428", {
                fptr("countries"),
                fptr("currencies"),
                fptr("currencyMap"),
                fstd("some_1"),
			});

            add("com.optionfair.client.model.LutResolverInitPackage/2546929785", {
                fptr("pspToPspType"),
                fptr("lutTypes"),
                fptr("paymentMethods"),
                fptr("assets"),
                fptr("someToPspCcTypeRels"),
                fptr("some_3"),
                fptr("blackWhiteList"),
                fptr("someList1"),
                fptr("brand"),
                fptr("brandCurrencyCountryResolver"),
                fstd("externalId"),
                fstd("some_5"),
                fstr("enfinium"),
                fptr("breandPspRels"),
                fptr("some_8"),
                fptr("brandRegulations"),
                fptr("primaryCurrencies"),
                fptr("limits"),
                fstd("basic"),
                fptr("date"),
                fptr("primaryDateFormat"),
                fstd("some_10"),
                fstd("some_11"),
                fptr("instrumentTypeFilters"),
                fptr("instrumentTypeSuperToInstrumentTypes"),
                fptr("instrumentTypeToInstrument"),
                fptr("superRels"),
                fstd("some_15"),
                fstd("some_16"),
                fstd("num_99999"),
                fptr("LoginDTO"),
                fstd("some_19"),
                fptr("markets"),
                fptr("messages"),
                fstd("some_20"),
                fptr("servers"),
                fptr("acceptedMethods"),
                fptr("pspCcTypeRels"),
                fstd("some_21"),
                fptr("some_22"),
                fptr("states"),
                fptr("some_23"),
                fptr("trackingInfo"),
                fptr("verifications"),
                fstd("some_24"),
                fstd("some_25")
            });

            add("com.optionfair.client.model.TradeLimits/516060049", {
                flng("some_1"),
                flng("some_2"),
                flng("some_3")
            });

            add("com.optionfair.client.ui.pages.tradePage.spread.High/1043352957", {
                fstr("some_1"),
                fdbl("some_2")
            });

			add("com.optionfair.client.ui.pages.tradePage.spread.NoSpread/2332320922", {
				//fint("some_1"),
				//fint("some_2"),
			});

            add("com.optionfair.client.ui.pages.tradePage.spread.Low/1807028581", {
                fstr("some_1"),
                fdbl("some_2")
            });

            add("com.optionfair.client.ui.pages.tradePage.spread.Touch/1022019549", {
                fstr("some_1"),
                fdbl("some_2")
            });

            add("com.optionfair.client.ui.pages.tradePage.spread.NoTouch/1845245280", {
                fstr("some_1"),
                fdbl("some_2")
            });

            add("com.optionfair.client.ui.pages.tradePage.spread.In/3813998762", {
                fstr("some_1"),
                fdbl("some_2")
            });

            add("com.optionfair.client.ui.pages.tradePage.spread.Out/2494560053", {
                fstr("some_1"),
                fdbl("some_2")
            });

            add("com.optionfair.client.ui.pages.tradePage.spread.Probability/1287790398", {
                fstr("some_1"),
                fdbl("some_2"),
            });

            add("com.optionfair.client.model.LoginDTO/3023568058", {
                fptr("account"),
                fstr("login"),
                fstr("some_1"),
                fptr("someMap"),
                fstd("some_2"),
                fstd("some_3"),
                fdbl("some_4"),
                fint("some_5"),
                fstd("some_6"),
                fptr("tradeLimits"),
                fstd("some_8"),
                fstd("some_9")
            });

            add("com.optionfair.client.spread.SpreadCalculator/3758543021", {
                flng("some_1"),
                fdbl("some_2"),
                fdbl("some_3"),
                fdbl("some_4"),
                fdbl("some_5"),
                fptr("spreads"),
                fstd("some_6"),
                fstd("some_7"),
                fdbl("some_8"),
                fdbl("some_9"),
                fstd("some_10"),
                fint("some_11"),
                fstd("some_12"),
                fstd("some_13"),
                fstd("some_14"),
                flng("some_15"),
                fstd("some_16"),
            });



            add("com.optionfair.client.model.Account/3668294554", {
                fstr("serialuid"),
                fstd("some_1"),
                fstr("regular"),
                fint("some_2"),
                fstd("some_3"),
                flng("some_4"),
                fstd("some_5"),
                fstd("some_6"),
                fstd("some_7"),
                fstd("some_8"),
                fstd("some_9"),
                flng("some_10"),
                flng("some_11"),
                fptr("some_date"),
                fptr("someList_1"),
                fptr("someList_2"),
                fdbl("some_12"),
                fint("externalId"),
                fptr("integer_1"),
                fstd("some_13"),
                fstd("some_14"),
                fint("some_15"),
                fstr("gotham"),
                fstr("guid"),
                fstr("name"),
                fstr("how"),
                fstd("some_16"),
                fstd("some_17"),
                fstr("email"),
                fstd("some_18"),
                fstr("login"),
                fstd("some_19"),
                fstr("fullName"),
                flng("some_20"),
                fstd("some_21"),
                flng("accountId"),
                fstd("some_23"),
                fstd("some_24"),
                fstd("some_25"),
                fstd("some_26"),
                fstd("some_27"),
                fstd("some_28"),
                fstd("some_29"),
                fstd("some_30"),
                fstd("some_31"),
                fstd("some_32"),
                fptr("lastVisit"),
                fstr("shortName"),
                fstd("some_33"),
                fstd("some_34"),
				fstd("some_345"),
                flng("some_35"),
                fstd("some_36"),
                flng("some_37"),
                flng("some_38"),
                flng("some_39"),
                fstd("some_40"),
                fstd("some_41"),
                fstd("some_42"),
                fstd("some_43"),
                fstd("some_44"),
                fstd("some_45"),
                fstd("some_46"),
                fstr("some_47"),
                fstd("some_48"),
                fptr("joinDate"),
                flng("some_49"),
                fint("some_50"),
                fptr("integer_2"),
                fstd("some_51"),
                fint("some_52"),
                fint("some_53"),
                flng("some_54"),
                fstd("some_55"),
                fstd("some_56"),
                fptr("some_57"),
                fstr("answer"),
                fstr("questing"),
                fstd("some_58"),
                fint("some_59"),
                fint("some_60"),
                fstr("some_61"),
				fstd("some_615"),
                fdbl("some_62"),
                fptr("spreadCalculator"),
                fstd("some_63"),
                fstd("some_64"),
                flng("some_65"),
            });

            add("com.optionfair.client.model.TrackingInfo/2864454600", {
                fstr("serialuid"),
                fstd("some_1"),
                flng("some_2"),
                flng("some_3"),
                fstd("some_4"),
                fstd("some_5"),
                fstd("some_6"),
                fstd("some_7"),
                fstd("some_8"),
                fstd("some_9"),
                flng("some_10"),
                fptr("date")
            });


            add("com.optionfair.client.verification.FieldVerification/2630919542", {
                fstd("some_1"),
                fstd("some_2"),
                fstd("some_3"),
                fstd("some_4"),
                fstd("some_5"),
                fint("some_6"),
                fstd("some_7"),
                fstr("name"),
                fstd("some_8"),
                fstr("regex"),
                fstd("some_9"),
                fstd("some_10")
            });

            add("com.optionfair.client.verification.FieldVerification/2808143106", {
                fstd("some_1"),
                fstr("some_2"),
                fstd("some_3"),
                fstd("some_4"),
                fstd("some_5"),
                fint("some_6"),
                fstd("some_7"),
                fstd("some_8"),
                fstr("name"),
                fstd("some_9"),
                fstr("regex"),
                fstd("some_10"),
                fstd("some_11"),
                fstd("some_12")
            });

            add("com.optionfair.client.model.State/1216850312", {
                fstr("serialuid"),
                fstd("country"),
                fint("order"),
                fstr("some_1"),
                fstr("name")
            });

            add("com.optionfair.client.model.PspCcTypeRel/300072509", {
                fstr("serialuid"),
                fint("creditCardTypeId"),
                flng("id"),
                fint("paymentMethodTypeId")
            });

            add("com.optionfair.client.model.MTServer/2069697836", {
                fstd("some_1"),
                fstd("some_2"),
                fstr("name")
            });

            add("com.optionfair.client.model.SubMarket/2699240068", {
                fstr("serialuid"),
                fstd("some_1"),
                fint("id"),
                fint("some_2"),
                fint("some_3"),
                fstr("name"),
                fint("some_4")
            });

            add("com.optionfair.client.model.Market/3553066736", {
                fstr("serialuid"),
                fint("id"),
                fint("some_2"),
                fstr("name"),
                fptr("subMarkets"),
                fstd("some_3"),
                fptr("parameters")
            });

            add("com.optionfair.client.model.InstrumentTypeFilter/1437805786", {
                fstr("serialuid"),
                fstd("some_1"),
                fstr("interval"),
                fint("order"),
                fstd("some_2"),
                fstr("title")
            });

            add("com.optionfair.client.model.InstrumentTypeInstrumentTypeSuperRel/1928340203", {
                fstr("serialuid"),
                fint("some_1"),
                fint("id"),
                fint("instrumentTypeId"),
                fint("instrumentTypeSuperId"),
                fint("order")
            });

            add("com.optionfair.client.model.CurrencyCountryLimit/3396769176", {
                fstr("serialuid"), // 6 [rO0ABXcEAAAAAA==]
                fstd("externalId"), // 37 [Admin]
                fstd("some_2"), // 1043 [XX]
                fint("currencyId"), // 1 [com.optionfair.client.model.LutResolverInitPackage/2546929785]
                fint("order"), // 1 [com.optionfair.client.model.LutResolverInitPackage/2546929785]
                fstd("some_3"), // 0
                flng("some_4"), // MNQ
                flng("some_5"), // B6Eg
                flng("some_6"), // MNQ
                flng("some_7"), // Ek$A
                flng("some_8"), // MNQ
                flng("some_9"), // MNQ
                flng("some_10"), // TEtA
                flng("some_11"), // MNQ
                flng("some_12"), // MNQ
                flng("some_13"), // D0JA
                flng("some_14"), // Yag
                flng("some_15"), // D0JA
                flng("some_16"), // HoSA+
                flng("some_17"), // Po
                flng("some_18"), // CcQ
                flng("some_19"), // Po
                flng("some_20"), // BOI
                flng("some_21"), // E4g
                flng("some_22"), // Po
                flng("some_23"), // Po
                flng("some_24"), // Po
				flng("some_25"), // Po
				flng("some_26"), // D0JA
                flng("some_27"), // D0JA
                flng("some_28"), // CcQ
			});

            add("com.optionfair.client.model.Currency/2484116626", {
                fstr("serialuid"),
                fint("id"),
                fstr("name"),
                fstr("symbol"),
                fstd("some_1"),
                fint("countryId"),
                fdbl("some_3"),
                fstd("some_4"),
                fstr("some_5"),
                fstr("abbr"),
                fstd("some_6"),
                flng("some_7")
            });

            add("com.optionfair.client.model.BreandPspRel/2784143874", {
                fstr("serialuid"),
                fint("externalId"),
                fstd("some_1"),
                fstd("some_2"),
                fstd("some_3"),
                fstd("some_4"),
                fstd("some_5"),
                fstd("some_6"),
                fstd("some_7"),
                fstd("some_8"),
                fstd("some_9"),
                fstd("some_10"),
                fstd("some_11")
            });

            add("com.optionfair.client.model.BrandRegulation/1013124097", {
                fstr("serialuid"), // 6 [rO0ABXcEAAAAAA==]
                fint("externalId"), // 37 [Admin]
                fstd("some_1"), // 1 [com.optionfair.client.model.LutResolverInitPackage/2546929785]
                fstd("some_2"), // 1680 [Swaziland]
                fstd("some_3"), // 1502 [PSE]
                fstd("some_4"), // 2800
                fstd("some_5"), // 0
                fstd("some_6"), // 3 [java.lang.Integer/3438268394]
                fstd("some_7"), // 3000
                fstd("some_8"), // 2000
                fstd("some_9"), // 1675 [SJM]
                fptr("some_10"), // -4104
                fstd("some_11"), // 141 [RU]
            });



			add("com.optionfair.client.model.DateTimeFormatOF/3059496483", {
				fstr("serialuid"),
				fstr("full_date"),
				fstr("mid_date"),
				fstr("short_date"),
				fstr("very_short"),
				fstd("some_1"),
				fstr("full_time"),
				fstr("mid_time"),
				fstr("short_time")
			});

            add("com.optionfair.client.model.Country/2926175992", {
                fstr("serialuid"),
                fptr("dateTimeFormat"),
                fstd("some_2"),
                fstd("language"),
                fint("some_3"),
                fint("some_4"),
                fstr("short_2"),
                fstr("short_3"),
                fstr("fullName"),
                fint("some_5"),
                fstr("some_6"),
                fstr("some_7")
            });

            add("com.optionfair.client.model.Asset/2513088007", {
				fstr("serialuid"),
				fdbl("some_1"),
				fdbl("some_2"),
				fdbl("some_3"),
				fstd("some_4"),
				fstd("some_5"),
				fdbl("some_6"),
				fstr("some_7"),
				fstd("some_8"),
				fdbl("some_9"),
				fint("some_10"),
				fint("some_11"),
				fint("id"),
				fdbl("some_13"),
				fint("some_14"),
				fint("some_15"),
				fptr("some_16"),
				fint("some_17"),
				fdbl("some_18"),
				fstd("some_19"),
				fstd("some_20"),
				fstd("some_21"),
				fint("some_22"),
				fdbl("some_225"),
				fint("some_23"),
				fint("some_24"),
				fdbl("some_25"),
				fint("some_26"),
				fdbl("some_27"),
				fstd("some_28"),
				fint("some_29"),
				fint("some_30"),
				fstd("some_31"),
				fint("subMarketId"),
				fstd("some_33"),
				fstr("name"),
				fstr("some_35"),
				fstr("some_36"),
				fstd("some_37"),
				fstr("some_375"),
				fdbl("some_38")
			});

            add("com.optionfair.client.model.LutType/4271068468", {
                fstr("serialuid"),
                fint("id"),
                fstr("name"),
                fptr("luts"),
                fstd("some_1"),
                fptr("parameters")
            }, "name");

            add("com.optionfair.client.model.Lut/3629739227", {
                fstr("serialuid"),
                fint("id"),
                fstr("name"),
                fint("lutTypeId")
            }, "name");

			add("com.optionfair.client.model.PaymentMethod/4059381229", {
				fstr("serialuid"),//6 [rO0ABXcEAAAAAA==]
				fint("commissionType"),//128 [UpBid]
				fdbl("some_2"),//3.000000e+02
				fint("id"),//93 [Tax]
				fint("some_4"),//1 [com.optionfair.client.model.LutResolverInitPackage/3212179674]
				fstd("some_5"),//0
				fint("some_6"),//8 [com.optionfair.client.model.Lut/3629739227]
				fstr("some_7"),//204 [Wire]
				fstd("some_8"),//0
				fstd("some_9"),//0
				fint("some_10"),//1 [com.optionfair.client.model.LutResolverInitPackage/3212179674]
				fint("regulationId"),//1502
				fstd("some_12"),//0
				fstd("some_13"),//267 [IdPay_cysec]
			});
			
			add("com.optionfair.client.model.TradeActionResponse/2831172712", {
                flng("accountId"), //ISw
                fint("some_2"),//0
                flng("money"),//Po
                fint("some_4"),//4
                fdbl("some_5", true),//0
				flng("some_55"),
                fdbl("some_6"),//1.13335
                fint("some_7"),//0
				fint("some_75"),//0
				fstr("some_8"),//2 [ACR+yRO4lkJ3cjauwhbeAQ==]
                fdbl("some_9", true),//0
                fint("some_10"),//1 [com.optionfair.client.model.TradeActionResponse/3233850074]
                flng("some_11"),//Xm$s
                fint("positionType"),//49
                flng("some_13"),//A
                fstr("some_14"),//3 [0]
                fdbl("some_15"),//1.13335
                fstr("some_16"),//4 [1.13335]
                fint("some_17"),//0
                flng("purchasedAt"),//VRdvfvC
                fint("some_19"),
				fint("some_20")
            });

            add("com.optionfair.client.model.Brand/160541779", {
                fstr("serialuid"),//        "6 [rO0ABXcEAAAAAA==]\n"
                fstr("some_1"), //        "778 [/]\n"
                fdbl("some_2"), //        "1.500000e+00\n"
                fint("some_3"), //        "1 [com.optionfair.client.model.LutResolverInitPackage/2546929785]\n"
                fint("some_4"), //        "0\n"
                fint("some_5"), //        "2835\n"
                fint("some_6"), //        "0\n"
                flng("some_7"), //        "Yag\n"
                fint("some_8"), //        "320 [Fee]\n"
                fint("some_9"), //        "0\n"
                fstr("some_10"), //        "776 [basic]\n"
                fint("some_11"), //        3 [java.lang.Integer/3438268394]\n"
                flng("some_12"), //        "MNQ\n"
                fint("some_13"), //        "0\n"
                flng("some_14"), //        "Yag\n"
                fstd("some_15"), //        "595 []\n"
                fstd("some_16"), //        "37 [Admin]\n"
                fstr("domain"), //        "779 [.*vantagefx.com]\n"
                fint("some_17"), //        "1 [com.optionfair.client.model.LutResolverInitPackage/2546929785]\n"
                fint("some_18"), //       "0\n"
                fint("some_19"), //        "1 [com.optionfair.client.model.LutResolverInitPackage/2546929785]\n"
                fint("some_20"), //        "0\n"
                fint("some_21"), //        "99999\n"
                fstr("some_22"), //        "780 [Enfinium]\n"
                fstd("some_23"), //        "10 [Mrs.]\n"
                fdbl("some_24"), //        "0.000000e+00")
                fptr("operator")
            });
            add("com.optionfair.client.model.Operator/1111544044", {
                fstr("serialuid"),//        "6 [rO0ABXcEAAAAAA==]\n"
                fptr("some_1"),
                fint("some_2"),
                fstr("some_3"),
                fint("some_4"),
                fint("some_5"),
                fptr("some_6"),
                fint("some_7"), //  1650
                fstr("some_8"), //  37 [Admin]
                fstr("some_9"), //  788 [none]
                fstr("some_10"),//  789 [http://nihul.4xp.com/digital_options/soap_do.php?wsdl]
                fint("some_11"), //                    1640
                fdbl("some_12"), //                    0.000000e+00
                fint("some_13"), //                    1666
                fdbl("some_14"), //                    1.500000e-01
                fdbl("some_15"), //                    2.000000e-01
                fint("some_16"), //                    1
                fptr("some_17"), //                    0
                fint("some_18"), //                    1
                fint("some_19"), //                    1
                fint("some_20"), //                    0
                fdbl("some_21"), //                    0.000000e+00
                fstd("some_22"), //                    77 [TransactionType]
                fdbl("some_23"), //                    0.000000e+00
                fstr("some_24"), //                    790 [0]
                fint("some_25"), //                    1 [com.optionfair.client.model.LutResolverInitPackage/2546929785]
                fint("some_26"), //                    0
                fint("some_27"), //                    0
                fint("some_28"), //                    0
                fint("some_29"), //                    0
                fint("some_30"), //                    0
            });

            add("com.optionfair.client.model.Position/1677341959", {
                fstr("serialuid"),//2 [rO0ABXcEAAAAAA==]
                flng("accountId"),//ISw
                fint("some_2"),//0
                flng("bet"),//Po
                flng("some_4"),//A
                fdbl("some_5"),//0.000000e+00
                fint("some_6"),//0
                fptr("some_7"),//0
                fdbl("some_8"),//0.000000e+00
                fdbl("some_9"),//0.000000e+00
                fptr("expiryDate"),//3 [java.util.Date/3385151746]
                flng("some_10"),//A
                flng("transactionId"),//KlvU
                flng("some_12"),//Po
                fint("some_13"),//1 [com.optionfair.client.model.Position/1677341959]
                fint("some_14"),//0
                fptr("openDate"),//3 [java.util.Date/3385151746]
                fstr("some_15"),//0
                flng("some_16"),//L2
                flng("some_17"),//L2
                flng("some_18"),//A
                flng("some_19"),//A
                fint("some_20"),//1 [com.optionfair.client.model.Position/1677341959]
                fint("positionOpenType"),//264
                fptr("option"),
                fptr("some_date_3"),
                fint("some_22"),
                fint("positionType"),//49
                flng("some_24"),//A
                flng("some_25"),//A
                fint("some_26"),//0
                fint("some_27"),//0
                flng("some_28"),//A
                flng("some_29"),//A
                fdbl("some_30"),//1.133350e+00
                fint("positionStatus"),//18
                fdbl("some_32"),//0.000000e+00
                fdbl("some_33"),//1.133350e+00
                fdbl("some_34"),//0.000000e+00
                fdbl("some_35"),//0.000000e+00
                fdbl("some_36"),//1.133350e+00
                fint("some_37"),//0
                fdbl("some_38"),//0.000000e+00
                flng("some_39"),//U7fV
                fdbl("some_40"),//0.000000e+00
            });

            add("com.optionfair.client.messages.MapMessageSource/3864621178", {
                    fptr("value")
            }, "value");

            add("net.sf.gilead.pojo.gwt.collection.ListParameter/3890945144", {
                    fptr("value")
            }, "value");

            add("net.sf.gilead.pojo.gwt.collection.MapParameter/1940144086", {
                    fptr("value")
            }, "value");

            add("net.sf.gilead.pojo.gwt.basic.IntegerParameter/3421404974", {
                    fptr("value")
            }, "value");

            add("net.sf.gilead.pojo.gwt.basic.StringParameter/2783524083", {
                    fstr("value")
            }, "value");

            add("net.sf.gilead.pojo.gwt.GwtSerializableId/2844303195", {
                    fstr("entityName"),
                    fptr("id"),
                    fstr("stringValue")
            });


            add(std::make_shared<GwtListType>("java.util.ArrayList/4159755760"));
            add(std::make_shared<GwtListType>("java.util.LinkedList/3953877921"));
            add(std::make_shared<GwtListType>("java.util.HashSet/3273092938"));
            add(std::make_shared<GwtMapType>("java.util.HashMap/1797211028", 0));
            add(std::make_shared<GwtMapType>("java.util.LinkedHashMap/3008245022", 1));
			add("java.lang.Integer/3438268394", { fint("value") }, "value");
			add("com.optionfair.client.luts.Instrument/2210327378", { fint("value") }, "value");
			add("java.lang.Long/4227064769", { flng("value") }, "value");
			add("java.lang.String/2004016611", { fstr("value") }, "value");
			add("java.util.Date/3385151746", { fdte("value") }, "value");
            add("java.lang.Double/858496421", { fdbl("value") }, "value");
        }
예제 #19
0
int main()
{
  //load ensemble list path and data path
  FILE *an_input_file=open_file("analysis_pars","r");
  int mode;
  char ens_list_path[1024],base_MZ_path[1024],obs_name[1024],meson_name[1024];
  read_formatted_from_file_expecting(ens_list_path,an_input_file,"%s","ens_list_path");
  read_formatted_from_file_expecting(base_MZ_path,an_input_file,"%s","base_MZ_path");
  read_formatted_from_file_expecting((char*)&mode,an_input_file,"%d","mode");
  read_formatted_from_file_expecting(obs_name,an_input_file,"%s","obs_name");
  read_formatted_from_file_expecting(meson_name,an_input_file,"%s","meson_name");
  fclose(an_input_file);
  
  //load ensembles list and parameters
  char **base_corrs_path,**ens_name;
  int nens,*T,*ibeta,*iml_un,*nlights,*nmass;
  double **mass;
  load_ensembles_list(base_corrs_path,ens_name,nens,T,ibeta,nmass,mass,iml_un,nlights,ens_list_path);
  
  //load all ensembles data
  bvec *aM,*Z;
  load_all_ensembles_MZ(aM,Z,nens,T,ibeta,nlights,nmass,base_MZ_path,obs_name,ens_name,base_corrs_path,mode);
  
  //compute f
  bvec *M=new bvec[nens];
  bvec *f=new bvec[nens];
  bvec *phi=new bvec[nens];
  bvec *x=new bvec[nens];
  for(int iens=0;iens<nens;iens++)
    {
      int b=ibeta[iens];
      
      M[iens]=aM[iens]/lat[b];
      
      if(string(obs_name)==string("P5P5"))
	{
	  f[iens]=sqrt(Z[iens])/(sinh(aM[iens])*aM[iens])/lat[b];
	  phi[iens]=f[iens]*sqrt(M[iens]);
	  int lim1;
	  if(mode==0) lim1=nmass[iens];
	  else lim1=nlights[iens];
	  for(int ims=0;ims<lim1;ims++)
	    for(int imc=ims;imc<nmass[iens];imc++)
	      {
		int ic=icombo(ims,imc,nmass[iens],nlights[iens],mode);
		f[iens].data[ic]*=mass[iens][ims]+mass[iens][imc];
		phi[iens].data[ic]*=mass[iens][ims]+mass[iens][imc];
		cout<<iens<<" ens, ic="<<ic<<", il="<<mass[iens][ims]<<", ih="<<mass[iens][imc]<<", Z="<<Z[iens][ic]<<", M="<<M[iens][ic]<<endl;
	      }
	    x[iens]=phi[iens];
	}
      
      if(string(obs_name)==string("S0S0"))
	{
	  f[iens]=(1/Zp_fr_Zs[b])*sqrt(Z[iens])/(sinh(aM[iens])*aM[iens])/lat[b];
	  phi[iens]=f[iens]*sqrt(M[iens]);
	  int lim1;
	  if(mode==0) lim1=nmass[iens];
	  else lim1=nlights[iens];
	  for(int ims=0;ims<lim1;ims++)
	    for(int imc=ims;imc<nmass[iens];imc++)
	      {
		int ic=icombo(ims,imc,nmass[iens],nlights[iens],mode);
		f[iens].data[ic]*=mass[iens][imc]-mass[iens][ims];
		phi[iens].data[ic]*=mass[iens][imc]-mass[iens][ims];
		cout<<iens<<" ens, ic="<<ic<<", il="<<mass[iens][ims]<<", ih="<<mass[iens][imc]<<", Z="<<Z[iens][ic]<<", M="<<M[iens][ic]<<endl;
	      }
	  x[iens]=phi[iens];
	}
      
      if(string(obs_name)==string("VKVK"))
	{
	  x[iens]=f[iens]=Za[b]*sqrt(Z[iens])/aM[iens]/lat[b];
	  x[iens]*=sqrt(aM[iens]/lat[b]);
	}
      
      if(string(obs_name)==string("AKAK"))
	{
	  x[iens]=f[iens]=Zv[b]*sqrt(Z[iens])/aM[iens]/lat[b];
	  x[iens]*=sqrt(aM[iens]/lat[b]);
	}
    }
  
  bvec fint(nens,nboot,njack);
  for(int iens=0;iens<nens;iens++)
    {
      if(string(meson_name)==string("Ds"))
	fint.data[iens]=interpolate_charm_strange(x[iens],nmass[iens],nlights[iens],mass[iens],ibeta[iens],mode,combine("out%02d",iens).c_str());
      
      if(string(meson_name)==string("D"))
	fint.data[iens]=interpolate_charm(x[iens],nmass[iens],nlights[iens],mass[iens],ibeta[iens],mode,combine("out%02d",iens).c_str())[iml_un[iens]];
      
      if(string(meson_name)==string("K"))
	fint.data[iens]=interpolate_unitary_light_strange(x[iens],nmass[iens],nlights[iens],iml_un[iens],mass[iens],ibeta[iens],mode,combine("out%02d",iens).c_str());
      
      if(string(meson_name)==string("Pi"))
	fint.data[iens]=x[iens][icombo(iml_un[iens],iml_un[iens],nmass[iens],nlights[iens],mode)];
      
      cout<<iens<<" "<<ibeta[iens]<<" "<<mass[iens][iml_un[iens]]<<" "<<fint.data[iens]<<endl;
    }
  
  if(string(meson_name)==string("K")) fint.write_to_binfile(combine("interpolated_f_%s",meson_name).c_str());
  else
    if(string(meson_name)==string("Pi")) fint.write_to_binfile(combine("f_%s",meson_name).c_str());
    else fint.write_to_binfile(combine("interpolated_phi_%s",meson_name).c_str());
  
  return 0;
}
예제 #20
0
 fint BB::IDOfDefAfter(PDefListElem* d) {
   PDefListElem* nextd = d->next();
   const fint BIG = ~(fint(1) << fint(sizeof(fint) * BitsPerByte  -  1));
   assert(sizeof(fint) != 4  ||  BIG == 0x7fffffff, "just checking");
   return nextd == NULL  ?  BIG  : nextd->data()->node->num();
 }