From owner-FreeBSD-users-jp@jp.freebsd.org  Thu Oct  2 23:29:56 1997
Received: by jaz.jp.freebsd.org (8.8.7+2.7Wbeta7/8.7.3) id XAA25612
	Thu, 2 Oct 1997 23:29:56 +0900 (JST)
Received: by jaz.jp.freebsd.org (8.8.7+2.7Wbeta7/8.7.3) with ESMTP id XAA25605
	for <FreeBSD-users-jp@jp.freebsd.org>; Thu, 2 Oct 1997 23:29:54 +0900 (JST)
Received: from singer.kt.rim.or.jp (ppp369.kt.rim.or.jp [202.247.140.69]) by mail.kt.rim.or.jp (8.8.5/3.4W3-rim1.1) with ESMTP id XAA08443; Thu, 2 Oct 1997 23:29:33 +0900 (JST)
Received: from localhost (localhost [127.0.0.1]) by singer.kt.rim.or.jp (8.8.5/3.4Wbeta3-96070610) with ESMTP id XAA00486; Thu, 2 Oct 1997 23:20:40 +0900 (JST)
To: FreeBSD-users-jp@jp.freebsd.org, ttakamat@nr.titech.ac.jp
In-Reply-To: Your message of "Thu, 02 Oct 1997 19:50:27 +0900"
	<19971002195027R.ttakamat@nr.titech.ac.jp>
References: <19971002195027R.ttakamat@nr.titech.ac.jp>
X-Mailer: Mew version 1.91 on Emacs 19.34 / Mule 2.3 (SUETSUMUHANA)
Mime-Version: 1.0
Content-Type: Text/Plain; charset=iso-2022-jp
Content-Transfer-Encoding: 7bit
Message-Id: <19971002232039H.yoshiaki@_domain_of_your_mail_address_>
Date: Thu, 02 Oct 1997 23:20:39 +0900
From: Yoshiaki Uchikawa <yoshiaki@_domain_of_your_mail_address_.rim.or.jp>
X-Dispatcher: imput version 970918
Lines: 84
Reply-To: FreeBSD-users-jp@jp.freebsd.org
Precedence: bulk
X-Distribute: distribute [version 2.1 (Alpha) patchlevel=20]
X-Sequence: FreeBSD-users-jp 18854
Subject: [FreeBSD-users-jp 18854] Re: =?iso-2022-jp?B?GyRCO01HXEA6RVkbKEI=?=(GCC)
Errors-To: owner-FreeBSD-users-jp@jp.freebsd.org
Sender: owner-FreeBSD-users-jp@jp.freebsd.org

$B$&$A$+$o(B $B!s(B $B!V%l%9!W$H$+!V%$%s%?%M!W$H$+Lu$N$o$+$i$J$$8@MU$r;H$C$FDL$V$k(B
            $B$N$O$_$C$H$b$J$$$H;W$&(B:-) $B$G$9(B

$B!t(B mew $B$N%F%9%H$b7s$M$F$$$^$9(B

> $B<B$O!"(Bgcc$B$K$D$$$F$*J9$-$7$?$$$3$H$,$"$j$^$9!#(B
> $B8=:_!"8&5f$G?tCM7W;;$r9T$C$F$$$k$N$G$9$,!"5^$-$g;MG\@:EY$G9T$J$&I,MW@-$,(B
> $B$G$-$^$7$?!#$=$3$G!"(Bgcc$B$G$G$-$J$$$b$N$+$HD4$Y$?=j(B(limits.h$B$J$I(B)$B$G$-$J$$(B
> ($B;d$,8+$?46$8$G$O(B)$B$h$&$G$7$?!#(B
> $B$^$?!"0J2<$NK\$r8+$^$7$?$,!";MB'1i;;$N%k!<%A%s$r:n$k$3$H$GB?G\D91i;;$r(B
> $B=hM}$7$F$*$j$^$7$?!#(B

$B9b>>$5$s$N%$%a!<%8$H$O>/$7(B($B$+$J$j(B?)$B0c$&$+$bCN$l$^$;$s$,!"(B4$BG\$J$I$H$$$&(B
$B$1$A$J$b$N$G$O$J$/I8=`@:EY(B 28 digit $B$G(B REAL$B$J$i(B 45 digit$B$^$G(B $B!"$b$A$m(B
$B$sJ#AG?t$b%5%]!<%H!"?t3X$N@lLg2H0J30$K$OJ9$$$?$3$H$b$J$$$h$&$J7?$^$G%5(B
$B%]!<%H$7$F$$$F!"1i;;B.EY$b$^$"$^$"B.$$(B pari $B$H$$$&D66/NO$J%i%$%V%i%j$,(B
$B$"$j$^$9!#(B

$BAH$_9~$_4X?t$bHs>o$K<oN`$,B?$/;d$J$I$O2?$K;H$&$+$o$+$i$J$$$b$N$NJ}$,B?(B
$B$$$/$i$$$G$9(B(^_^;

$B@NI,MW$KGw$i$l$F!";H$C$?$3$H$b$"$C$F(B ports $B$K$7$F$"$j$^$9!#(B 
ports-current $B$K$"$j$^$9$N$G$b$7$h$+$C$?$i8+$F$/$@$5$$!#(B /math/pari $B$G(B
$B$9!#(B

$B!t$a$s$I$&$@$C$?$N$G6&M-%i%$%V%i%j$K$7$F$^$;$s!#$3$N$?$a$3$N%i%$%V%i%j(B
$B!t$r%j%s%/$9$k$H%P%$%J%j$,!V$H$F$b!WBg$-$/$J$j$^$9(B

$B$^$"!"$3$s$J46$8$G$9(B
-------
  static GEN gmat;
  /* n:$B%^%H%j%/%9$NBg$-$5(B */
  pari_init(n)
  int n;
  {
     int i;

    init(16000000,2);
    gmat = cgetg(n+1,19);
    for(i = 1; i <= n; i++)
      gmat[i] = lgetg(n+1,17);
  }

  rgauss(ab,eps,mlen,da,coef)
  double ab[MATRIX3][MATRIX3],eps,*da,*coef;
  int mlen;
  {
    int k,l;
    long ltop,lbtm;
    GEN vec,tmp;

    vec = cgetg(mlen+1,17);
    for (k = 1; k <= mlen; k++)
      vec[k] = (GENC)cgetg(MSIZE,2);

    for(k = 0; k < mlen; k++){
      ltop = avma;
      for(l = 0; l < mlen; l++){
        tmp = gclone(dbltor(ab[l][k]));
        setlg(tmp,MSIZE);
        lbtm = avma;
        vec[l+1] = (GENC)gclone(tmp);
      }
      /* $B%k!<%WCf$G>CHq$5$l$?JQ?tNN0h$N2s<}(B */
      vec = gerepile(ltop,lbtm,vec); 
      gmat[k+1] = (GENC)gcopy(vec);
    }
------

$B;H$&>e$G$N:GBg$N7gE@$O!"JQ?tNN0h$N4IM}$r<+J,$G$d$i$J$1$l$P$J$i$J$$$H$$(B
$B$&$3$H$G$9!#$3$N$?$a$K$A$g$C$HL}CG$9$k$H4JC1$KJQ?tNN0h$r;H$$2L$?$7$F$7(B
$B$^$$$^$9!#$=$&$$$&0UL#$G$O$+$J$j;H$&$N$,Fq$7$$$H$$$($^$9!#(B

$B!tMW$9$k$K%,%Y!<%8%3%l%/%7%g%s$rA4It<+J,$G$d$i$J$/$F$O$J$i$J$$(B Lisp $B$r(B
$B!t;H$&$h$&$J$b$N$G$9!#(B-- $B$-$C$HFbIt=hM}$O(B Lisp $BIw$J$s$@$m$&(B :-)

$B$?$@!"%I%-%e%a%s%H$,$+$J$j$7$C$+$j$7$F$$$k$3$H!"$=$l$J$j$KMxMQ2ACM$N9b(B
$B$$%5%s%W%k%"%W%j%1!<%7%g%s$,IUB0$7$F$$$k$3$H$,$"$j$^$9$N$G!"$8$C$/$j;H$C(B
$B$F$_$k2ACM$,$"$k$+$b$7$l$^$;$s!#(B
----
  yoshiaki@kt.rim.or.jp
 ($B$&(B)


