/*
* This file is part of rasdaman community.
*
* Rasdaman community is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* Rasdaman community is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with rasdaman community. If not, see .
*
* Copyright 2003, 2004, 2005, 2006, 2007, 2008, 2009 Peter Baumann /
rasdaman GmbH.
*
* For more information please see
* or contact Peter Baumann via .
/
/**
* SOURCE: test_bmark_stat.cc
*
* MODULE: rasodmg
*
* PURPOSE: Load database for statistic tiling benchmark
*
* COMMENTS:
* None
*/
/*
ATENTION: The format of the input file for using with this program is:
border_threshold accesses_threshold tilesize domain
access1
access2
accessXXX
...
Example:
50 0.20 1000 [0:799, 0:599]
[10:20, 30:40]
[12:20, 35:39]
[100:300, 300:400]
[120:300, 310:410]
[200:500, 350:500]
[200:510, 350:500]
*/
#include
#include
#include
#include "raslib/oid.hh"
#include "raslib/dlist.hh"
#include "raslib/minterval.hh"
#include "raslib/sinterval.hh"
#include "rasodmg/ref.hh"
#include "rasodmg/set.hh"
#include "rasodmg/marray.hh"
#include "rasodmg/tiling.hh"
#include "rasodmg/database.hh"
#include "rasodmg/odmgtypes.hh"
#include "rasodmg/stattiling.hh"
#include "rasodmg/transaction.hh"
#include "rasodmg/storagelayout.hh"
#include "RGBCube.hh"
#define S_32K (32 * 1024L)
#define S_64K (64 * 1024L)
#define S_128K (128 * 1024L)
#define S_256K (256 * 1024L)
#define TOTAL_CUBES 8
#define SIZE_X 120L
#define SIZE_Y 159L
#define SIZE_Z 119L
char* filename;
char* server_name;
char* dbase_name;
char* colect_name;
int cube_i;
DList stat_info;
unsigned int border_threshold;
double interesting_threshold;
unsigned long tile_size;
r_Minterval* domain;
void parse(int argc, char* argv[])
{
if (argc != 6)
{
cout << "Usage: " << argv[0] << " server db collec ptfile cube" << endl;
cout << endl;
cout << " server ............. server to use" << endl;
cout << " db ................. database" << endl;
cout << " collec ............. collection" << endl;
cout << " ptfile ............. pattern file" << endl;
cout << " cube ............... datacube to create" << endl;
exit(1);
}
server_name = argv[1];
dbase_name = argv[2];
colect_name = argv[3];
filename = argv[4];
cube_i = atoi(argv[5]);
if ((cube_i<0) || (cube_i>7))
{
cout << "Invalid datacube. Must be in 0..7" << endl;
exit(0);
}
}
void read_data()
{
const unsigned int BUF_SIZE = 256;
char buf[BUF_SIZE], buf2[BUF_SIZE];
cout << "Opening " << filename << " for reading... ";
ifstream is(filename, ios::in);
if (!is)
{
cout << "Couldn't open!!!" << endl;
exit(1);
}
else
cout << "done." << endl;
cout << "Reading parameters... ";
is >> border_threshold;
is >> interesting_threshold;
is >> tile_size;
is.getline(buf, BUF_SIZE);
domain = new r_Minterval(buf);
cout << "done." << endl;
cout << "Geting the accesses... ";
unsigned long count = 0;
while (!is.eof())
{
is.getline(buf, BUF_SIZE);
if (sscanf(buf, "%s", buf2) == 1)
{
r_Minterval inter(buf);
stat_info += inter;
++count;
cout << "*";
}
}
is.close();
cout << endl;
cout << "Geting the accesses... done." << endl << endl;
cout << "Border threshold = " << border_threshold << endl;
cout << "Interesting threshold = " << interesting_threshold << endl;
cout << "Tile size = " << tile_size << endl;
cout << "Domain = " << *domain << endl;
cout << "Number of accesses = " << count << endl << endl;
}
void insert_datacube()
{
r_Ref< r_Set< r_Ref< r_Marray > > > cube_set;
r_Minterval domain;
domain = r_Minterval(3);
domain << r_Sinterval(0L, SIZE_X)
<< r_Sinterval(0L, SIZE_Y)
<< r_Sinterval(0L, SIZE_Z);
// For alligned tiling (Regular tiling)
r_Minterval block_config(3);
block_config << r_Sinterval(0L, SIZE_X)
<< r_Sinterval(0L, SIZE_Y)
<< r_Sinterval(0L, SIZE_Z);
unsigned long ts;
switch (cube_i)
{
case 0:
ts = S_32K;
break;
case 1:
ts = S_64K;
break;
case 2:
ts = S_128K;
break;
case 3:
ts = S_256K;
break;
default: // Irrelevant, not used.
ts = S_32K;
break;
}
r_Aligned_Tiling til_reg(block_config, ts);
r_Stat_Tiling til_stat(border_threshold, interesting_threshold, tile_size);
if (cube_i > 3)
{
til_stat.update_stat_information(stat_info);
}
// Domain storage layout
// This is a hack due to problems with the pointers
r_Domain_Storage_Layout* dsl[2];
r_Domain_Storage_Layout* use;
dsl[0] = new r_Domain_Storage_Layout(domain, &til_reg);
dsl[1] = new r_Domain_Storage_Layout(domain, &til_stat);
if (cube_i<4)
use = dsl[0];
else
use = dsl[1];
// Create cube
r_Database db;
db.set_servername(server_name);
r_Transaction trans;
r_Ref< r_Marray > cube;
try
{
cout << "Opening database " << dbase_name << " on " << server_name
<< "... " << flush;
db.open(dbase_name);
cout << "Ok" << endl;
cout << "Starting transaction... " << flush;
trans.begin();
cout << "Ok" << endl;
cout << "Opening the set... " << flush;
try
{
cube_set = db.lookup_object(colect_name);
}
catch (...)
{
cout << "*Failed*" << endl;
cout << "Creating the set... " << flush;
cube_set =
new(&db, "RGB_3D_Set") r_Set< r_Ref< r_Marray > >;
db.set_object_name(*cube_set, colect_name);
}
cout << "Ok" << endl;
cout << "Creating the datacube... " << flush;
cube =
new(&db, "RGB_3D_Cube") r_Marray(domain, use);
cube_set->insert_element(cube);
cout << "*" << flush;
cout << " ... Ok" << endl;
cout << "Commiting transaction... " << flush;
trans.commit();
cout << "Ok" << endl;
cout << "Closing database... " << flush;
db.close();
cout << " Ok" << endl << flush;
}
catch (r_Error& e)
{
cout << e.what() << endl;
exit(0);
}
catch (...)
{
cout << "Undefined error..." << endl;
exit(0);
}
}
int main(int argc, char* argv[])
{
parse(argc, argv);
read_data();
insert_datacube();
return 0;
}