void Look_traj()
{
//make_plot_crd_diff 110 1605 0
vector<double> X;
Load_traj(110, 1605, 0, X);
int szX = X.size();
if (szX == 0)
return;
vector<double> diffX(szX-1);
for (int i = 0; i < szX-1; ++i)
{
diffX[i] = X[i+1] - X[i];
}
{
//------------------------------------------------------------------
char filter[] = TEXT("Plot Text File (*.txt)\0*.txt\0")
TEXT("All Files (*.*)\0*.*\0");
char outfile_name[1048];
sprintf(outfile_name, "\0");
DWORD nFilterIndex;
if (SaveFileDlg(0, outfile_name, filter, nFilterIndex) == S_OK)
{
cout << "now saving file " << outfile_name << endl;
ofstream ofile;
ofile.open(outfile_name, ios::out);
for (int i = 0; i < szX-1; ++i)
{
ofile << X[i] << "," << diffX[i] << endl;
}
ofile.close();
}
}
for (int N = 10; N < 100; N += 5)
Statistika (diffX, "diffX", N, calcMathOzh);
for (int N = 10; N < 100; N += 5)
Statistika (X, "linearTrendParameterA", N, calcLinearTrendParameterA);
}
int main()
{
int pil;
gotoxy(20,0); printf("Data Mahasiswa JTK Polban Kelas 1-B 2011\n\n");
printf("1. Input data mahasiswa\n");
printf("2. Tampilkan data mahasiwa\n");
printf("3. Daftar Mahasiswa yang tinggal di Bandung\n");
printf("4. Daftar Mahasiswa laki-laki yang tingginya lebih dari 170cm\n");
printf("5. Statistika Mahasiswa berdasarkan kota\n\n");
printf("Masukkan pilihan anda = "); scanf("%d",&pil);
switch(pil)
{
case 1 : Input(); break;
case 2 : Tampil(); break;
case 3 : Bandung(); break;
case 4 : CariMhs(); break;
case 5 : Statistika(); break;
default : printf("Nomor yang anda masukkan salah!\n");
}
return 0;
}
int Look_txt()
{
TCHAR filter[] = TEXT("Ghemical MD results File (*.txt)\0*.txt\0")
TEXT("All Files (*.*)\0*.*\0");
TCHAR fpath[1024];
TCHAR filename[1024];
sprintf(filename, "\0");
{
DWORD nFilterIndex;
vector<string> names;
vector<string> *pnames = &names;
vector<vector<double> > vectors;
vectors.reserve(2000000);
while (OpenFileDlg(0, filter, fpath, nFilterIndex) == S_OK)
{
ReadDatFile(NULL, fpath, filename, &vectors, pnames);
pnames = NULL;
printf("\nfilename %s\n\n", filename);
int cols = names.size();
int rows = vectors.size();
#if WRITE_LOCKED_FORCES
int cMom = 4 - 1;
int cVx = 5 - 1;
int cFxup = 14 - 1;
int cFxdw = 17 - 1;
int cVxup = 8 - 1;
int cVxdw = 11 - 1;
#endif
#if WRITE_WORKED_FORCES
int cMom = 4 - 1;
int cVx = 5 - 1;
int cVxup = 14 - 1;
int cVxdw = 17 - 1;
int cVx_wk_up = 8 - 1;
int cVx_wk_dw = 11 - 1;
int cFx_wk_up = 20 - 1;
int cFx_wk_dw = 23 - 1;
#endif
vector<double> means(cols, 0.0);
printf("vectors.size() = %d\n",rows);
printf("names.size() = %d\n", cols);
for (vector<vector<double> >::iterator it = vectors.begin();
it != vectors.end(); it++)
{
for (int c = 0; c < cols; c++)
{
means[c] += (*it).operator [](c);
}
}
for (int c = 0; c < cols; c++)
{
means[c] /= rows;
printf("mean(%s) = %f\n", names[c].c_str(), means[c]);
}
#if WRITE_LOCKED_FORCES || WRITE_WORKED_FORCES
int r0 = 0;
cout << "enter r0\n";
cin >> r0;
#endif
#if WRITE_LOCKED_FORCES
vector<double> dF(rows-r0);
for (int r = r0; r < rows; r++)
{
dF[r-r0] = vectors[r][cFxup] - vectors[r][cFxdw];
}
Statistika (dF, "dF");
vector<double> Mom(rows-r0);
for (r = r0; r < rows; r++)
{
Mom[r-r0] = vectors[r][cMom];
}
Statistika (Mom, "Mom");
vector<double> dV(rows-r0);
for (r = r0; r < rows; r++)
{
dV[r-r0] = vectors[r][cVxup] - vectors[r][cVxdw];
}
Statistika (dV, "dV");
vector<double> Vx(rows-r0);
for (r = r0; r < rows; r++)
{
Vx[r-r0] = vectors[r][cVx];
}
Statistika (Vx, "Vx");
#endif
#if WRITE_WORKED_FORCES
vector<double> dF_wk(rows-r0);
for (int r = r0; r < rows; r++)
{
dF_wk[r-r0] = vectors[r][cFx_wk_up] - vectors[r][cFx_wk_dw];
}
Statistika (dF_wk, "dF_wk");
vector<double> dV_wk(rows-r0);
for (r = r0; r < rows; r++)
{
dV_wk[r-r0] = vectors[r][cVx_wk_up] - vectors[r][cVx_wk_dw];
}
Statistika (dV_wk, "dV_wk");
//if (!worked[n1])
vector<double> Mom(rows-r0);
for (r = r0; r < rows; r++)
{
Mom[r-r0] = vectors[r][cMom];
}
Statistika (Mom, "Mom");
vector<double> dV(rows-r0);
for (r = r0; r < rows; r++)
{
dV[r-r0] = vectors[r][cVxup] - vectors[r][cVxdw];
}
Statistika (dV, "dV");
vector<double> Vx(rows-r0);
for (r = r0; r < rows; r++)
{
Vx[r-r0] = vectors[r][cVx];
}
Statistika (Vx, "Vx");
#endif
}
}
/*else
{
DWORD nFilterIndex;
if (SaveFileDlg(0, filename, filter, nFilterIndex) == S_OK)
{
SetDlgItemText(ref->hDlg,IDC_EDIT_TRAJFILE2, filename);
}
}*/
printf("Hello World!\n");
return 0;
}
