padic.h

/*=============================================================================

    This file is part of FLINT.

    FLINT is free software; you can redistribute it and/or modify
    it under the terms of the GNU General Public License as published by
    the Free Software Foundation; either version 2 of the License, or
    (at your option) any later version.

    FLINT is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU General Public License for more details.

    You should have received a copy of the GNU General Public License
    along with FLINT; if not, write to the Free Software
    Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA  02110-1301 USA

=============================================================================*/
/******************************************************************************

    Copyright (C) 2011 Sebastian Pancratz
 
******************************************************************************/

#ifndef PADIC_H
#define PADIC_H

#undef ulong /* interferes with system includes */
#include <stdlib.h>
#include <stdio.h>
#define ulong unsigned long

#include <mpir.h>

#include "flint.h"
#include "fmpz.h"
#include "fmpq.h"
#include "fmpz_vec.h"
#include "ulong_extras.h"

typedef struct {
    fmpz u;
    long v;
} padic_struct;

typedef padic_struct padic_t[1];

#define padic_unit(x)  (&((x)->u))
#define padic_val(x)   ((x)->v)

enum padic_print_mode
{
    PADIC_TERSE, 
    PADIC_SERIES, 
    PADIC_VAL_UNIT
};

typedef struct {

    fmpz_t p;
    long N;

    double pinv;

    fmpz *pow;
    long min;
    long max;

    enum padic_print_mode mode;

} padic_ctx_struct;

typedef padic_ctx_struct padic_ctx_t[1];

typedef struct {
    long n;
    fmpz *pow;
    fmpz *u;
} padic_inv_struct;

typedef padic_inv_struct padic_inv_t[1];

/* Context *******************************************************************/

void padic_ctx_init(padic_ctx_t ctx, const fmpz_t p, long N, 
                    enum padic_print_mode mode);

void padic_ctx_clear(padic_ctx_t ctx);

static __inline__ int 
_padic_ctx_pow_ui(fmpz_t rop, ulong e, const padic_ctx_t ctx)
{
    if (ctx->min <= e && e < ctx->max)
    {
        *rop   = *(ctx->pow + (e - ctx->min));
        return 0;
    }
    else
    {
        fmpz_init(rop);
        fmpz_pow_ui(rop, ctx->p, e);
        return 1;
    }
}

/* Memory management *********************************************************/

void _padic_init(padic_t rop);

void padic_init(padic_t rop, const padic_ctx_t ctx);

void _padic_clear(padic_t rop);

void padic_clear(padic_t rop, const padic_ctx_t ctx);

static __inline__ void _padic_canonicalise(padic_t rop, const padic_ctx_t ctx)
{
    if (!fmpz_is_zero(padic_unit(rop)))
    {
        padic_val(rop) += _fmpz_remove(padic_unit(rop), ctx->p, ctx->pinv);
    }
    else
    {
        padic_val(rop) = 0;
    }
}

void _padic_reduce(padic_t rop, const padic_ctx_t ctx);

void padic_reduce(padic_t rop, const padic_ctx_t ctx);

/* Randomisation *************************************************************/

void padic_randtest(padic_t rop, flint_rand_t state, const padic_ctx_t ctx);

void padic_randtest_not_zero(padic_t rop, flint_rand_t state, 
                             const padic_ctx_t ctx);

/* Assignments and conversions ***********************************************/

void _padic_set(padic_t rop, const padic_t op);
void padic_set(padic_t rop, const padic_t op, const padic_ctx_t ctx);

void _padic_set_si(padic_t rop, long op, const padic_ctx_t ctx);
void padic_set_si(padic_t rop, long op, const padic_ctx_t ctx);

void _padic_set_ui(padic_t rop, ulong op, const padic_ctx_t ctx);
void padic_set_ui(padic_t rop, ulong op, const padic_ctx_t ctx);

void _padic_set_fmpz(padic_t rop, const fmpz_t op, const padic_ctx_t ctx);
void padic_set_fmpz(padic_t rop, const fmpz_t op, const padic_ctx_t ctx);

void padic_set_fmpq(padic_t rop, const fmpq_t op, const padic_ctx_t ctx);

void _padic_set_mpz(padic_t rop, const mpz_t op, const padic_ctx_t ctx);
void padic_set_mpz(padic_t rop, const mpz_t op, const padic_ctx_t ctx);

void padic_set_mpq(padic_t rop, const mpq_t op, const padic_ctx_t ctx);

void _padic_get_fmpz(fmpz_t rop, const padic_t op, const padic_ctx_t ctx);
void padic_get_fmpz(fmpz_t rop, const padic_t op, const padic_ctx_t ctx);

void _padic_get_fmpq(fmpq_t rop, const padic_t op, const padic_ctx_t ctx);
void padic_get_fmpq(fmpq_t rop, const padic_t op, const padic_ctx_t ctx);

void _padic_get_mpz(mpz_t rop, const padic_t op, const padic_ctx_t ctx);
void padic_get_mpz(mpz_t rop, const padic_t op, const padic_ctx_t ctx);

void _padic_get_mpq(mpq_t rop, const padic_t op, const padic_ctx_t ctx);
void padic_get_mpq(mpq_t rop, const padic_t op, const padic_ctx_t ctx);

static __inline__ void padic_swap(padic_t op1, padic_t op2)
{
    long t;

    fmpz_swap(padic_unit(op1), padic_unit(op2));
    t              = padic_val(op1);
    padic_val(op1) = padic_val(op2);
    padic_val(op2) = t;
}

static __inline__ void padic_zero(padic_t rop)
{
    fmpz_zero(padic_unit(rop));
    padic_val(rop) = 0;
}

static __inline__ void _padic_one(padic_t rop)
{
    fmpz_one(padic_unit(rop));
    padic_val(rop) = 0;
}

static __inline__ void padic_one(padic_t rop, const padic_ctx_t ctx)
{
    if (ctx->N > 0)
        _padic_one(rop);
    else
        padic_zero(rop);
}

/* Arithmetic operations *****************************************************/

void _padic_add(padic_t rop, const padic_t op1, const padic_t op2, 
                const padic_ctx_t ctx);

void padic_add(padic_t rop, const padic_t op1, const padic_t op2, 
               const padic_ctx_t ctx);

void _padic_sub(padic_t rop, const padic_t op1, const padic_t op2, 
                const padic_ctx_t ctx);

void padic_sub(padic_t rop, const padic_t op1, const padic_t op2, 
               const padic_ctx_t ctx);

void _padic_neg(padic_t rop, const padic_t op);

void padic_neg(padic_t rop, const padic_t op, const padic_ctx_t ctx);

void _padic_mul(padic_t rop, const padic_t op1, const padic_t op2);

void padic_mul(padic_t rop, const padic_t op1, const padic_t op2, 
               const padic_ctx_t ctx);

void padic_shift(padic_t rop, const padic_t op, long v, const padic_ctx_t ctx);

void padic_div(padic_t rop, const padic_t op1, const padic_t op2, 
               const padic_ctx_t ctx);

void _padic_inv_precompute(padic_inv_t S, const fmpz_t p, long N);

void _padic_inv_clear(padic_inv_t S);

void _padic_inv_precomp(fmpz_t rop, const fmpz_t op, padic_inv_t S);

void _padic_inv(fmpz_t rop, const fmpz_t op, const fmpz_t p, long N);

void padic_inv(padic_t rop, const padic_t op, const padic_ctx_t ctx);

int padic_sqrt(padic_t rop, const padic_t op, const padic_ctx_t ctx);

void padic_pow_si(padic_t rop, const padic_t op, long e, 
                  const padic_ctx_t ctx);

/* Comparison ****************************************************************/

static __inline__ int _padic_is_zero(const padic_t op)
{
    return fmpz_is_zero(padic_unit(op));
}

static __inline__ int padic_is_zero(const padic_t op, const padic_ctx_t ctx)
{
    return fmpz_is_zero(padic_unit(op)) || (padic_val(op) >= ctx->N);
}

static __inline__ int _padic_is_one(const padic_t op)
{
    return fmpz_is_one(padic_unit(op)) && (padic_val(op) == 0);
}

static __inline__ int padic_is_one(const padic_t op, const padic_ctx_t ctx)
{
    if (ctx->N > 0)
        return _padic_is_one(op);
    else
        return (padic_val(op) >= ctx->N);
}

int _padic_equal(const padic_t op1, const padic_t op2);

int padic_equal(const padic_t op1, const padic_t op2, const padic_ctx_t ctx);

/* Special functions *********************************************************/

void padic_teichmuller(padic_t rop, const padic_t op, const padic_ctx_t ctx);

int padic_exp(padic_t rop, const padic_t op, const padic_ctx_t ctx);

int padic_log(padic_t rop, const padic_t op, const padic_ctx_t ctx);

ulong padic_val_fac_ui2(ulong N);

ulong padic_val_fac_ui(ulong N, const fmpz_t p);

void padic_val_fac(fmpz_t rop, const fmpz_t op, const fmpz_t p);

/* Input and output **********************************************************/

char * _padic_get_str(char * str, const padic_t op, const padic_ctx_t ctx);

char * padic_get_str(char * str, const padic_t op, const padic_ctx_t ctx);

int __padic_fprint(FILE * file, const fmpz_t u, long v, const padic_ctx_t ctx);

int _padic_fprint(FILE * file, const padic_t op, const padic_ctx_t ctx);

int padic_fprint(FILE * file, const padic_t op, const padic_ctx_t ctx);

static __inline__ int 
__padic_print(const fmpz_t u, long v, const padic_ctx_t ctx)
{
    return __padic_fprint(stdout, u, v, ctx);
}

static __inline__ int _padic_print(const padic_t op, const padic_ctx_t ctx)
{
    return _padic_fprint(stdout, op, ctx);
}

static __inline__ int padic_print(const padic_t op, const padic_ctx_t ctx)
{
    return padic_fprint(stdout, op, ctx);
}

static __inline__ void _padic_debug(const padic_t op)
{
    printf("{u = ");
    fmpz_print(padic_unit(op)); 
    printf(", v = %ld}", padic_val(op));
}

#endif