caml-list - the Caml user's mailing list
 help / color / mirror / Atom feed
From: Christophe Raffalli <christophe.raffalli@univ-savoie.fr>
To: Thomas Gazagnaire <thomas.gazagnaire@irisa.fr>,
	caml-list <caml-list@inria.fr>
Subject: Re: [Caml-list] Integral solutions of rational linear equations
Date: Wed, 30 Nov 2005 09:49:38 +0100	[thread overview]
Message-ID: <438D67A2.4050500@univ-savoie.fr> (raw)
In-Reply-To: <438CA20C.2080200@irisa.fr>

[-- Attachment #1: Type: text/plain, Size: 1722 bytes --]

Thomas Gazagnaire a écrit :
> Hello,
> 
> in my one of my Ocaml programs, I need to find all integers solutions of
> a rational equation systems. This algo uses Gauss reduction and Hermite
> normal form, and need to know if a rational is an integer or not (ie. I
> don't want to use numerical approximation : (1/3) * 3 is an integer but
> 0,333333*3 we don't know). I didn't find any integer algebra package for
> ocaml, so I tried to implement Gauss elimination (easy) and Hermite
> normal form (more difficult...). But I didn't implement optimized
> version of these algorithms...
> 
> So my question is : do you know if exists a native ocaml module or an
> interface with a C library which is able to do integer/rational matrix
> manipulation (essentialy the Hermite normal form) in an efficient way ?
> 

you did not look very hard ? OCaml is distributed with the bignum
library in standard which does both int and rationnal ...

$ ocaml nums.cma
        Objective Caml version 3.09.0

# open Num

let rec fact n = if n = 0 then Int 1 else fact (n-1) */ Int n
      ;;
val fact : int -> Num.num = <fun>
# fact 100;;
- : Num.num = Big_int <abstr>
# string_of_num (fact 100);;
- : string =
"93326215443944152681699238856266700490715968264381621468592963895217599993229915608941463976156518286253697920827223758251185210916864000000000000000000000000"

> Thomas
> 
> 
> _______________________________________________
> Caml-list mailing list. Subscription management:
> http://yquem.inria.fr/cgi-bin/mailman/listinfo/caml-list
> Archives: http://caml.inria.fr
> Beginner's list: http://groups.yahoo.com/group/ocaml_beginners
> Bug reports: http://caml.inria.fr/bin/caml-bugs


[-- Attachment #2: OpenPGP digital signature --]
[-- Type: application/pgp-signature, Size: 894 bytes --]

  parent reply	other threads:[~2005-11-30  8:49 UTC|newest]

Thread overview: 17+ messages / expand[flat|nested]  mbox.gz  Atom feed  top
2005-09-02 12:02 Unsafe features Florian Weimer
2005-09-03  9:24 ` [Caml-list] " Damien Bobillot
2005-09-03  9:40   ` Florian Weimer
2005-09-03 11:19     ` Jon Harrop
2005-09-03 12:07     ` yoann padioleau
2005-09-03 14:28       ` Florian Weimer
2005-09-03 14:35         ` yoann padioleau
2005-09-03 14:47           ` Florian Weimer
2005-09-03 14:51             ` Florian Weimer
2005-09-03 14:55               ` yoann padioleau
2005-09-04  1:58         ` Jacques Garrigue
2005-09-03  9:29 ` Erik de Castro Lopo
2005-11-29 18:46 ` Integral solutions of rational linear equations Thomas Gazagnaire
2005-11-30  8:32   ` [Caml-list] " Sebastian Egner
2005-11-30  8:49   ` Christophe Raffalli [this message]
2005-11-30  9:06     ` Christophe Raffalli
2005-11-30  9:08     ` Christophe Raffalli

Reply instructions:

You may reply publicly to this message via plain-text email
using any one of the following methods:

* Save the following mbox file, import it into your mail client,
  and reply-to-all from there: mbox

  Avoid top-posting and favor interleaved quoting:
  https://en.wikipedia.org/wiki/Posting_style#Interleaved_style

* Reply using the --to, --cc, and --in-reply-to
  switches of git-send-email(1):

  git send-email \
    --in-reply-to=438D67A2.4050500@univ-savoie.fr \
    --to=christophe.raffalli@univ-savoie.fr \
    --cc=caml-list@inria.fr \
    --cc=thomas.gazagnaire@irisa.fr \
    /path/to/YOUR_REPLY

  https://kernel.org/pub/software/scm/git/docs/git-send-email.html

* If your mail client supports setting the In-Reply-To header
  via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for NNTP newsgroup(s).