#! /bin/bash -e
# SPDX-License-Identifier: EPL-2.0 OR GPL-2.0-or-later
# SPDX-FileCopyrightText: Bradley M. Bell <bradbell@seanet.com>
# SPDX-FileContributor: 2003-22 Bradley M. Bell
# ----------------------------------------------------------------------------
gigabytes='1.0'   # memory limit in gigabytes
# -----------------------------------------------------------------------------
kilobytes=`echo "($gigabytes * 10^9) / 1024" | bc`
ulimit -Sv $kilobytes
cat << EOF
This is testing the idea that
   vector< std::set<size_t> > sparsity;
is inefficient compared to
   vector<size_t> row, col;
The results on one system for this script are:
   ./sparsity set 1e7 elapsed_seconds = 3.62206
   ./sparsity set 2e7 std::bad_alloc
   ./sparsity vec 2e7 elapsed_seconds = 1.34243
   ./sparsity vec 4e7 std::bad_alloc
EOF
cat << EOF > bug.$$
# include <cppad/cppad.hpp>

int main(int argc, char *argv[])
{  using CppAD::elapsed_seconds;
   if( argc != 3 )
   {  std::cerr << "usage: $0 (set|vec) n" << std::endl;
      return 1;
   }
   bool set  = std::strcmp(argv[1], "set") == 0;
   bool vec  = std::strcmp(argv[1], "vec") == 0;
   bool ok   = vec || set;
   if( ! ok )
   {  std::cerr << "usage: $0 (set|vec) n" << std::endl;
      return 1;
   }
   size_t n = size_t( std::atof( argv[2] ) );
   const char* label;
   elapsed_seconds();
   try
   {
      if( set )
      {  std::vector< std::set<size_t> > my_set(n);
         for(size_t i = 0; i < n; i++)
            my_set[i].insert(i);
      }
      else
      {  std::vector<size_t> row;
         std::vector<size_t> col;
         for(size_t i = 0; i < n; i++)
         {  row.push_back(i);
            col.push_back(i);
         }
      }
      for(int i = 0; i < argc; i++)
         std::cout << argv[i] << " ";
      std::cout << "elapsed_seconds = " << elapsed_seconds() << std::endl;
   }
   catch( std::bad_alloc& ba )
   {
      for(int i = 0; i < argc; i++)
         std::cout << argv[i] << " ";
      std::cout << ba.what() << std::endl;
   }
   return 0;
}
EOF
# -----------------------------------------------------------------------------
if [ ! -e build ]
then
   mkdir build
fi
cd build
echo "$0"
name=`echo $0 | sed -e 's|.*/||' -e 's|\..*||'`
mv ../bug.$$ $name.cpp
echo "g++ -I../.. -DNDEBUG --std=c++11 $name.cpp -o $name"
g++ -I../.. -DNDEBUG --std=c++11 $name.cpp -o $name
#
./$name set 1e7
./$name set 2e7
./$name vec 2e7
./$name vec 4e7
