//unix program, reading SCTstrips_Data.txt(FEM solutions) and make ExArray, EyArray
//2010.1.18, 4.13, 9.16 (T. Kondo)

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


int main()
{
ifstream fin("temp.DAT");
int n=115;
int m=17;
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]; 

//----------- in case of under depletion ---------------
if (i == 0 && Y[k][i] != 0.) {
    int j = int ( Y[k][i] /2.5 );
     i = j;
     E[k][i] = E[k][0]/100.;
     Ex[k][i] = Ex[k][0]/100.;
     Ey[k][i] = Ey[k][0]/100.;
     for (int jj=0 ; jj<j ; jj++ ) {
       Ex[k][jj] = 0.;
       Ey[k][jj] = 0.;
     }
} else { 
  E[k][i] = E[k][i]/100.;
  Ex[k][i] = Ex[k][i]/100.;
  Ey[k][i] = Ey[k][i]/100.;
}
// cout <<"k,i="<<k<<" "<<i<<" x,y="<<X[k][i] << " "<<Y[k][i]<<" "<<E[k][i]<<" Ex="<<Ex[k][i]<< " Ey=" <<Ey[k][i]<<endl; 
//-----------special cases for under full depletion -------------------------
  if (Ey[k][i] > 0.) continue;
  Ex[k][i] = 0.;
  Ey[k][i] = 0.;
}
}
fin.close();
//---------------------------------------------------------------------

cout<<"#include <math.h>"<<endl;
cout<< std::setprecision(2)<<std::setiosflags(std::ios::fixed)<<std::endl;
cout<<"double EyValueXXX(int& ix, int& iy) {"<<endl;
cout<<"/"<<"/"<<"  EyArray for VD=70V(rho=0.1540[C/m**3]), VB=XXXV" << endl;
cout << "double EyArray[17][115] = {"<< endl; 
for (int k=0; k<m; k++){
  cout << "{"<<endl;
   for (int n=0; n<115; n+=5) {
     cout<<std::setw(10)<<Ey[k][n]<<","<<std::setw(10)<<Ey[k][n+1]<<","<<std::setw(10)<<Ey[k][n+2]<<","
     <<std::setw(10)<<Ey[k][n+3]<<","<<std::setw(10)<<Ey[k][n+4]<<","<<endl;
   }
   cout<<"}, "<<endl;
}
cout<<" };" << endl;
cout<<"/"<<"/"<<endl;
cout<<"return EyArray[ix][iy];"<<endl;
cout<<" } " << endl;
cout<<"                                      "<<endl;
cout<<"                                      "<<endl;
cout<<"                                      "<<endl;
cout<<"double ExValueXXX(int& ix, int& iy) {"<<endl;
cout<<"/"<<"/"<<"  ExArray for VD=70V(rho=0.1540[C/m**3]), VB=XXXV" << endl;
cout << "double ExArray[17][115] = {"<< endl; 
for (int k=0; k<m; k++){

   cout << "{"<<endl;
   for (int n=0; n<115; n+=5) {
     cout<<std::setw(10)<<Ex[k][n]<<","<<std::setw(10)<<Ex[k][n+1]<<","<<std::setw(10)<<Ex[k][n+2]<<","
     <<std::setw(10)<<Ex[k][n+3]<<","<<std::setw(10)<<Ex[k][n+4]<<","<<endl;
   }
   cout<<"}, "<<endl;
}
cout<<" };" << endl;
cout<<"/"<<"/"<<endl;
cout<<"return ExArray[ix][iy];"<<endl;
cout<<" } " << endl;
cout<<"                                      "<<endl;
cout<<"                                      "<<endl;
cout<<"                                      "<<endl;
return 0;
}