diff options
Diffstat (limited to 'rasodmg/test/test_bmark_pet.cc')
| -rw-r--r-- | rasodmg/test/test_bmark_pet.cc | 235 |
1 files changed, 235 insertions, 0 deletions
diff --git a/rasodmg/test/test_bmark_pet.cc b/rasodmg/test/test_bmark_pet.cc new file mode 100644 index 0000000..8fa94fb --- /dev/null +++ b/rasodmg/test/test_bmark_pet.cc @@ -0,0 +1,235 @@ +/* +* 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_pet.cc + * + * MODULE: rasodmg + * + * PURPOSE: benchmark interest tiling + * + * COMMENTS: + * This program is used to load the database with information + * for benchmarking interesting tiling + * +*/ + +#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 "rasodmg/odmgtypes.hh" +#include "rasodmg/storagelayout.hh" +#include "rasodmg/tiling.hh" +#include "rasodmg/dirtiling.hh" +#include "rasodmg/dirdecompose.hh" +#include "rasodmg/storagelayout.hh" +#include "raslib/oid.hh" +#include "include/basictypes.hh" + +#define S_32K (32 * 1024L) +#define S_64K (64 * 1024L) +#define S_128K (128 * 1024L) +#define S_256K (256 * 1024L) + +#define TOTAL_CUBES 7 +#define SIZE_X 185L +#define SIZE_Y 150L +#define SIZE_Z 141L + +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]; +} + +void insert_datacube( ) +{ + + r_Ref< r_Set< r_Ref< r_Marray<r_Short> > > > cube_set; + r_Minterval domain[TOTAL_CUBES]; + r_Domain_Storage_Layout* dsl[TOTAL_CUBES]; + r_OId oid[TOTAL_CUBES]; + + for (int i = 0; i < TOTAL_CUBES; i++) + { + domain[i] = r_Minterval(3); + domain[i] << r_Sinterval(0L, SIZE_X - 1 ) + << r_Sinterval(0L, SIZE_Y - 1) + << r_Sinterval(0L, SIZE_Z - 1); + } + + // For aligned 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); + + 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); + r_Aligned_Tiling* til_reg_64k1 = new r_Aligned_Tiling(block_config, S_64K); + + + // For areas of interest tiling + r_Minterval interest1_1("[23:170,4:75,12:138]"); + r_Minterval interest1_2("[61:149,70:135,8:130]"); + r_Minterval interest2_1("[28:175,2:73,8:134]"); + r_Minterval interest2_2("[63:151,75:140,11:133]"); + + + DList<r_Minterval> areas1; + areas1 += interest1_1; + areas1 += interest1_2; + DList<r_Minterval> areas2; + areas2 += interest2_1; + areas2 += interest2_2; + + r_Interest_Tiling* til_int1 = + new r_Interest_Tiling(areas1, r_Interest_Tiling::NO_LIMIT); + r_Interest_Tiling* til_int2 = + new r_Interest_Tiling(areas2, r_Interest_Tiling::NO_LIMIT); + + // Domain storage layouts + + dsl[0] = new r_Domain_Storage_Layout(domain[0], til_reg_32k); + dsl[1] = new r_Domain_Storage_Layout(domain[1], til_reg_64k); + dsl[2] = new r_Domain_Storage_Layout(domain[2], til_reg_128k); + dsl[3] = new r_Domain_Storage_Layout(domain[3], til_reg_256k); + + dsl[4] = new r_Domain_Storage_Layout(domain[4], til_int1); + dsl[5] = new r_Domain_Storage_Layout(domain[5], til_int2); + dsl[6] = new r_Domain_Storage_Layout(domain[6], til_reg_64k1); + + + // Create cubes + + r_Database db; + db.set_servername(server_name); + + for ( i=0; i< TOTAL_CUBES ; i++) + { + r_Transaction trans; + + r_Ref< r_Marray<r_Short> > 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, "ShortSet3") r_Set< r_Ref< r_Marray<r_Short> > >; + + db.set_object_name(*cube_set, colect_name); + } + + cout << "Ok" << endl; + cout << "Creating the datacube... " << flush; + r_Minterval newDomain( domain[i]); + cube = + new(&db, "ShortCube") r_Marray<r_Short>(newDomain, dsl[i]); + + cube_set->insert_element(cube); + + cout << "Cube[" << i+1 << "]: " << cube->get_oid() << endl; + cout << "Spatial domain: " << cube->spatial_domain( ) <<endl; + cout << "Storage Layout " << endl; + dsl[i]->print_status( ); + + cout << "*" << flush; + cout << " ... Ok" << endl; + cout << "Commiting transaction... " << flush; + + trans.commit(); + + cout << "Ok" << endl; + // cout << "Destroying cube... " <<flush; + // cube.destroy( ); + cout << "Closing database... " << flush; + db.close(); + } + 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); + insert_datacube( ); + return 0; +} + + + + + + + + |