lattice.cc File Reference

#include "coeffs/bigintmat.h"
#include "lattice.h"
#include "coeffs/numbers.h"
#include "coeffs/coeffs.h"
#include "Singular/ipid.h"
#include <iostream>
#include <vector>
#include <utility>
#include "kernel/numeric/mpr_numeric.h"

Go to the source code of this file.

Macros
#define	DEBUG_PRINTS   1
#define	DEBUG_BLOCK(x)   bool debug_block = x;
#define	DEBUG_PRINT(x)   do {if(debug_block) {Print x ;}} while(0)
#define	DEBUG_CMD(x)   do {if(debug_block) {x;}} while(0)
#define	DEBUG_VAR(x)   do {if(debug_block) {std::cout<<#x<<": "<<x<<std::endl;}} while(0)
#define	DEBUG_N(x)   do {if(debug_block) {Print(#x);Print(": ");n_Print(x,coef);Print("\n");}} while(0)
#define	DEBUG_BIM(x)   do {if(debug_block) {Print(#x);Print(": ");x->Print();Print("\n");}} while(0)

Functions
number	scalarproduct (bigintmat *a, bigintmat *b)
bigintmat *	minkowksi (bigintmat **elementarray, int size_elementarray, number *poly, int deg, coeffs coef, int precision)
bool	IsReal (number a, coeffs coef)
bool	ImagGreaterZero (number a, coeffs coef)
number	squareroot (number a, coeffs coef, int prec)

Variables
STATIC_VAR bool	debug_block = false

Macro Definition Documentation

DEBUG_BIM

#define DEBUG_BIM

(

x

)

do {if(debug_block) {Print(#x);Print(": ");x->Print();Print("\n");}} while(0)

Definition at line 37 of file lattice.cc.

DEBUG_BLOCK

#define DEBUG_BLOCK

(

x

)

bool debug_block = x;

Definition at line 29 of file lattice.cc.

DEBUG_CMD

#define DEBUG_CMD

(

x

)

do {if(debug_block) {x;}} while(0)

Definition at line 34 of file lattice.cc.

DEBUG_N

#define DEBUG_N

(

x

)

do {if(debug_block) {Print(#x);Print(": ");n_Print(x,coef);Print("\n");}} while(0)

Definition at line 36 of file lattice.cc.

DEBUG_PRINT

#define DEBUG_PRINT

(

x

)

do {if(debug_block) {Print x ;}} while(0)

Definition at line 33 of file lattice.cc.

DEBUG_PRINTS

#define DEBUG_PRINTS   1

Definition at line 25 of file lattice.cc.

DEBUG_VAR

#define DEBUG_VAR

(

x

)

do {if(debug_block) {std::cout<<#x<<": "<<x<<std::endl;}} while(0)

Definition at line 35 of file lattice.cc.

Function Documentation

ImagGreaterZero()

bool ImagGreaterZero	(	number	a,
		coeffs	coef
	)

Definition at line 1104 of file lattice.cc.

                                           { //Im(a)>0
    number imag = n_ImPart(a, coef);
    bool out = n_GreaterZero(imag,coef);
    n_Delete(&imag,coef);
    return out;
}

IsReal()

bool IsReal	(	number	a,
		coeffs	coef
	)

Definition at line 1097 of file lattice.cc.

                                  { //Im(a)==0
    number imag = n_ImPart(a, coef);
    bool out = n_IsZero(imag,coef);
    n_Delete(&imag,coef);
    return out;
}

minkowksi()

bigintmat * minkowksi	(	bigintmat **	elementarray,
		int	size_elementarray,
		number *	poly,
		int	deg,
		coeffs	coef,
		int	precision
	)

Definition at line 948 of file lattice.cc.

                                                                                                                         {
    DEBUG_BLOCK(true);
    DEBUG_PRINT(("Begin Minkowski map\n"));
    DEBUG_PRINT(("Input check\n"));
    if(elementarray == NULL || poly == NULL || coef != elementarray[0]->basecoeffs()){
        WerrorS("wrong input!\n");
        return NULL;
    }
 
    for(int i=1;i<size_elementarray;i++){
        if(coef != elementarray[0]->basecoeffs()){
            WerrorS("wrong input!\n");
            return NULL;
        }
    }
 
    //char = 0
    if ( !(nCoeff_is_Z(coef) || nCoeff_is_R(coef) || nCoeff_is_Q(coef) ||
            nCoeff_is_long_R(coef) || nCoeff_is_long_C(coef)) ){
        WerrorS("Ground field not implemented!\n");
        return NULL;
    }
 
    if(deg<2){
        WerrorS("degree of polynomial to small\n");
        return NULL;
    }
    //check and define precision for Q
    if(precision<6){
        precision = 6;
    }
    if ( !(nCoeff_is_R(coef) || nCoeff_is_long_R(coef) || nCoeff_is_long_C(coef)) ){
        setGMPFloatDigits( precision+6,precision+6);
    }
 
    DEBUG_PRINT(("find roots\n"));
    ring CurrOldRing = rCopy(currRing);//need to change currRing, because rootContainer uses the coef of it
    char* n[] = {(char*)"i"};
    ring newring = rDefault(coef, 1, n);
    rChangeCurrRing(newring);
    DEBUG_PRINT(("initialize rootContainer\n"));
    rootContainer * rootcont= new rootContainer();
    rootcont->fillContainer( poly, NULL, 1, deg, rootContainer::onepoly, 1 );///
    rootcont->solver( precision+12);
    int number_roots = rootcont ->getAnzRoots();
    if(number_roots != deg){
        WerrorS("something went wrong: \n\tnot all roots found\n");
        return NULL;
    }
    LongComplexInfo paramComp;
    paramComp.float_len = si_min (precision+6, 32767);
    paramComp.float_len2 = si_min (precision+8, 32767);
    paramComp.par_name=(const char*)"i";
 
    coeffs comp = nInitChar(n_long_C, &paramComp);
 
    number* roots = new number[deg+1];
    number* complexroots = new number[deg+1];
    int r1 = 0;
    int r2 = 0;
    for(int j=0; j<deg; j++){
        number a = n_Copy((number)(rootcont->getRoot(j)),comp);
        if( IsReal(a,comp)){
            roots[r1] = n_Copy(a,comp);
            r1++;
        }else if(ImagGreaterZero(a, comp)){
            complexroots[r2] = n_Copy(a,comp);
            r2++;
        }
        n_Delete(&a,comp);
    }
    rChangeCurrRing(CurrOldRing);
    DEBUG_PRINT(("delete some variables\n"));
    rDelete(newring);
    delete &n;
    //delete rootcont;DEBUG_PRINT(("can't delete rootcontainer rootcont\n"));
    DEBUG_VAR(r1);
    DEBUG_VAR(r2);
    for(int j=0;j<r2;j++){
        roots[r1+j]= n_Copy(complexroots[j],comp);
        n_Delete(&complexroots[j],comp);
    }
    delete complexroots;
    DEBUG_PRINT(("map elementarray to complex\n"));
    bigintmat ** elements = new bigintmat*[size_elementarray];
    for(int i=0;i<size_elementarray;i++){
        elements[i] = bimChangeCoeff(elementarray[i],comp);
    }
    DEBUG_PRINT(("generate output matrix\n"));
    DEBUG_PRINT(("real part\n"));
    bigintmat * complexmat = new bigintmat(r1+2*r2,size_elementarray,comp);
    for(int i=1; i<= r1; i++){
        number pot = n_Init(1,comp);
        for(int l=0; l< deg; l++){
            for(int j=0; j<size_elementarray;j++){
                number mult = n_Mult(pot,elements[j]->view(l),comp);
                complexmat->rawset(i,j+1,n_Add(complexmat->view(i,j+1),mult,comp),comp);
                n_Delete(&mult,comp);
            }
            n_InpMult(pot, roots[i-1],comp);
        }
        n_Delete(&pot,comp);
    }
    DEBUG_PRINT(("imaginary part\n"));
    number sqrt2 = n_Init(1,comp);
    if(r2>0){
        number two = n_Init(2,comp);
        number sqrt2 = squareroot(two,comp,precision+10);
        n_Delete(&two,comp);
        for(int i=1; i<= r2; i++){
            number pot = n_Init(1,comp);
            for(int l=0; l< deg; l++){
                for(int j=0; j<size_elementarray;j++){
                    number mult = n_Mult(pot,elements[j]->view(l),comp);
                    complexmat->rawset(r1+2*i,j+1,n_Add(complexmat->view(r1+2*i,j+1),mult,comp),comp);
                    n_Delete(&mult,comp);
                }
                n_InpMult(pot, roots[i-1],comp);
            }
            n_Delete(&pot,comp);
            for(int j=1;j<=size_elementarray;j++){
                complexmat->set(r1+2*i,j,n_Mult(complexmat->view(r1+2*i,j),sqrt2,comp),comp);
                complexmat->set(r1+2*i-1,j,n_RePart(complexmat->view(r1+2*i,j),comp),comp);
                complexmat->set(r1+2*i,j,n_ImPart(complexmat->view(r1+2*i,j),comp),comp);
            }
        }
    }
    DEBUG_PRINT(("delete Variables\n"));
    for(int i=0;i<size_elementarray;i++){
        delete elements[i];
    }
    delete elements;
    for(int i=0;i<r1+r2;i++){
        n_Delete(&roots[i],comp);
    }
    delete roots;
 
    DEBUG_PRINT(("to real\n"));
    LongComplexInfo paramReal;
    paramReal.float_len = si_min (precision, 32767);
    paramReal.float_len2 = si_min (precision, 32767);
    paramComp.par_name=(const char*)"i";
    coeffs real = nInitChar(n_long_R, &paramReal);
    //setGMPFloatDigits( precision, precision);
    bigintmat * realmat = bimChangeCoeff(complexmat,real);
    delete complexmat;
    return realmat;
}

scalarproduct()

number scalarproduct ( bigintmat *  a, bigintmat *  b  )

inline

Definition at line 915 of file lattice.cc.

                                                   {
    if(a->cols()!=1) {
        Werror("a->cols()!=1 in scalarproduct(a,b)\n");
        return NULL;
    }
    if(b->cols()!=1) {
        Werror("b->cols()!=1 in scalarproduct(a,b)\n");
        return NULL;
    }
    if(a->rows()!=b->rows()) {
        Werror("b->cols()!=1 in scalarproduct(a,b)\n");
        return NULL;
    }
    if(a->basecoeffs()!=b->basecoeffs()) {
        Werror("a->basecoeffs()!=b->basecoeffs() in scalarproduct(a,b)\n");
        return NULL;
    }
 
    coeffs coef = a->basecoeffs();
    number p = n_Init(0,coef);
    for(int i = 1; i <= b->rows(); i++){
        number prod = n_Mult(a->view(i,1), b->view(i,1), coef);
        n_InpAdd(p, prod, coef);
        n_Delete(&prod,coef);
    }
    return p;
}

squareroot()

number squareroot	(	number	a,
		coeffs	coef,
		int	prec
	)

Definition at line 1111 of file lattice.cc.

                                                  {
    if(n_IsZero(a,coef)){
        return a;
    }
    if(!n_GreaterZero(a,coef)){
        return n_Init(0,coef);
    }
    number two = n_Init(2,coef);
    number xn = n_Copy(a,coef);
    number xn1,xn2;
    for(int i=0;i<prec;i++){
        xn1 = n_Div(a,xn,coef);
        xn2 = n_Add(xn,xn1,coef);
        n_Delete(&xn,coef);
        xn = n_Div(xn2,two,coef);
        n_Delete(&xn1,coef);
        n_Delete(&xn2,coef);
    }
    return xn;
}

Variable Documentation

debug_block

STATIC_VAR bool debug_block = false

Definition at line 31 of file lattice.cc.