matrix-con.cpp

Go to the documentation of this file.

#include "matrix-con.hpp"
#include "matrix.hpp"
#include "debug.hpp"
 
#include <iostream>
 
MatrixConstructor::MatrixConstructor()
    : R(nullptr), rows(nullptr), cols(nullptr), cols_frozen(false), deg(nullptr)
{
}
 
MatrixConstructor::MatrixConstructor(const FreeModule *target, int ncols)
    : R(target->get_ring()), rows(target), cols(nullptr), cols_frozen(false), deg(nullptr)
{
  cols = R->make_FreeModule(ncols);
 
  entries.reserve(ncols);
  for (int i = 0; i < ncols; i++) entries.push_back(nullptr);
 
  deg = R->degree_monoid()->make_one();
}
 
MatrixConstructor::MatrixConstructor(const FreeModule *target,
                                     const FreeModule *source,
                                     const_monomial deg0)
    : R(target->get_ring()), rows(target), cols(source), cols_frozen(true)
{
  entries.reserve(source->rank());
  for (int i = 0; i < source->rank(); i++) entries.push_back(nullptr);
 
  deg = (deg0 == nullptr ? R->degree_monoid()->make_one()
                   : R->degree_monoid()->make_new(deg0));
}
 
void MatrixConstructor::append(vec v)
{
  assert(!cols_frozen);
  monomial d = R->degree_monoid()->make_one();
  if (v != nullptr) R->vec_multi_degree(rows, v, d);
  append(v, d);
  R->degree_monoid()->remove(d);
}
void MatrixConstructor::append(vec v, const_monomial deg0)
{
  if (cols_frozen)
    {
      INTERNAL_ERROR("trying to append to an immutable free module.");
    }
  entries.push_back(v);
  FreeModule *mutable_cols = const_cast<FreeModule *>(cols);
  mutable_cols->append(deg0);
}
 
void MatrixConstructor::set_entry(int r, int c, ring_elem a)
{
  R->set_entry(entries[c], r, a);
}
void MatrixConstructor::set_column(int c, vec v) { entries[c] = v; }
void MatrixConstructor::compute_column_degrees()
/* Takes into account the matrix degree */
{
  if (cols_frozen)
    {
      INTERNAL_ERROR("trying to append to an immutable free module.");
    }
  for (int i = 0; i < cols->rank(); i++) compute_column_degree(i);
}
 
void MatrixConstructor::set_column_degree(int i, const_monomial deg0)
{
  if (cols_frozen)
    {
      INTERNAL_ERROR("trying to append to an immutable free module.");
    }
  FreeModule *mutable_cols = const_cast<FreeModule *>(cols);
  mutable_cols->change_degree(i, deg0);
}
 
void MatrixConstructor::compute_column_degree(int i)
{
  if (cols_frozen)
    {
      INTERNAL_ERROR("trying to append to an immutable free module.");
    }
  monomial d = R->degree_monoid()->make_one();
  const vec v = entries[i];
  if (v != nullptr) R->vec_multi_degree(rows, v, d);
  FreeModule *mutable_cols = const_cast<FreeModule *>(cols);
  mutable_cols->change_degree(i, d);
  R->degree_monoid()->remove(d);
}
 
void MatrixConstructor::set_matrix_degree(const_monomial deg0) { deg = deg0; }
Matrix *MatrixConstructor::to_matrix()
{
  //  if (!will_be_mutable && !cols->is_frozen)
  //    cols->freeze(hashval);
  return new Matrix(rows, cols, deg, entries);
}
 
void MatrixConstructor::debugDisplay() const
{
  std::cout << "MatrixConstructor: ring = " << R << std::endl;
  std::cout << "MatrixConstructor: rows = " << rows << std::endl;
  std::cout << "MatrixConstructor: cols = " << cols << std::endl;
  std::cout << "Matrixconstructor: #entries = " << entries.size() << std::endl;
}
 
// Local Variables:
// compile-command: "make -C $M2BUILDDIR/Macaulay2/e "
// indent-tabs-mode: nil
// End: