//unix program, reading Ramo_table.txt(FEM solutions) and make PotentialArray 
//2010.1.21 9.19 (T. Kondo)

#include <iostream.h>
#include <stdio.h>
#include <fstream.h>
#include <iomanip>
using namespace std;


int main()
{
ifstream fin("temp.DAT");
ofstream fout("GetPotential_new.cxx"); 
int m=81;
int n=115;
double X[m][n], Y[m][n], Potential[m][n], E[m][n],j[m][n], PDens[m][n], Ey[m][n], Ex[m][n], Jx[m][n], Jy[m][n], Epsi[m][n], Rho[m][n]; 
for (int k=0; k<m; k++){
for (int i=0; i<n; i++){
fin >> X[k][i]>>Y[k][i]>>Potential[k][i]>>E[k][i]>>j[k][i]>>PDens[k][i]>>Ex[k][i]>>Ey[k][i]>>Jx[k][i]>>Jy[k][i]>>Epsi[k][i]>>Rho[k][i]; 

//cout <<"k,i="<<k<<" "<<i<<" x,y="<<X[k][i] << " "<<Y[k][i]<<"  Potential="<<Potential[k][i]<<endl; 
}
}
fin.close();
//
//
fout << "double GetPotentialValue(int& ix, int& iy) {"<<endl;
fout << "double PotentialArray[81][115] = {"<< endl; 
fout<<"//--- PotentialArray[81][115] "<<endl;
fout<<"//      ix = 0 -  80 : 5um step, 40um - 440 um"<<endl;
fout<<"//      iy = 0 - 114 : 2.5um step, 0um - 285 um"<<endl;
fout<<"//"<<endl;
fout<< std::setprecision(6)<<std::setiosflags(std::ios::fixed)<<endl;

for (int k=0; k<m; k++){
  double x = 40. + k*5.;
  fout<<"//--- PotentialArray["<<k<<"][115] , x = "<<x<<" um --------"<<endl;
  fout << "{"<<std::endl;
  for(int j=0;j<115;j+=5)
     fout<<setw(10)<<Potential[k][j]<<",   "<<setw(10)<<Potential[k][j+1]<<",   "<<setw(10)<<Potential[k][j+2]<<",   "
     <<setw(10)<<Potential[k][j+3]<<",   "<<setw(10)<<Potential[k][j+4]<<","<<std::endl;
  fout<<"},"<<std::endl;
}
fout<<" };" << endl;
fout<<"return PotentialArray[ix][iy];" << endl;
fout<<"}"<<std::endl;
fout.close();

return 0;
}