derive.H

#ifndef _DERIVE_H_
#define _DERIVE_H_

#include "algebras.H"
#include "multipoly.H"
#include "multigrobner.H"
#include <map>
#include <string>
#include <sstream>
using namespace std;


class DerivationAnalyzer {
public:
  GradedAlgebra<F_2> *source;
  GradedAlgebra<F_2> squares;
  AffineHomomorphism<F_2> sq_inclusion;
  Tuple< multiPolynomial<F_2> > d;
  GradedAlgebra<F_2> constants;
  AffineHomomorphism<F_2> cs_inclusion;


  Tuple< Polynomial<F_2> > mod_gens;
  map< PowProd, long > index_of;

  void setup(GradedAlgebra<F_2> *thesource);

  void write_gens(list< Polynomial<F_2> >& result, 
                  long ext, long first, long last);

  virtual multiPolynomial<F_2> coords_of(PowProd P);
  virtual multiPolynomial<F_2> coords_of(const Polynomial<F_2>& P);

  void constant_subalgebra(string prefix, long offset);
  virtual ~DerivationAnalyzer(){};
};


class MilnorAnalyzer : public DerivationAnalyzer {
public:
  UnstableAlgebra* A;  
  GradedAlgebra<F_2> ans;

  AffineHomomorphism<F_2> inclusion;

  AffineHomomorphism<F_2> ans_inclusion;


  bool done();

  AbstractHomomorphism<F_2> newvar_squared;

  long n_cs_vars;
  GradedAlgebra<F_2> alg;
  virtual multiPolynomial<F_2> coords_of(PowProd P);
  virtual multiPolynomial<F_2> coords_of(const Polynomial<F_2>& P);
  void start(UnstableAlgebra *A);

  void setup_for_Sq1();

  void re_setup(long k);

};


#endif

---