src/blad_session.cpp

Go to the documentation of this file.

/******************************************************************************
* MODULE     : blad_session.cpp
* DESCRIPTION: wrapper of the blad library
* COPYRIGHT  : (C) 2012  Gregoire Lecerf
*******************************************************************************
* This software falls under the GNU general public license and comes WITHOUT
* ANY WARRANTY WHATSOEVER. See the file $TEXMACS_PATH/LICENSE for more details.
* If you don't have this file, write to the Free Software Foundation, Inc.,
* 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
******************************************************************************/

#include <basix/mmx_syntax.hpp>
#include <mblad/bav.hpp>
#include <bad.h>

namespace mmx {
#define Coordinate multivariate_coordinate<>

static bav_symbol*
bav_independent_symbol (const string& s, nat i, nat j) {
  bav_symbol* x= ::bav_new_symbol ();
  char* c= as_charp (s);
  x -> ident= ::ba0_strdup (c); 
  free_charp (c);
  x -> type= ::bav_independent_symbol;
  x -> index= i;
  x -> derivation_index= j;
  return x; }

static bav_symbol*
bav_dependent_symbol (const string& s, nat i) {
  bav_symbol* x= ::bav_new_symbol ();
  char* c= as_charp (s);
  x -> ident= ::ba0_strdup (c); 
  free_charp (c);
  x -> type= ::bav_dependent_symbol;
  x -> index= i;
  //x -> derivation_index= -1; // FIXME ???
  return x; }

static bav_variable*
bav_independent_variable (const string& s, nat i, nat j) {
  bav_variable* x= ::bav_new_variable ();
  x -> root= bav_independent_symbol (s, i, j);
  x -> index= i;
  ba0_insert_table ((ba0_table*) &(x -> number), 0, (void*) 2);
  return x; }

static bav_variable*
bav_dependent_variable (const string& s, nat i, nat r) {
  bav_variable* x= ::bav_new_variable ();
  x -> root= bav_dependent_symbol (s, i);
  x -> index= i;
  ba0_insert_table ((ba0_table*) &(x -> number), 0, (void*) 2);
  for (nat k= 0; k < r; k++)
    ba0_insert_table ((ba0_table*) &(x -> order), k, (void*) 0);
  return x; }

blad_session::blad_session (const vector<Coordinate>& indeps2,
                            const vector<Coordinate>& deps2)
  : indeps (indeps2), deps (deps2) {
  bad_restart (0, 0);
  bav_R_init ();
  
  nat r= N(indeps), s= N(deps);
  bav_ordering* ord= ::bav_new_ordering ();
  bav_block* blk= ::bav_new_block ();
  bav_variable* var;
  for (nat i= 0; i < r; i++) {
    var= bav_independent_variable (as_mmx (**indeps[i]), i, i);
    bav_indeps << (void*) var;
    ba0_insert_table ((ba0_table*) &(bav_R.strs), i,
                      (void*) var -> root -> ident);
    ba0_insert_table ((ba0_table*) &(bav_R.syms), i,
                      (void*) var -> root);
    ba0_insert_table ((ba0_table*) &(bav_R.vars), i,
                      (void*) var);
    ba0_insert_table ((ba0_table*) &(bav_R.ders), i,
                      (void*) i);
    ba0_insert_table ((ba0_table*) &(ord -> ders), i,
                      (void*) var-> root);
  }
  for (nat i= 0; i < s; i++) {
    var= bav_dependent_variable (as_mmx (**deps[i]), r + i, r);
    bav_deps << (void*) var;
    ba0_insert_table ((ba0_table*) &(bav_R.strs), r + i,
                      (void*) var -> root -> ident);
    ba0_insert_table ((ba0_table*) &(bav_R.syms), r + i,
                      (void*) var -> root);
    ba0_insert_table ((ba0_table*) &(bav_R.vars), r + i,
                      (void*) var);
    ba0_insert_table ((ba0_table*) &(bav_R.deps), i, (void*) (r + i));
    ba0_insert_table ((ba0_table*) &(blk -> indices), i, (void*) (r + i));
    ba0_insert_table ((ba0_table*) &(blk -> strs), i,
                      (void*) var -> root -> ident);
  }
  ba0_insert_table ((ba0_table*) &(ord -> blocks), 0, (void*) blk);
  ba0_insert_table ((ba0_table*) &(bav_R.ords), 0, (void*) ord);
  ba0_insert_table ((ba0_table*) &(bav_R.ord_stack), 0, (void*) 0);
  bav_R.opra= 0;
  bav_R.empty= false;
  // ba0_printf ((char*) "current_ordering= %ordering\n", 0);
}
  
blad_session::~blad_session () {
  // mmout << "deleting blad session" << lf;
  bad_terminate (ba0_init_level); }

} // namespace mmx

---