Re: sparse matrix

Karim Belabas on Tue, 31 Oct 2006 09:01:01 +0100

[Date Prev] [Date Next] [Thread Prev] [Thread Next] [Date Index] [Thread Index]

Re: sparse matrix

To: pari-users list <pari-users@list.cr.yp.to>
Subject: Re: sparse matrix
From: "Karim Belabas" <Karim.Belabas@math.u-bordeaux1.fr>
Date: Tue, 31 Oct 2006 07:38:45 +0100
Delivery-date: Tue, 31 Oct 2006 09:01:01 +0100
In-reply-to: <Pine.LNX.4.64.0610292306060.18040@prism.math.mcgill.ca>
Mail-followup-to: pari-users list <pari-users@list.cr.yp.to>
Mailing-list: contact pari-users-help@list.cr.yp.to; run by ezmlm
References: <Pine.LNX.4.64.0610292306060.18040@prism.math.mcgill.ca>
User-agent: Mutt/1.5.9i

* Mak Trifkovic [2006-10-30 05:24]:
> I have a big (~600x60, very sparse matrix (~3 non-zero entries per row). 
> Is there a clever way to get it into Hilbert normal form in real time?

What is Hilbert normal form ? If you meant Hermite Normal Form, then
your best bet using only public interfaces is mathnfmod(M, matdetint(M)).

Much more efficient internal functions are available, but require
installing the undocumented function mathnfspec().

Cheers,

    K.B.
-- 
Karim Belabas                  Tel: (+33) (0)5 40 00 26 17
Universite Bordeaux 1          Fax: (+33) (0)5 40 00 69 50
351, cours de la Liberation    http://www.math.u-bordeaux.fr/~belabas/
F-33405 Talence (France)       http://pari.math.u-bordeaux.fr/  [PARI/GP]

References:
- sparse matrix
  - From: Mak Trifkovic <mak@math.mcgill.ca>

Prev by Date: Re: confluent hypergeometric function
Previous by thread: sparse matrix
Index(es):
- Date
- Thread