read_input_file.cpp

利用c++语言写的三维FDTD

			ReadFileVar<long>(inClientFile, "n_z_W_b", &Inp_D->n_z_W_b);

			ReadFileVar<long>(inClientFile, "Save_nr_W_DIV_x", &Inp_D->Save_nr_W_DIV_x);
			ReadFileVar<long>(inClientFile, "Save_nr_W_DIV_y", &Inp_D->Save_nr_W_DIV_y);
			ReadFileVar<long>(inClientFile, "Save_nr_W_DIV_z", &Inp_D->Save_nr_W_DIV_z);
			
		#endif

		ReadFileVar<long>(inClientFile, "source_type", &Inp_D->source_type);
		ReadFileVar<long>(inClientFile, "jel_plane_wave", &Inp_D->jel_plane_wave);

		ReadFileVar<long>(inClientFile, "n_kxkykz", &Inp_D->n_kxkykz);

		if ( Inp_D->jel_plane_wave) //plane wave
		{
			ReadFileVar<long>(inClientFile, "n_TS", &Inp_D->n_TS);
			
			ReadFileVar<int>(inClientFile, "jel_TS_xa", &Inp_D->jel_TS_xa);
			ReadFileVar<int>(inClientFile, "jel_TS_xb", &Inp_D->jel_TS_xb);
			ReadFileVar<int>(inClientFile, "jel_TS_ya", &Inp_D->jel_TS_ya);
			ReadFileVar<int>(inClientFile, "jel_TS_yb", &Inp_D->jel_TS_yb);
			ReadFileVar<int>(inClientFile, "jel_TS_za", &Inp_D->jel_TS_za);
			ReadFileVar<int>(inClientFile, "jel_TS_zb", &Inp_D->jel_TS_zb);

			ReadFileVar<long>(inClientFile, "nPML_1D", &Inp_D->nPML_1D);
			ReadFileVar<double>(inClientFile, "PML_eps_r_1D", &Inp_D->PML_eps_r_1D);
			ReadFileVar<double>(inClientFile, "PML_mu_r_1D", &Inp_D->PML_mu_r_1D);

			ReadFileVar<double>(inClientFile, "teta", &Inp_D->teta);
			ReadFileVar<double>(inClientFile, "phi", &Inp_D->phi);
			ReadFileVar<double>(inClientFile, "gamma", &Inp_D->gamma);
			
			ReadFileVar<double>(inClientFile, "X0", &Inp_D->X0);
			switch ( Inp_D->source_type)
			{
				case 1: //Gaussian source
					ReadFileVar<double>(inClientFile, "t0", &Inp_D->t0);
					ReadFileVar<double>(inClientFile, "tw", &Inp_D->tw);
					break;
				case 2: //Sin source
					ReadFileVar<double>(inClientFile, "omega", &Inp_D->omega);
					ReadFileVar<double>(inClientFile, "const_alfa", &Inp_D->const_alfa);
					break;
				case 3://Gaussian-Sin source
					ReadFileVar<double>(inClientFile, "t0", &Inp_D->t0);
					ReadFileVar<double>(inClientFile, "tw", &Inp_D->tw);
					ReadFileVar<double>(inClientFile, "omega", &Inp_D->omega);
					break;
			}
		}
		else
		{
			ReadFileVar<long>(inClientFile, "n_Coord", &Inp_D->n_Coord);
			ReadFileVarString(inClientFile, "path_name_CoordPointSource", &Inp_D->path_name_CoordPointSource);
			ReadFileVar<int>(inClientFile, "readFile_Param_Pt_Source", &Inp_D->readFile_Param_Pt_Source);
			if (Inp_D->readFile_Param_Pt_Source)
			{
				ReadFileVarString(inClientFile, "path_name_ParamPointSource", &Inp_D->path_name_ParamPointSource);
			}
			else
			{
				ReadFileVar<double>(inClientFile, "X0", &Inp_D->X0);
				switch ( Inp_D->source_type)
				{
					case 1: //Gaussian source
						ReadFileVar<double>(inClientFile, "t0", &Inp_D->t0);
						ReadFileVar<double>(inClientFile, "tw", &Inp_D->tw);
						break;
					case 2: //Sin source
						ReadFileVar<double>(inClientFile, "omega", &Inp_D->omega);
						ReadFileVar<double>(inClientFile, "const_alfa", &Inp_D->const_alfa);
						break;
					case 3://Gaussian-Sin source
						ReadFileVar<double>(inClientFile, "t0", &Inp_D->t0);
						ReadFileVar<double>(inClientFile, "tw", &Inp_D->tw);
						ReadFileVar<double>(inClientFile, "omega", &Inp_D->omega);
						break;
				}
			}
			ReadFileVar<int>(inClientFile, "pt_source_Ex", &Inp_D->pt_source_Ex);
			ReadFileVar<int>(inClientFile, "pt_source_Ey", &Inp_D->pt_source_Ey);
			ReadFileVar<int>(inClientFile, "pt_source_Ez", &Inp_D->pt_source_Ez);
			ReadFileVar<int>(inClientFile, "pt_source_Hx", &Inp_D->pt_source_Hx);
			ReadFileVar<int>(inClientFile, "pt_source_Hy", &Inp_D->pt_source_Hy);
			ReadFileVar<int>(inClientFile, "pt_source_Hz", &Inp_D->pt_source_Hz);

			ReadFileVarString(inClientFile, "Path_Kx_Ky_Kz", &Inp_D->Path_Kx_Ky_Kz);
		}

		ReadFileVarString(inClientFile, "path_name_Index", &Inp_D->path_name_Index);

		ReadFileVar<long>(inClientFile, "n_Mat", &Inp_D->n_Mat);
		ReadFileVarString(inClientFile, "path_name_Mat", &Inp_D->path_name_MatParam);

		//the followed points
		ReadFileVar<int>(inClientFile, "load_foll_field_points", &Inp_D->load_foll_field_points);
		if (Inp_D->load_foll_field_points > 0)
		{
			ReadFileVar<long>(inClientFile, "n_Ind_F", &Inp_D->n_Ind_F);
			ReadFileVarString(inClientFile, "path_name_foll_field_points", &Inp_D->path_name_foll_field_points);
		}

		//the PML parameters in x,y and z direction
		ReadFileVar<double>(inClientFile, "PML_eps_r_x1", &Inp_D->PML_eps_r_x1);
		ReadFileVar<double>(inClientFile, "PML_mu_r_x1", &Inp_D->PML_mu_r_x1);
		ReadFileVar<double>(inClientFile, "PML_eps_r_x2", &Inp_D->PML_eps_r_x2);
		ReadFileVar<double>(inClientFile, "PML_mu_r_x2", &Inp_D->PML_mu_r_x2);
	
		ReadFileVar<double>(inClientFile, "PML_eps_r_y1", &Inp_D->PML_eps_r_y1);
		ReadFileVar<double>(inClientFile, "PML_mu_r_y1", &Inp_D->PML_mu_r_y1);
		ReadFileVar<double>(inClientFile, "PML_eps_r_y2", &Inp_D->PML_eps_r_y2);
		ReadFileVar<double>(inClientFile, "PML_mu_r_y2", &Inp_D->PML_mu_r_y2);
	
		ReadFileVar<double>(inClientFile, "PML_eps_r_z1", &Inp_D->PML_eps_r_z1);
		ReadFileVar<double>(inClientFile, "PML_mu_r_z1", &Inp_D->PML_mu_r_z1);
		ReadFileVar<double>(inClientFile, "PML_eps_r_z2", &Inp_D->PML_eps_r_z2);
		ReadFileVar<double>(inClientFile, "PML_mu_r_z2", &Inp_D->PML_mu_r_z2);
	
		//save the data
		ReadFileVarString(inClientFile, "path_save_data", &Inp_D->path_save_data);
	}
	catch( char * str )
    {
        cout << "Exception raised - " << str << ut_nev << '\n';
		return 1;
    }

	if (inClientFile)
	{
		delete [] inClientFile;
		inClientFile = NULL;
	}
	
	return 0;	

	
}

///////////////////////////////////////////////////////////////////////////////////////
//Reads from a previously opened file the value of a variable
///////////////////////////////////////////////////////////////////////////////////////
template <class T>int ReadFileVar(ifstream *inClientFile, char var_name[], T *var)
{
	//inClientFile - pointer to the previously opened File stream
	//var_name - contains the name of the variable
	//var - pointer to a long, the function will return the value of the variable in this
	
	int length_var_name = (int) strlen(var_name);
	char line[512];
	int i, j;
		
	while (inClientFile->getline(line,512))
	{
		if (line[0] != '/' && line[1] != '/')
		{
			i = 0;
			while (line[i] != '\0')
			{
				if (!strncmp(&line[i],var_name,length_var_name))
				{
					j = i + length_var_name;
					while (line[j] != '\0')
					{
						if ( line[j] >= '0' && line[j] <= '9')
						{
							*var = (T) atof(&line[j]);
							inClientFile->seekg( 0, ios_base::beg ); //back to the beginning of the file
							return 0;
						}
						j++;
					}
				}
				i++;
			}
		}
	}

	*var = NULL;
	cerr << var_name << " - cannot be found" << endl;
	throw "error reading input data from: ";

	return 1; //the specified variable was not found in the file
}

///////////////////////////////////////////////////////////////////////////////////////
//Reads from a previously opened file the value of a variable
///////////////////////////////////////////////////////////////////////////////////////
int ReadFileVarString(ifstream *inClientFile, char var_name[], char **var)
{
	//inClientFile - pointer to the previously opened File stream
	//var_name - contains the name of the variable
	//var - pointer to a long, the function will return the value of the variable in this
	int length_var_name = (int) strlen(var_name);
	char line[512];
	int i, j, k;
	
	while (inClientFile->getline(line,512))
	{
		if (line[0] != '/' && line[1] != '/')
		{
			i = 0;
			while (line[i] != '\0')
			{
				if (!strncmp(&line[i],var_name,length_var_name))
				{
					j = i + length_var_name;
					while (line[j] != '\0')
					{
						if ( ( (line[j] >= 48 && line[j] <= 59 ) || // "1-9"
							   (line[j] >= 65 && line[j] <= 90 ) || // "A-Z"
							   (line[j] >= 97 && line[j] <= 122) || // "a-z"
							    line[j] == 47 || line[j] == 92 ) )    // "/" and "\"
						{  
							k = j+1;
							while (line[k] != ' ' && line[k] != '\0' && line[k] != ';') //find the end
							{
								k++;
							}
							long n_array = k-j;
							*var = new char[n_array+1];
							strncpy(*var,&line[j],n_array);//does not copy null character
							var[0][n_array] = '\0';
							inClientFile->seekg( 0, ios_base::beg ); //back to the beginning of the file
							return 0;
						}
						j++;
					}
				}
				i++;
			}
		}
	}
	
	*var = NULL;
	cerr << var_name << " - cannot be found" << endl;
	throw "error reading input data from: ";
	return 1; //the specified variable was not found in the file
}