diff options
Diffstat (limited to 'rasodmg/test/test_bmark_dir1.cc')
-rw-r--r-- | rasodmg/test/test_bmark_dir1.cc | 318 |
1 files changed, 318 insertions, 0 deletions
diff --git a/rasodmg/test/test_bmark_dir1.cc b/rasodmg/test/test_bmark_dir1.cc new file mode 100644 index 0000000..c185858 --- /dev/null +++ b/rasodmg/test/test_bmark_dir1.cc @@ -0,0 +1,318 @@ +/* +* 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 <http://www.gnu.org/licenses/>. +* +* Copyright 2003, 2004, 2005, 2006, 2007, 2008, 2009 Peter Baumann / +rasdaman GmbH. +* +* For more information please see <http://www.rasdaman.org> +* or contact Peter Baumann via <baumann@rasdaman.com>. +/ +/** + * SOURCE: test_bmark_dir1.cc + * + * MODULE: rasodmg + * + * PURPOSE: benchmark dirtiling + * + * COMMENTS: + * This program creates 3D and 5D datacubes for benchmarking + * the directional tiling method + * +*/ + +#include <iostream> +#include <stdio.h> +#include <stdlib.h> +#include "rasodmg/ref.hh" +#include "rasodmg/transaction.hh" +#include "rasodmg/database.hh" +#include "rasodmg/set.hh" +#include "rasodmg/marray.hh" +#include "raslib/odmgtypes.hh" +#include "rasodmg/tiling.hh" +#include "rasodmg/dirtiling.hh" +#include "rasodmg/dirdecompose.hh" +#include "rasodmg/storagelayout.hh" +#include "raslib/oid.hh" + +#define YEARS (2L) + +#define TOTAL_DAYS (365L * 2L) +#define TOTAL_PRODUCTS (60L) +#define TOTAL_STORES (100L) + + +#define MAX_SALES 500 + +#define S_32K (32 * 1024L) +#define S_64K (64 * 1024L) +#define S_128K (128 * 1024L) +#define S_256K (256 * 1024L) + +#define TOTAL_CUBES 10 + + +char* server_name; +char* dbase_name; +char* colect_name; + +void parse(int argc, char* argv[]) +{ + if (argc != 4) + { + cout << "Usage: " << argv[0] << " [server name] [db name] [colection name]" + << endl; + + exit(0); + } + + server_name = argv[1]; + dbase_name = argv[2]; + colect_name = argv[3]; +} + +r_ULong init(const r_Point& pnt ) +{ + return (long)(rand() % MAX_SALES); +} + +void insert_datacube() +{ + + r_Ref< r_Set< r_Ref< r_Marray<r_ULong> > > > cube_set; + // r_Ref< r_Marray<r_ULong> > cube[TOTAL_CUBES]; + r_Minterval domain, block_config; + r_Domain_Storage_Layout* dsl[TOTAL_CUBES]; + r_OId oid[TOTAL_CUBES]; + + domain = r_Minterval(3); + domain << r_Sinterval(1L, TOTAL_DAYS) + << r_Sinterval(1L, TOTAL_PRODUCTS) + << r_Sinterval(1L, TOTAL_STORES); + + + block_config = r_Minterval(3); + block_config << r_Sinterval(0L, TOTAL_DAYS) + << r_Sinterval(0L, TOTAL_PRODUCTS) + << r_Sinterval(0L, TOTAL_STORES); + + // Each storage object must have an own dynamic tiling obj or else the client + // ( and server ) crashes because memory is released for a non heap memory free + // r_Marray become responsible for managing the memory allocated for the + // tiling object. + r_Aligned_Tiling* til_reg_32k = new r_Aligned_Tiling(block_config, S_32K); + r_Aligned_Tiling* til_reg_64k = new r_Aligned_Tiling(block_config, S_64K); + r_Aligned_Tiling* til_reg_128k = new r_Aligned_Tiling(block_config, S_128K); + r_Aligned_Tiling* til_reg_256k = new r_Aligned_Tiling(block_config, S_256K); + + + // For directional tiling + + r_Dir_Decompose decomp[3]; + + // Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dez + unsigned int daysMonth[] = { 31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31}; + int ix = 0; + int year; + int month; + decomp[0] << 1; + for ( year = 0; year < YEARS ; year++) + { + for (month = 0; month < 12; month++ ) + { + decomp[0] << ix + daysMonth[month]; + ix += daysMonth[month]; + } + } + + decomp[1] << 1 << 60; + + // << 1 << 1 + 26 << 27 + 8 << 35 + 6 << 41 + 18 << 59 + 14 << 73 + 16 << 89 + 8 << 97 + 3; + decomp[2] << 1 << 27 << 35 << 41 << 59 << 73 << 89 << 97 << 100; + + + r_Dir_Tiling* til_dir_32k = new r_Dir_Tiling(3, decomp, S_32K); + r_Dir_Tiling* til_dir_64k = new r_Dir_Tiling(3, decomp, S_64K); + r_Dir_Tiling* til_dir_128k = new r_Dir_Tiling(3, decomp, S_128K); + r_Dir_Tiling* til_dir_256k = new r_Dir_Tiling(3, decomp, S_256K); + + + r_Dir_Decompose decomp1[3]; + + ix = 0; + decomp1[0] << 1; + for ( year = 0; year < YEARS ; year++) + { + for (month = 0; month < 12; month++ ) + { + decomp1[0] << ix + daysMonth[month]; + ix += daysMonth[month]; + } + } + + // Products + // << 1<< 1+ 26 << 27 + 15 << 42 + 28 + decomp1[1] << 1 << 27 << 42 << 60; + + // Stores + // << 1 << 1 + 26 << 27 + 8 << 35 + 6 << 41 + 18 << 59 + 14 << 73 + 16 << 89 + 8 << 97 + 3; + decomp1[2] << 1 << 27 << 35 << 41 << 59 << 73 << 89 << 97 << 100; + + r_Dir_Tiling* til_dir1_32k = new r_Dir_Tiling(3, decomp1, S_32K); + r_Dir_Tiling* til_dir1_64k = new r_Dir_Tiling(3, decomp1, S_64K); + + // Domain storage layouts + + dsl[0] = new r_Domain_Storage_Layout(domain, til_reg_32k); + dsl[1] = new r_Domain_Storage_Layout(domain, til_reg_64k); + dsl[2] = new r_Domain_Storage_Layout(domain, til_reg_128k); + dsl[3] = new r_Domain_Storage_Layout(domain, til_reg_256k); + + dsl[4] = new r_Domain_Storage_Layout(domain, til_dir_32k); + dsl[5] = new r_Domain_Storage_Layout(domain, til_dir_64k); + dsl[6] = new r_Domain_Storage_Layout(domain, til_dir_128k); + dsl[7] = new r_Domain_Storage_Layout(domain, til_dir_256k); + + dsl[8] = new r_Domain_Storage_Layout(domain, til_dir1_32k); + dsl[9] = new r_Domain_Storage_Layout(domain, til_dir1_64k); + + +for (int i= 0 ; i< TOTAL_CUBES ; i++) +{ + r_Database db; + r_Transaction trans; + // The main phase of the database creation + r_Ref< r_Marray<r_ULong> > cube1; + + db.set_servername(server_name); + + 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, "ULong_3D_Set") r_Set< r_Ref< r_Marray<r_ULong> > >; + db.set_object_name(*cube_set, colect_name); + } + + cout << "Ok" << endl; + cout << "Creating the datacube... " << flush; + + // for (int i=0; i<TOTAL_CUBES; i++) + // { + // cube[i] = + cout << "domain == " << domain << endl; + cube1 = + new(&db, "ULong_3D_Cube") r_Marray<r_ULong>(domain, 1L/* &init */, dsl[i]); + + // cube_set->insert_element(cube[i]); + // oid[i] = cube[i].get_oid(); + cube_set->insert_element(cube1); + oid[i] = cube1->get_oid(); + + cout << "*" << flush; + + + cout << " ... Ok" << endl; + + cout << " Cube[" << i+1 << "]: " << oid[i] << endl; + cout << " Spatial domain: " << cube1->spatial_domain( ) << endl; + cout << " Type length: " << cube1->get_type_length( ) << endl; + cout << " Storage Layout: "; + if( i < 4 ) + cout << "regular; tile size " << dsl[i]->get_tile_size( ) << endl; + else + { + cout << "directional; tile size " << dsl[i]->get_tile_size( ) << endl; + cout << "Dir decompose: "; + for ( int j = 0; j < 3 ; j++ ) + { + for ( int k = 0; k < 3 ; k++) + { + if ( i < 8 ) + decomp[k].print_status(cout ); + else + decomp1[k].print_status( cout ); + } + } + cout << endl; + } + + cout << "Commiting transaction... " << flush; + + trans.commit(); + + cout << "Ok" << endl; + + cout << "Closing database... " << flush; + + db.close(); + + } + catch (r_Error& e) + { + cout << e.what() << endl; + exit(0); + } + catch (...) + { + cout << "Undefined error..." << endl; + exit(0); + } +} + + cout << "Ok [******************]" << endl << flush; + + cout << endl; + cout << "Inserted data resume" << endl; + cout << "====================" << endl; +} + +int main(int argc, char* argv[]) +{ + parse(argc, argv); + insert_datacube(); + + return 0; +} + + + + + + + + |