Exemplo n.º 1
0
int main(int argc, char **argv) {
    if (argc != 2) {
        std::cout << "Usage: HitoriSolver <table.txt>\n";
        return 0;
    }

    std::ifstream in(argv[1]);
    if (!in.good()) {
        std::cout << "Unable to open " << argv[1] << " for reading\n";
        return 1;
    }

    try {
        Table t = readTable(in);
        std::cout << solve(t);

    } catch (std::exception e) {
        std::cout << e.what() << std::endl;
        return 2;
    }

    in.close();

    return 0;
}
Exemplo n.º 2
0
QueryEngine::QueryEngine(const char * filepath){
	
	//Create filepaths
	std::string basepath(filepath);
	std::string path_to_hashtable = basepath + "/probing_hash.dat";
	std::string path_to_data_bin = basepath + "/binfile.dat";
	std::string path_to_vocabid = basepath + "/vocabid.dat";

	//Read config file
	std::string line;
	std::ifstream config ((basepath + "/config").c_str());
	getline(config, line);
	int tablesize = atoi(line.c_str()); //Get tablesize.

	getline(config, line);
	largest_entry = atoi(line.c_str()); //Set largest_entry.
	config.close();

	//Mmap binary table
	struct stat filestatus;
	stat(path_to_data_bin.c_str(), &filestatus);
	binary_filesize = filestatus.st_size;
	binary_mmaped = read_binary_file(path_to_data_bin.c_str(), binary_filesize);

	//Read hashtable
	size_t table_filesize = Table::Size(tablesize, 1.2);
	mem = readTable(path_to_hashtable.c_str(), table_filesize);
	Table table_init(mem, table_filesize);
	table = table_init;

	//Read vocabid
	read_map(&vocabids, path_to_vocabid.c_str());

	std::cout << "Initialized successfully! " << std::endl;
}
Exemplo n.º 3
0
QueryEngine::QueryEngine(const char * filepath) : decoder(filepath){
    
    //Create filepaths
    std::string basepath(filepath);
    std::string path_to_hashtable = basepath + "/probing_hash.dat";
    std::string path_to_data_bin = basepath + "/binfile.dat";
    std::string path_to_source_vocabid = basepath + "/source_vocabids";

    ///Source phrase vocabids
    read_map(&source_vocabids, path_to_source_vocabid.c_str());

    //Target phrase vocabIDs
    vocabids = decoder.get_target_lookup_map();

    //Read config file
    std::string line;
    std::ifstream config ((basepath + "/config").c_str());
    getline(config, line);
    int tablesize = atoi(line.c_str()); //Get tablesize.
    config.close();

    //Mmap binary table
    struct stat filestatus;
    stat(path_to_data_bin.c_str(), &filestatus);
    binary_filesize = filestatus.st_size;
    binary_mmaped = read_binary_file(path_to_data_bin.c_str(), binary_filesize);

    //Read hashtable
    size_t table_filesize = Table::Size(tablesize, 1.2);
    mem = readTable(path_to_hashtable.c_str(), table_filesize);
    Table table_init(mem, table_filesize);
    table = table_init;

    std::cerr << "Initialized successfully! " << std::endl;
}
Exemplo n.º 4
0
void NetworkClient::init(char* address)
{
    open(address);
    writeHeader(1, 0x0);
    writeBuffer();
    readTable();
}
void TestApplication::SetUp() {
    table_.clear();

    std::stringstream ss;
    ss << kTableSource;
    readTable(&ss);
}
Exemplo n.º 6
0
int SP_NKEndPointTableConfig :: init( SP_NKIniFile * iniFile )
{
	if( NULL != mTable ) delete mTable, mTable = NULL;

	mTable = readTable( iniFile );

	return NULL != mTable ? 0 : -1;
}
Exemplo n.º 7
0
Foam::interpolationTable<Type>::interpolationTable(const dictionary& dict)
:
    List<Tuple2<scalar, Type> >(),
    boundsHandling_(wordToBoundsHandling(dict.lookup("outOfBounds"))),
    fileName_(dict.lookup("fileName")),
    reader_(tableReader<Type>::New(dict))
{
    readTable();
}
Exemplo n.º 8
0
Foam::interpolationTable<Type>::interpolationTable(const fileName& fName)
:
    List<Tuple2<scalar, Type> >(),
    boundsHandling_(interpolationTable::WARN),
    fileName_(fName),
    reader_(new openFoamTableReader<Type>())
{
    readTable();
}
Exemplo n.º 9
0
SP_NKEndPointTable * SP_NKEndPointTableConfig :: readTable( const char * configFile )
{
	SP_NKEndPointTable * table = NULL;

	SP_NKIniFile iniFile;

	if( 0 == iniFile.open( configFile ) ) table = readTable( &iniFile );

	return table;
}
Exemplo n.º 10
0
void medXMLToLUTReaderPrivate::readLUT()
{
    Q_ASSERT(xml.isStartElement() && xml.name() == "medLUTs");

    while (xml.readNextStartElement()) {
        if (xml.name() == "table")
            readTable();
        else
            xml.skipCurrentElement();
    }

}
Exemplo n.º 11
0
void DummyBird::initParameters()
{ 
  if (!nh_private_.getParam("gravity_constant", g))
    g = 9.81;
  if (!nh_private_.getParam("heli_mass",heli_mass))
    heli_mass = 615;
  if (!nh_private_.getParam("ROS_thrust_cmd_table",strThrustCmd))
    strThrustCmd = "0 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5";
  if (!nh_private_.getParam("actual_thrust_table",strActThrust))
    strActThrust = "57.8 95.3 138.3 192.3 265.3 350.3 430.3 540.3 665.3 790.3 932.3";
  if (!nh_private_.getParam("quadrotor_type",quadrotorType))
    quadrotorType = "Hummingbird";

  thrustCmd = readTable(strThrustCmd);
  actualThrust = readTable(strActThrust);
  ROS_INFO("Gravity constant initialized %.3f",g);
  ROS_INFO("Helicopter mass initialized %.3f",heli_mass);
  ROS_INFO("Quadrotor Type %s",quadrotorType.c_str());
  ROS_INFO("Thrust table: ");
  for (int i = 0; (unsigned)i<thrustCmd.size(); i++)
    ROS_INFO("%.3f : %.3f",thrustCmd[i],actualThrust[i]);
}
Exemplo n.º 12
0
int TLKFile::init(int n)
{
	if(n == -1) n = initstages;
	if((n > 0) && !inited[0] && readHeader()) return errcode;
	if((n > 1) && !inited[1] && readTable()) return errcode;
	if((n > 2) && !inited[2] && readStrings()) return errcode;
	if((n > 3) && !inited[3])
	{
		close();
		inited.set(3);
	}
	return errcode = 0;
}
Exemplo n.º 13
0
int Decoder(List* head, char* filename, char* outname)
{
	FILE* fIn  = NULL;
	FILE* fOut = NULL;

	Node* node = NULL;

	char chname[128]	  = {0, };
	char inname[128]	  = {0, };
	char line_buffer[256] = {0, };
	
	char* line_ptr		  = NULL;
	
	int len = 0;
	int i;

	readTable(head, filename);

	for(i=0; filename[i]!='.'; i++)	chname[i] = filename[i];
	sprintf(inname, "%s_encoded.txt", chname);
	if((fIn = fopen(inname, "rt")) == NULL)	return -1;

	sprintf(outname, "%s_decoded.txt", chname);
	if((fOut = fopen(outname, "wt")) == NULL)	return -1;

	while(feof(fIn) == 0)
	{
		fgets(line_buffer, 256, fIn);
		line_ptr = line_buffer;

		node = head->pHead;
		while(node)
		{
			len = strlen(node->Code);
			if(strncmp(line_ptr, node->Code, len) == 0)
			{
				fprintf(fOut, "%c", node->ascii);
				line_ptr += len;
				node = head->pHead;
			}
			else	node = node->pNext;
		}

	}

	Destroy(head);

	fclose(fIn);
	fclose(fOut);
	return 1;
}
Exemplo n.º 14
0
/** Reads a file containing a list of Julian days and measurements
 * 
 * @param[in] fileName the name of a file to be read. The file 
 *	is assumed to be formatted as a 2xN comma-delimited table, 
 *	with the first column a floating point Julian date and the second 
 *	a floating point measurement. The file may also contain comment lines 
 *	preceded by '#'.
 * @param[out] timeVec a vector containing the times of each 
 *	observation
 * @param[out] dataVec a vector containing the measurement (typically flux 
 *	or magnitude) observed at each time
 *
 * @post @p timeVec is sorted in ascending order
 * @post @p timeVec.size() = @p dataVec.size()
 * @post for all i, @p dataVec[i] is the measurement taken at @p timeVec[i]
 *
 * @exception std::bad_alloc Thrown if there is not enough memory to store 
 *	the data.
 * @exception kpfutils::except::FileIo Thrown if any file operation fails.
 *
 * @exceptsafe The function arguments are unchanged in the event of an exception.
 */
void readCsvLightCurve(const string& fileName, DoubleVec &timeVec, DoubleVec &dataVec) {
	// copy-and-swap
	vector<double> tempTimes, tempData;
	
	boost::shared_ptr<FILE> hInput = fileCheckOpen(fileName, "r");

	readTable(hInput.get(), " %lf , %lf", tempTimes, tempData);
	sortByTime(tempTimes, tempData);

	// IMPORTANT: no exceptions beyond this point
	
	swap(timeVec, tempTimes);
	swap(dataVec, tempData );
}
Exemplo n.º 15
0
/** Reads a file containing a list of obsids, Julian days, measurements, errors, and limits
 * 
 * @param[in] fileName the name of a file to be read. The file 
 *	is assumed to be formatted as a 5xN space-delimited table, with 
 *	the first column a running index, the second a floating point 
 *	Julian date, the third a floating point measurement, the fourth 
 *	a floating point error, and the fifth a detection limit in the same 
 *	units as the measurement. The file may also contain comment lines 
 *	preceded by '#'.
 * @param[in] errMax the maximum error to tolerate in a data point. Any 
 *	points with an error exceeding @p errMax are ignored.
 * @param[out] timeVec a vector containing the times of each 
 *	observation
 * @param[out] dataVec a vector containing the measurement (typically flux 
 *	or magnitude) observed at each time
 * @param[out] errVec a vector containing the error on each measurement
 *
 * @post @p timeVec is sorted in ascending order
 * @post @p timeVec.size() = @p dataVec.size() = @p errVec.size()
 * @post for all i, @p dataVec[i] &plusmn; @p errVec[i] is the 
 *	measurement taken at @p timeVec[i]
 * @post for all i, @p errVec[i] &le; @p errMax
 *
 * @exception std::bad_alloc Thrown if there is not enough memory to store 
 *	the data.
 * @exception kpfutils::except::FileIo Thrown if any file operation fails.
 *
 * @exceptsafe The function arguments are unchanged in the event of an exception.
 *
 * @bug Current implementation ignores limits.
 */
void readWg2LightCurve(const string& fileName, double errMax, DoubleVec &timeVec, 
		DoubleVec &dataVec, DoubleVec &errVec) {
	// copy-and-swap
	vector<double> tempTimes, tempData, tempErrs;
	
	boost::shared_ptr<FILE> hInput = fileCheckOpen(fileName, "r");

	readTable(hInput.get(), " %*i %lf %lf %lf %*lf", tempTimes, tempData, tempErrs);
	errorFilter(errMax, tempTimes, tempData, tempErrs);
	sortByTime(tempTimes, tempData, tempErrs);

	// IMPORTANT: no exceptions beyond this point
	
	swap(timeVec, tempTimes);
	swap(dataVec, tempData );
	swap( errVec, tempErrs );
}
Exemplo n.º 16
0
int main(void) {
	DIR *dir = opendir("/sys/acpi");
	if(!dir)
		error("Unable to open /sys/acpi");

	size_t i = 0;
	struct dirent e;
	while(readdirto(dir,&e)) {
		sRSDT table;
		if(strcmp(e.d_name,".") == 0 || strcmp(e.d_name,"..") == 0)
			continue;
		readTable(e.d_name,&table);
		printf("ACPI Table %zu:\n",i);
		printf("\tSignature: %.4s\n",(char*)&table.signature);
		printf("\tLength: %u\n",table.length);
		printf("\tRevision: %u\n",table.revision);
		printf("\tOEMID: %.6s\n",table.oemId);
		printf("\tOEMTableID: %.8s\n",table.oemTableId);
		i++;
	}

	closedir(dir);
	return EXIT_SUCCESS;
}
Exemplo n.º 17
0
void initwwwPot(struct parameters *par,struct systemPos pos){
  int error,i,j,k;
  char filename[20];
  struct bondList *blist;
  getBondList(&blist);
  
  strcpy(filename,"wwwPot.par");
  openParameterFile(filename);
  error=0;
  
  error+=getParValue("includeKeatingPot",&wwwPar.includeKeatingPot,"%d");
  error+=getParValue("includeRepPot",&wwwPar.includeRepPot,"%d");
  error+=getParValue("includeSnnRep",&wwwPar.includeSnnRep,"%d");
  error+=getParValue("includeSuboxidePot",&wwwPar.includeSuboxidePot,"%d");
  error+=getParValue("includeDanglbondPot",&wwwPar.includeDanglbondPot,"%d");
  error+=getParValue("includeRingPot",&wwwPar.includeRingPot,"%d");
  
  if(wwwPar.includeKeatingPot==1){/*original keatin */
    error+=getParValue("okp1_kStretch_SiSi",&wwwPar.kStretch[1][1],"%lf");
    error+=getParValue("okp1_kStretch_SiO",&wwwPar.kStretch[1][0],"%lf");
    error+=getParValue("okp1_kBend_SiSiSi",&wwwPar.kBend[1][1][1],"%lf");
    error+=getParValue("okp1_kBend_SiOSi",&wwwPar.kBend[1][0][1],"%lf");
    error+=getParValue("okp1_kBend_OSiO",&wwwPar.kBend[0][1][0],"%lf");
    error+=getParValue("okp1_r0_SiSi",&wwwPar.r0[1][1],"%lf");
    error+=getParValue("okp1_r0_SiO",&wwwPar.r0[1][0],"%lf");
    error+=getParValue("okp1_cos0_Si",&wwwPar.cos0[1],"%lf");
    error+=getParValue("okp1_cos0_O",&wwwPar.cos0[0],"%lf");
  }
  else { /*simplified keatin */
    error+=getParValue("skp2_kStretch_SiSi",&wwwPar.kStretch[1][1],"%lf");
    error+=getParValue("skp2_kStretch_SiO",&wwwPar.kStretch[1][0],"%lf");
    error+=getParValue("skp2_kBend_SiSiSi",&wwwPar.kBend[1][1][1],"%lf");
    error+=getParValue("skp2_kBend_SiOSi",&wwwPar.kBend[1][0][1],"%lf");
    error+=getParValue("skp2_kBend_OSiO",&wwwPar.kBend[0][1][0],"%lf");
    error+=getParValue("skp2_r0_SiSi",&wwwPar.r0[1][1],"%lf");
    error+=getParValue("skp2_r0_SiO",&wwwPar.r0[1][0],"%lf");
    error+=getParValue("skp2_cos0_Si",&wwwPar.cos0[1],"%lf");
    error+=getParValue("skp2_cos0_O",&wwwPar.cos0[0],"%lf");
	 
  }
  if(wwwPar.includeKeatingPot==3){/*stixrudes potential, no keating at all...*/
	 error+=getParValue("stix_D",&wwwPar.stix_D,"%lf");
	 error+=getParValue("stix_beta",&wwwPar.stix_beta,"%lf");
	 error+=getParValue("stix_r0",&wwwPar.stix_r0,"%lf");
	 error+=getParValue("stix_galpha",&wwwPar.stix_galpha,"%lf");
	 error+=getParValue("stix_alpha0",&wwwPar.stix_alpha0,"%lf");
	 error+=getParValue("stix_gL",&wwwPar.stix_gL,"%lf");
	 error+=getParValue("stix_L0",&wwwPar.stix_L0,"%lf");
	 error+=getParValue("stix_A",&wwwPar.stix_A,"%lf");
	 error+=getParValue("stix_b",&wwwPar.stix_b,"%lf");
	 error+=getParValue("stix_rc",&wwwPar.stix_rc,"%lf");
	 error+=getParValue("stix_gamma",&wwwPar.stix_gamma,"%lf");
 
  }

  
  
  error+=getParValue("repulsive_r0",&wwwPar.repR0,"%lf");
  error+=getParValue("repulsive_k",&wwwPar.repK,"%lf");
  
  error+=getParValue("ringPenalty3",&wwwPar.ringPenalty3,"%lf");
  error+=getParValue("ringPenalty2",&wwwPar.ringPenalty2,"%lf");

  error+=getParValue("SubOxidePenalty+0",&wwwPar.subOxidePenalty[0],"%lf");
  error+=getParValue("SubOxidePenalty+1",&wwwPar.subOxidePenalty[1],"%lf");
  error+=getParValue("SubOxidePenalty+2",&wwwPar.subOxidePenalty[2],"%lf");
  error+=getParValue("SubOxidePenalty+3",&wwwPar.subOxidePenalty[3],"%lf");
  error+=getParValue("SubOxidePenalty+4",&wwwPar.subOxidePenalty[4],"%lf");
  closeParameterFile();
  
  for(i=0;i<2;i++) 
    for(j=0;j<2;j++){
      wwwPar.kStretch[i][j]*=0.5; 
      for(k=0;k<2;k++)
		  wwwPar.kBend[i][j][k]*=0.5;
    }  
  wwwPar.repK*=0.5;
  
  wwwPar.kBend[1][1][0]=sqrt(wwwPar.kBend[1][1][1]*wwwPar.kBend[0][1][0]);
  wwwPar.kBend[0][1][1]=wwwPar.kBend[1][1][0];
  
  
  /*make the matrices symmetric*/
  wwwPar.kStretch[0][1]=wwwPar.kStretch[1][0];
  wwwPar.r0[0][1]=wwwPar.r0[1][0];

  for(i=0;i<2;i++)
    for(j=0;j<2;j++)
      wwwPar.r02[i][j]=POW2(wwwPar.r0[i][j]);
  
  /*change the units*/
  for(i=0;i<2;i++)
    for(j=0;j<2;j++){
      wwwPar.kStretch[i][j]*=EV;
      for(k=0;k<2;k++)
		  wwwPar.kBend[i][j][k]*=EV;
    }
  for(i=0;i<5;i++)
    wwwPar.subOxidePenalty[i]*=EV;

  wwwPar.ringPenalty3*=EV;  
  wwwPar.ringPenalty2*=EV;  
  wwwPar.repK*=EV;
  wwwPar.stix_D*=EV;
  wwwPar.stix_A*=EV;
  wwwPar.stix_galpha*=EV;
  wwwPar.stix_gL*=EV;
  if(wwwPar.includeKeatingPot==3)
	 readTable("stixrep.dat",&wwwPar.stix_repTable);

  wwwPar.rcut=wwwPar.repR0;
  wwwPar.rcut2=POW2(wwwPar.rcut);
  wwwPar.repR02=POW2(wwwPar.repR0);
  
  if(wwwPar.includeSnnRep){
	 readTable("ooPottable.dat",&wwwPar.ooPottable);
	 readTable("ssPottable.dat",&wwwPar.ssPottable);
	 readTable("soPottable.dat",&wwwPar.soPottable);
  }
  if(error!=0){
    if(getRank()==0) printf("could not read the required parameters from file %s\n",filename);
    exit(0);
  }

  if(wwwPar.includeRepPot && wwwPar.includeKeatingPot==3){
	 printf("stix and reppot cannot be on at the same time");
	 exit(0);
  }

  /*init ngbrs */
  initNgbrs(&nd,pos,par,wwwPar.rcut);
  initNgbrs(&ndLastStep,pos,par,wwwPar.rcut);
  
  updateNgbrs(&nd,pos,par);
  if(wwwPar.includeRepPot)
	 nonNnNgbrs(&nd,pos.nAtoms); 
  
  
  prevInteractionHead=malloc(sizeof(int)*pos.nAtoms);
  prevInteractionList=malloc(sizeof(int)*pos.nAtoms*(blist->nMax+nd.nMax)); /* if these lists are later on larger there might  occur a sigsegv */
  reinitWwwAfterBondChange(par, pos);
  reinitWwwPotAfterStep(par,pos,1); /* this updates some stuff.. (prevInteractionHead,prevInteractionList,ndLastStep) */
}
TEST_F(TestApplication, can_read_table_from_nullptr_source) {
    // Arrange, Act, Assert
    ASSERT_NO_THROW(readTable(nullptr));
}
Exemplo n.º 19
0
void THDMcache::read(){
    

    std::string tablepath="/Users/Roma1/Desktop/HEPfit/THDM/tabs/";
    std::stringstream br1,br2,br3,br4,br5,br6,br7;
    std::stringstream pc1,pc2,pc3,pc4,pc5;
    std::stringstream dw1;
    std::stringstream cs1,cs2,cs3,cs4,cs5,cs6;
    std::stringstream ex1,ex2,ex3,ex4,ex5,ex6,ex7,ex8,ex9,ex10,ex11,ex12,ex13,ex14,ex15,ex16;

    std::cout<<"reading tables"<<std::endl;

    br1 << tablepath << "GridSM1.dat";
    array1 = readTable(br1.str(),19861);
    br2 << tablepath << "GridSM2.dat";
    array2 = readTable(br2.str(),19861);
    br3 << tablepath << "GridSM3.dat";
    array3 = readTable(br3.str(),19861); 
    br4 << tablepath << "GridSM4.dat";
    array4 = readTable(br4.str(),19861);
    br5 << tablepath << "GridSM5.dat";
    array5 = readTable(br5.str(),19861);
    br6 << tablepath << "GridSM6.dat";
    array6 = readTable(br6.str(),19861);
    br7 << tablepath << "GridSM7.dat";
    array7 = readTable(br7.str(),19861);
    pc1 << tablepath << "GridSM11.dat";
    array11 = readTable(pc1.str(),1971);
    pc2 << tablepath << "GridSM12.dat";
    array12 = readTable(pc2.str(),1971);
    pc3 << tablepath << "GridSM13.dat";
    array13 = readTable(pc3.str(),1971);
    pc4 << tablepath << "GridSM14.dat";
    array14 = readTable(pc4.str(),1971);
    pc5 << tablepath << "GridSM15.dat";
    array15 = readTable(pc5.str(),1971);
    ex1 << tablepath << "GridSM16.dat";
    array16 = readTable(ex1.str(),9851);
    ex2 << tablepath << "GridSM17.dat";
    array17 = readTable(ex2.str(),9851);
    dw1 << tablepath << "GridSM18.dat";
    array18 = readTable(dw1.str(),19861);
    cs1 << tablepath << "GridSM19.dat";
    array19 = readTable(cs1.str(),186);
    cs2 << tablepath << "GridSM20.dat";
    array20 = readTable(cs2.str(),186);
    cs3 << tablepath << "GridSM21.dat";
    array21 = readTable(cs3.str(),186);
    cs4 << tablepath << "GridSM22.dat";
    array22 = readTable(cs4.str(),186);
    cs5 << tablepath << "GridSM26.dat";
    array26 = readTable(cs5.str(),992);
    cs6 << tablepath << "GridSM27.dat";
    array27 = readTable(cs6.str(),185);
    ex3 << tablepath << "GridSM29.dat";
    array29 = readTable(ex3.str(),986);
    ex4 << tablepath << "GridSM31.dat";
    array31 = readTable(ex4.str(),985);
    ex5 << tablepath << "GridSM32.dat";
    array32 = readTable(ex5.str(),496);
    ex6 << tablepath << "GridSM33.dat";
    array33 = readTable(ex6.str(),496);
    ex7 << tablepath << "GridSM34.dat";
    array34 = readTable(ex7.str(),991);
    ex8 << tablepath << "GridSM35.dat";
    array35 = readTable(ex8.str(),496);
    ex9 << tablepath << "GridSM36.dat";
    array36 = readTable(ex9.str(),496);
    ex10 << tablepath << "GridSM37.dat";
    array37 = readTable(ex10.str(),100);
    ex11 << tablepath << "GridSM38.dat";
    array38 = readTable(ex11.str(),100);
    ex12 << tablepath << "GridSM39.dat";
    array39 = readTable(ex12.str(),986);
    ex13 << tablepath << "GridSM44.dat";
    array44 = readTable(ex13.str(),986);
    ex14 << tablepath << "GridSM45.dat";
    array45 = readTable(ex14.str(),986);
    ex15 << tablepath << "Grid_X_bb.dat";
    arrayX_bb = readTable(ex15.str(),199);
    ex16 << tablepath << "Grid_X_tt.dat";
    arrayX_tt = readTable(ex16.str(),198);
}    
Exemplo n.º 20
0
// This is the first prototype version of HSpiceRead function for reading HSpice
// output files.
// TODO:
//   ascii format support
//   different vector types support (like voltage, current ..., although I do not
//                                   know what it would be good for)
//   scale monotonity check
static PyObject *HSpiceRead(PyObject *self, PyObject *args)
{
    const char *fileName;
    char *token, *buf = NULL, **name = NULL;
    int debugMode, num, numOfVectors, numOfVariables, type, sweepSize = 1,
        i = dateStartPosition - 1, offset = 0;
    struct FastArray faSweep, *faPtr = NULL;

    int postVersion = 0;

    FILE *f = NULL;
    PyObject *date = NULL, *title = NULL, *scale = NULL, *sweep = NULL,
             *sweepValues = NULL, *dataList = NULL, **tmpArray = NULL, *sweeps = NULL,
             *tuple = NULL, *list = NULL;

    // Get hspice_read() arguments.
    if(!PyArg_ParseTuple(args, "si", &fileName, &debugMode)) return Py_None;

    if(debugMode) fprintf(debugFile, "HSpiceRead: reading file %s.\n", fileName);

    f = fopen(fileName, "rb");    // Open the file.
    if(f == NULL)
    {
        if(debugMode)
            fprintf(debugFile, "HSpiceRead: cannot open file %s.\n", fileName);
        goto failed;
    }

    num = getc(f);
    ungetc(num, f);
    if(num == EOF)    // Test if there is data in the file.
    {
        if(debugMode)
            fprintf(debugFile, "HSpiceRead: file %s is empty.\n", fileName);
        goto failed;
    }
    if((num & 0x000000ff) >= ' ')    // Test if the file is in binary format.
    {
        if(debugMode) fprintf(debugFile,
                "HSpiceRead: file %s is in ascii format.\n",
                fileName);
        goto failed;
    }

    // Read file header blocks.
    do num = readHeaderBlock(f, debugMode, fileName, &buf, &offset);
    while(num == 0);
    if(num > 0) goto failed;

    if(strncmp(&buf[postStartPosition1], "9007", numOfPostCharacters) == 0)
        postVersion = 9007;
    else if(strncmp(&buf[postStartPosition1], "9601", numOfPostCharacters) == 0)
        postVersion = 9601;
    else if(strncmp(&buf[postStartPosition2], "2001", numOfPostCharacters) == 0)
        postVersion = 2001;
    else
    {
        if(debugMode) fprintf(debugFile, "HSpiceRead: unknown post format.\n");
        goto failed;
    }

    if(debugMode) fprintf(debugFile, "POST_VERSION=%d\n", postVersion);

    buf[dateEndPosition] = 0;
    date = PyString_FromString(&buf[dateStartPosition]);    // Get creation date.
    if(date == NULL)
    {
        if(debugMode)
            fprintf(debugFile, "HSpiceRead: failed to create date string.\n");
        goto failed;
    }

    while(buf[i] == ' ') i--;
    buf[i + 1] = 0;
    title = PyString_FromString(&buf[titleStartPosition]);    // Get title.
    if(title == NULL)
    {
        if(debugMode)
            fprintf(debugFile, "HSpiceRead: failed to create title string.\n");
        goto failed;
    }

    buf[numOfSweepsEndPosition] = 0;    // Check number of sweep parameters.
    num = atoi(&buf[numOfSweepsPosition]);
    if(num < 0 || num > 1)
    {
        if(debugMode) fprintf(debugFile,
                "HSpiceRead: only onedimensional sweep supported.\n");
        goto failed;
    }

    buf[numOfSweepsPosition] = 0;    // Get number of vectors (variables and probes).
    numOfVectors = atoi(&buf[numOfProbesPosition]);
    buf[numOfProbesPosition] = 0;
    numOfVariables = atoi(&buf[numOfVariablesPosition]);    // Scale included.
    numOfVectors = numOfVectors + numOfVariables;

    // Get type of variables. Scale is always real.
    token = strtok(&buf[vectorDescriptionStartPosition], " \t\n");
    type = atoi(token);
    if(type == frequency) type = complex_var;
    else type = real_var;

    for(i = 0; i < numOfVectors; i++) token = strtok(NULL, " \t\n");
    if(token == NULL)
    {
        if(debugMode) fprintf(debugFile,
                "HSpiceRead: failed to extract independent variable name.\n");
        goto failed;
    }

    scale = PyString_FromString(token);    // Get independent variable name.
    if(scale == NULL)
    {
        if(debugMode) fprintf(debugFile,
                "HSpiceRead: failed to create independent variable name string.\n");
        goto failed;
    }

    // Allocate space for pointers to vector names.
    name = (char **)PyMem_Malloc((numOfVectors - 1) * sizeof(char *));
    if(name == NULL)
    {
        if(debugMode) fprintf(debugFile,
                "HSpiceRead: cannot allocate pointers to vector names.\n");
        goto failed;
    }

    for(i = 0; i < numOfVectors - 1; i++)    // Get vector names.
    {
        name[i] = strtok(NULL, " \t\n");
        if(name[i] == NULL)
        {
            if(debugMode) fprintf(debugFile,
                    "HSpiceRead: failed to extract vector names.\n");
            goto failed;
        }
    }

    // Process vector names: make name lowercase, remove v( in front of name
    for(i=0; i < numOfVectors - 1; i++) {
        int j;
        for(j=0;name[i][j];j++) {
            if (name[i][j]>='A' && name[i][j]<='Z') {
                name[i][j]-='A'-'a';
            }
        }
        if (name[i][0]=='v' && name[i][1]=='(') {
            for(j=2;name[i][j];j++) {
                name[i][j-2]=name[i][j];
            }
            name[i][j-2]=0;
        }
    }

    if(num == 1)    // Get sweep information.
    {
        int num = getSweepInfo(debugMode, postVersion,
                &sweep, buf, &sweepSize, &sweepValues,
                &faSweep);
        if(num) goto failed;
    }

    dataList = PyList_New(0);    // Create an empty list for data dictionaries.
    if(dataList == NULL)
    {
        if(debugMode)
            fprintf(debugFile, "HSpiceRead: failed to create data list.\n");
        goto failed;
    }

    // Allocate space for pointers to arrays.
    tmpArray = (PyObject **)PyMem_Malloc(numOfVectors * sizeof(PyObject *));
    if(tmpArray == NULL)
    {
        if(debugMode)
            fprintf(debugFile, "HSpiceRead: cannot allocate pointers to arrays.\n");
        goto failed;
    }

    // Allocate space for fast array pointers.
    faPtr =
        (struct FastArray *)PyMem_Malloc(numOfVectors * sizeof(struct FastArray));
    if(faPtr == NULL)
    {
        if(debugMode) fprintf(debugFile, "HSpiceRead: failed to create array.\n");
        goto failed;
    }

    if(debugMode) fprintf(debugFile,"numOfVectors=%d\n", numOfVectors);

    for(i = 0; i < sweepSize; i++)    // Read i-th table.
    {
        num = readTable(f, debugMode, postVersion,
                fileName, sweep, numOfVariables, type,
                numOfVectors, &faSweep, tmpArray, faPtr, token, name,
                dataList);
        if(num) goto failed;
    }
    fclose(f);
    f = NULL;
    PyMem_Free(faPtr);
    faPtr = NULL;
    PyMem_Free(buf);
    buf = NULL;
    PyMem_Free(name);
    name = NULL;
    PyMem_Free(tmpArray);
    tmpArray = NULL;

    // Create sweeps tuple.
    if(sweep == NULL) sweeps = PyTuple_Pack(3, Py_None, Py_None, dataList);
    else sweeps = PyTuple_Pack(3, sweep, sweepValues, dataList);
    if(sweeps == NULL)
    {
        if(debugMode) fprintf(debugFile,
                "HSpiceRead: failed to create tuple with sweeps.\n");
        goto failed;
    }
    Py_XDECREF(sweep);
    Py_XDECREF(sweepValues);
    Py_XDECREF(dataList);

    // Prepare return tuple.
    tuple = PyTuple_Pack(6, sweeps, scale, Py_None, title, date, Py_None);
    if(tuple == NULL)
    {
        if(debugMode) fprintf(debugFile,
                "HSpiceRead: failed to create tuple with read data.\n");
        goto failed;
    }
    Py_XDECREF(date);
    date = NULL;
    Py_XDECREF(title);
    title = NULL;
    Py_XDECREF(scale);
    scale = NULL;
    Py_XDECREF(sweeps);
    sweeps = NULL;

    list = PyList_New(0);    // Create an empty list.
    if(list == NULL)
    {
        if(debugMode)
            fprintf(debugFile, "HSpiceRead: failed to create return list.\n");
        goto failed;
    }

    num = PyList_Append(list, tuple);    // Insert tuple into return list.
    if(num)
    {
        if(debugMode) fprintf(debugFile,
                "HSpiceRead: failed to append tuple to return list.\n");
        goto failed;
    }
    Py_XDECREF(tuple);

    return list;

failed:    // Error occured. Close open file, relese memory and python references.
    if (f)
        fclose(f);
    PyMem_Free(buf);
    Py_XDECREF(date);
    Py_XDECREF(title);
    Py_XDECREF(scale);
    PyMem_Free(name);
    Py_XDECREF(sweep);
    Py_XDECREF(sweepValues);
    Py_XDECREF(dataList);
    PyMem_Free(tmpArray);
    PyMem_Free(faPtr);
    Py_XDECREF(sweeps);
    Py_XDECREF(tuple);
    Py_XDECREF(list);
    return Py_None;
}
Exemplo n.º 21
0
int
videoGetRects(lua_State *L, int index, Array *rects)
{
    return readTable(L, index, rects, sizeof (SDL_Rect), readRects);
}
Exemplo n.º 22
0
bool Archive::open( const char * archivePath, Access access )
{
	lock();

	// close this archive if previously open
	close();

	// set the access
	m_Access = access;
	// attempt to open the file
	if (! m_pFile->open( archivePath, access == READ ? FileDisk::READ : FileDisk::READ_WRITE ) || !readTable() )
	{
		unlock();
		return false;	// failed
	}

	m_Open = true;

	unlock();

	// success
	return true;
}
Exemplo n.º 23
0
int
videoGetPoints(lua_State *L, int index, Array *points)
{
    return readTable(L, index, points, sizeof (SDL_Point), readPoints);
}
Exemplo n.º 24
0
Decoder::Decoder(const char filename[], double prune_threshold, unsigned int prune_count):
	flex(new Lexicon(french)), elex(new Lexicon(english)), schwarz(new PTree<PTree<Cost>>){
	readTable(filename, prune_threshold, prune_count);
}