1 files changed, 331 insertions, 0 deletions
diff --git a/rasodmg/test/test_bmark_stat.cc b/rasodmg/test/test_bmark_stat.cc
new file mode 100644
index 0000000..2c5b776
--- /dev/null
+++ b/rasodmg/test/test_bmark_stat.cc
@@ -0,0 +1,331 @@
+/*
+* 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_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 <iostream>
+
+#include <stdio.h>
+#include <stdlib.h>
+#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<r_Access>  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<r_ULong> > > >  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<r_ULong> > 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<RGBPixel> > >;
+
+      db.set_object_name(*cube_set, colect_name);
+    } 
+
+    cout << "Ok" << endl;
+    cout << "Creating the datacube... " << flush;
+
+    cube =
+        new(&db, "RGB_3D_Cube") r_Marray<RGBPixel>(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;
+}
+
+
+
+
+
+
+
+
+
+
+