[mdk-re] использование gsl и gmp

[mdk-re] использование gsl и gmp

[Pavel Marakhovsky]

Привет!

Кто-нибудь gsl и gmp использовал? Интересует скорость расчетов для одномерной и многомерной
минимизации и численного интегрирования.

Возможен ли симбиоз этих двух библиотек, т.е типы данных из gmp, а функции из gsl?
А то gmp только предоставляет основные функции, а нужна минимизация (одномерная) и интегрирование,
ну или хотя бы интегрирование.

И еще один вопросик:) Пытался сделать документацию из через texi2dvi->dvips. PS получился,
только почему-то на 242 странице gv выдает вот это:

Error: /rangecheck in --get--
Operand stack:
   --nostringval--   --nostringval--   --nostringval--   descender   0   --nostringval--   1
Execution stack:
   %interp_exit   .runexec2   --nostringval--   --nostringval--   --nostringval--   2   %stopped_push
--nostringval--   --nostringval--   --nostringval--   false   1   %stopped_push   1   3   %oparray_pop   .runexec2
--nostringval--   --nostringval--   --nostringval--   2   %stopped_push   --nostringval--   --nostringval--
--nostringval--   --nostringval--
Dictionary stack:
   --dict:906/941(G)--   --dict:0/20(G)--   --dict:52/200(L)--   --dict:138/300(L)--   --dict:44/200(L)--
   --dict:39/50(L)--   --dict:2/17(L)--   --dict:5/17(L)--   --dict:1/3(L)--   --dict:14/17(L)--
Current allocation mode is local
Current file position is 594446
GNU Ghostscript: Unrecoverable error, exit code 1

Это глюк чего? ps или gv?

--
WBW, Unatine
e-mail: unatine@narod.ru
JID: unatine@jabber.org

Now playing: Horie Yui - Happy happy rice shower

[mdk-re] Re: [mdk-re] использование gsl и gmp

[Michael Shigorin]

On Fri, Sep 14, 2001 at 01:11:01PM +0400, Pavel Marakhovsky wrote:
> Кто-нибудь gsl и gmp использовал? Интересует скорость расчетов для одномерной и многомерной
> минимизации и численного интегрирования.
gmp использовали 2.0.2 и 3.0.1 и чрезвычайно довольны. Чрезвычайно.

> Возможен ли симбиоз этих двух библиотек, т.е типы данных из gmp, а функции из gsl?
> А то gmp только предоставляет основные функции, а нужна минимизация (одномерная) и интегрирование,
> ну или хотя бы интегрирование.
Вот на gsl посмотреть не успел, но я бы копал в сторону "возможности
сборки gsl с [поддержкой] gmp" -- чего тут еще придумаешь?

-- 
WBR, Michael Shigorin, webmaster of www.chem.univ.kiev.ua
>Home Page:  http://visa.chem.univ.kiev.ua/~mike/  ICQ: 113344029
>Brainbench: http://www.brainbench.com/transcript.jsp?pid=2434729

[mdk-re] Re: [mdk-re] Re: [mdk-re] использование gsl и gmp

[Sergey Krivulja]

Hello Michael,

Monday, September 17, 2001, 2:33:03 PM, you wrote:

MS> gmp использовали 2.0.2 и 3.0.1 и чрезвычайно довольны. Чрезвычайно.
gsl знаю, а gmp незнаю. Почему?
Ссылочку киньте плиз :)
-- 
Best regards,
 Sergey                            mailto:skrivulja@erec.ru

[mdk-re] Re: [mdk-re] Re: [mdk-re] Re: [mdk-re] использование gsl и gmp

[Michael Shigorin]

On Mon, Sep 17, 2001 at 05:28:31PM +0400, Sergey Krivulja wrote:
> MS> gmp использовали 2.0.2 и 3.0.1 и чрезвычайно довольны. Чрезвычайно.
> gsl знаю, а gmp незнаю. Почему?
Мне по науке нужно было здоровенные цифири без потери точности ворочать --
см. http://home.chem.univ.kiev.ua/~mike/works/BSc.ps.gz -- а оно это умеет
и _отменно_. Вплоть до asm на надцати архитектурах и CPU.

> Ссылочку киньте плиз :)
ftp.altlinux.ru 8-)

-- 
WBR, Michael Shigorin, webmaster of www.chem.univ.kiev.ua
>Home Page:  http://visa.chem.univ.kiev.ua/~mike/  ICQ: 113344029
>Brainbench: http://www.brainbench.com/transcript.jsp?pid=2434729

Re: [mdk-re] Re: [mdk-re] Re: [mdk-re] Re: [mdk-re] использование gsl и gmp

[Pavel Marakhovsky]

On Mon, 17 Sep 2001 17:21:28 +0300
Michael Shigorin <mike@lic145.kiev.ua> wrote:

> On Mon, Sep 17, 2001 at 05:28:31PM +0400, Sergey Krivulja wrote:
> > MS> gmp использовали 2.0.2 и 3.0.1 и чрезвычайно довольны. Чрезвычайно.
> > gsl знаю, а gmp незнаю. Почему?
> Мне по науке нужно было здоровенные цифири без потери точности ворочать --
> см. http://home.chem.univ.kiev.ua/~mike/works/BSc.ps.gz -- а оно это умеет
> и _отменно_. Вплоть до asm на надцати архитектурах и CPU.
Интересно:) А что за библиотечка saml? Где ее найти можно? а то freshmeat нет.


--
WBW, Unatine
e-mail: unatine@narod.ru
JID: unatine@jabber.org

Now playing: Horie Yui - Happy happy rice shower

[mdk-re] Re: [mdk-re] использование gsl и gmp

[Michael Shigorin]

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

On Tue, Sep 18, 2001 at 06:36:02AM +0400, Pavel Marakhovsky wrote:
> > > MS> gmp использовали 2.0.2 и 3.0.1 и чрезвычайно довольны. Чрезвычайно.
> > > gsl знаю, а gmp незнаю. Почему?
> > Мне по науке нужно было здоровенные цифири без потери точности ворочать --
> > см. http://home.chem.univ.kiev.ua/~mike/works/BSc.ps.gz -- а оно это умеет
> > и _отменно_. Вплоть до asm на надцати архитектурах и CPU.
> Интересно:) А что за библиотечка saml? Где ее найти можно? а то freshmeat нет.
LSM в аттаче. А в доке -- 
   The Simple Algebraic Math Library (SAML) is a low-level library to 
   handle the classic objects of symbolic calculus: arbitrary big     
   integers, reals of arbitrary precision, rationals, polynomials,
   tensors, et cetera. This document describes the internals of the
   library and other implementation details.

Оно умеет, например, так (за что и пользовали):
    v1 = mref_new(); v2 = mref_new();       
    mref_build(v2, ST_INTEGER, "0");
    mref_cast(v2, ST_APOLY);        
    /* Now v2 contains the null polynomial */
//                              ^^^^^^^^^^
    saml_init_lexer_mem("3*xy+1", 6);        
    saml_parse(v1, v2);
    /* Now v1 contains 3*xy+1, v2 is unchanged */
    saml_init_lexer_mem("x-y", 3);
    saml_parse(v2, v2);           

Только я обленился писать на сях вокруг, взял имеющийся в комплекте
интерпретатор и ваял на нем.

Насколько я понимаю, проект все [saml-970418] :(

И смутно припоминаю, что нашел что-то ну никак не хуже, но использующее
libgmp (SAML -- нет).

А, кстати -- имеет смысл обратить внимание на yacas:

Name        : yacas                        Relocations: (not relocateable)
Version     : 1.0.45                            Vendor: (none)
Release     : 1                             Build Date: Mon Apr 23 17:00:10 2001
Install date: Mon Apr 23 18:33:16 2001      Build Host: work.fair.net
Group       : Applications/Math             Source RPM: yacas-1.0.45-1.src.rpm
Size        : 1161000                          License: LGPL
URL         : http://www.xs4all.nl/~apinkus/yacas.html
Summary     : Yet Another Computer Algebra System
Description :
Yacas (Yet Another Computer Algebra System) is a small and highly flexible
computer algebra language. The syntax uses a infix-operator grammar
parser. The distribution contains a small library of mathematical
functions, but its real strength is in the language in which you can
easily write your own symbolic manipulation algorithms. It supports
arbitrary precision arithmetic.

-- 
WBR, Michael Shigorin, webmaster of www.chem.univ.kiev.ua
>Home Page:  http://visa.chem.univ.kiev.ua/~mike/  ICQ: 113344029
>Brainbench: http://www.brainbench.com/transcript.jsp?pid=2434729

[-- Attachment #2: saml.lsm --]
[-- Type: text/plain, Size: 1189 bytes --]

Begin3
Title:		SAML, the Simple Algebraic Math Library
Version:	970418
Entered-date:	18APR97
Description:	SAML is a C library for symbolic calculations, accompanied
		by some application programs (samuel, factorint, induce),
		and Python bindings.
		The library provides an object-oriented framework for
		defining and handling mathematical types, and implements
		the most common data types of computer algebra: integers,
		reals, fractions, complex numbers, polynomials, tensors,
		matrices, etc.
		The application programs consist of an interactive symbolic
		calculator (samuel), a programming language (induce) and
		a program to factorize integers (factorint).
Keywords:	Computer algebra, C, Python, library, symbolic calculus,
		multiple-precision arithmetic, quadratic sieve
Author:		bousch@topo.math.u-psud.fr (Thierry Bousch)
Maintained-by:	bousch@topo.math.u-psud.fr (Thierry Bousch)
Primary-site:	topo.math.u-psud.fr /pub/bousch
		243kB saml-970418.tgz
Platforms:	Requires a decent Un*x system, gcc 2.5.0 or later, gdbm,
		perl, gmake, and TeXinfo for the documentation.
		It is strongly recommended to have Python and the
		readline library installed.
Copying-policy:	GPL
End

Re: [mdk-re] Re: [mdk-re] Re: [mdk-re] использование gsl и gmp

[Pavel Marakhovsky]

On Mon, 17 Sep 2001 17:28:31 +0400
Sergey Krivulja <skrivulja@erec.ru> wrote:

> MS> gmp использовали 2.0.2 и 3.0.1 и чрезвычайно довольны. Чрезвычайно.
> gsl знаю, а gmp незнаю. Почему?
> Ссылочку киньте плиз :)

Есть и в Спринге и в Сизифе: libgmp
Вот урла: http://www.swox.com/gmp/
--
WBW, Unatine
e-mail: unatine@narod.ru
JID: unatine@jabber.org

Now playing: Sakamoto Maya - Aoi hitomi

Re: [mdk-re] Re: [mdk-re] использование gsl и gmp

[Pavel Marakhovsky]

On Mon, 17 Sep 2001 13:33:03 +0300
Michael Shigorin <mike@lic145.kiev.ua> wrote:

> > Кто-нибудь gsl и gmp использовал? Интересует скорость расчетов для одномерной и многомерной
> > минимизации и численного интегрирования.
> gmp использовали 2.0.2 и 3.0.1 и чрезвычайно довольны. Чрезвычайно.
Ага.... начал переписывать программу... потом обнаружилось что для mpf_t не функции экспоненты.
И вообще такой нет. А mpf_pow_ui только в целую степень:(((
Мелькнула, конечно, мысля написать функцию для приближенного вычисления exp, но как подумал
сколько тогда эта прога работать будет, так мысля сразу и исчезла:( Блин... ну такую функцию и не сделать...

> > Возможен ли симбиоз этих двух библиотек, т.е типы данных из gmp, а функции из gsl?
> > А то gmp только предоставляет основные функции, а нужна минимизация (одномерная) и интегрирование,
> > ну или хотя бы интегрирование.
> Вот на gsl посмотреть не успел, но я бы копал в сторону "возможности
> сборки gsl с [поддержкой] gmp" -- чего тут еще придумаешь?
Не могет оно такое. Только стандартные параметры для configure (даже внутри смотрел).

Так что похоже придется одно gsl использовать, может получится выигрыш в скорости работы программ
 хотябы.
--
WBW, Unatine
e-mail: unatine@narod.ru
JID: unatine@jabber.org

Now playing: Sakamoto Maya - Aoi hitomi