From help-octave-request@bevo.che.wisc.edu Fri Jul 3 17:07:14 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id RAA14755 for help-octave-dist@bevo.che.wisc.edu; Fri, 3 Jul 1998 17:07:14 -0500 (CDT) Resent-Date: Fri, 3 Jul 1998 17:07:14 -0500 (CDT) Date: Fri, 3 Jul 1998 18:02:56 -0400 (EDT) From: Andy Adler X-Sender: aadler@pluto.dewbio.ca Reply-To: Andy Adler To: help-octave@bevo.che.wisc.edu Subject: accessing structures from oct files Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"yBTqq8_XFYK.A.I2.SYVn1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Hi, I'm trying to write *.oct files which access structures (ideally to hack a way of handling sparse matrices into octave) The idea being that a sparse matrix S would be represented S.nrow -> number of rows S.ncol -> number of cols S.nnz -> number of nonzero elements S.nzval -> vector of nonzero elements S.rowind-> vector index to rows S.colptr-> vector pointer to rows S.val -> values I could then call the superlu package (which already has mex file examples for matlab) However, I'm wondering how to access these structures from within the *oct file. I'd like to get octave_value A_arg = args(0); Matrix Aval= A_arg.val.matrix_value(); ^^^ But (of course) this doesn't work. Thanks for any help, _____________________________________________________________________ Andy Adler | American Biometric Company | Tel:(613)736-5100x154 aadler@abio.com | 3429 Hawthorne Rd,Ottawa,ON | Fax:(613)736-1348 From help-octave-request@bevo.che.wisc.edu Mon Jul 6 16:36:15 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id QAA00617 for help-octave-dist@bevo.che.wisc.edu; Mon, 6 Jul 1998 16:36:15 -0500 (CDT) Resent-Date: Mon, 6 Jul 1998 16:36:15 -0500 (CDT) Message-ID: <35A1423E.D1BCE80E@frb.gov> Date: Mon, 06 Jul 1998 17:31:42 -0400 From: "Steven W. Sumner" X-Mailer: Mozilla 4.05 [en] (X11; U; SunOS 5.4 sun4c) MIME-Version: 1.0 To: help-octave@bevo.che.wisc.edu Subject: Memory usage in octave Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Resent-Message-ID: <"BrBeKi8me0C.A.4rG.PNUo1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu I am trying to understand how octave's memory usage works. I have a simple program and function file which I have attached (see below) that seems to be leaking memory. The same program called in matlab leaks no memory. Could someone explain the logic of memory usage in this simple program. I have a much larger and more complex program with calls to .oct files and functions, which also leaks memory. The main program file (testmem.m) has the following lines of code: for i=1:500, x=rand(250); xy=testfun(x); clear end ___________________________________________________________ The function file (testfun.m) contains the following code: function xx = testfun(p) xx=p+p; xx(1:100,:)=[]; xx(:,1:100)=[]; return When this program is run in matlab it takes up 8488K (as measured by top running on a unix workstation Solaris 2.5.1). In Octave (2.0.13) it starts out using 9208K but grows by the 500th time through the loop to 9416K. What is causing this memory leak? Also, I thought it was odd that the clear command did not also delete the counter "i" and cause the program to bomb. Why does this not occur. Any insights into memory usage in octave would be appreciated. Thanks, Steve Sumner Federal Reserve Board From help-octave-request@bevo.che.wisc.edu Thu Jul 9 09:55:27 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id JAA04464 for help-octave-dist@bevo.che.wisc.edu; Thu, 9 Jul 1998 09:55:27 -0500 (CDT) Resent-Date: Thu, 9 Jul 1998 09:55:27 -0500 (CDT) Message-Id: <3.0.2.32.19980709165533.0080b460@fabi.vub.ac.be> X-Sender: borosy@fabi.vub.ac.be X-Mailer: QUALCOMM Windows Eudora Light Version 3.0.2 (32) Date: Thu, 09 Jul 1998 16:55:33 +0200 To: help-octave@bevo.che.wisc.edu From: =?iso-8859-1?Q?Andr=E1s?= Borosy Subject: path Mime-Version: 1.0 Content-Type: text/plain; charset="us-ascii" Resent-Message-ID: <"VwBaQuLtcKM.A.EEB.fnNp1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu I have sevaral Matlab .m files in a direcrory. How can I set Octave to access all these files, I think I should set something which is similar to matlabarc file. Many thank, Andras Borosy From help-octave-request@bevo.che.wisc.edu Tue Jul 14 11:17:52 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id LAA30662 for help-octave-dist@bevo.che.wisc.edu; Tue, 14 Jul 1998 11:17:52 -0500 (CDT) Resent-Date: Tue, 14 Jul 1998 11:17:52 -0500 (CDT) Message-ID: <35AB848B.5974C10F@hrp.no> Date: Tue, 14 Jul 1998 18:17:16 +0200 From: Terje Tverberg Organization: OECD Halden Reactor Project X-Mailer: Mozilla 4.05 [en] (X11; I; HP-UX B.10.20 9000/780) MIME-Version: 1.0 To: help-octave@bevo.che.wisc.edu Subject: subscribe Content-Type: multipart/alternative; boundary="------------EFB456E19E6B1DEB7ED9B2C9" Resent-Message-ID: <"xQ8aJ7_k7eC.A.ZeH.wS4q1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu --------------EFB456E19E6B1DEB7ED9B2C9 Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit I wish to subscribe to this mailing list. Regards Terje Tverberg -- ------------------------------------------------ Terje Tverberg Institutt for energiteknikk OECD Halden Reactor Project P.O. Box 175, N-1751 Halden, Norway Tel: +47 69 21 23 43 Fax: +47 69 21 22 01 E-mail: Terje.Tverberg@hrp.no http://www.ife.no/ ------------------------------------------------ --------------EFB456E19E6B1DEB7ED9B2C9 Content-Type: text/html; charset=us-ascii Content-Transfer-Encoding: 7bit I wish to subscribe to this mailing list.

Regards

Terje Tverberg

-- 

------------------------------------------------
 Terje Tverberg 
 Institutt for energiteknikk
 OECD Halden Reactor Project
 P.O. Box 175, N-1751 Halden, Norway
 Tel:    +47 69 21 23 43  
 Fax:    +47 69 21 22 01
 E-mail: Terje.Tverberg@hrp.no
 http://www.ife.no/
------------------------------------------------
  --------------EFB456E19E6B1DEB7ED9B2C9-- From help-octave-request@bevo.che.wisc.edu Tue Jul 14 13:32:18 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id NAA17313 for help-octave-dist@bevo.che.wisc.edu; Tue, 14 Jul 1998 13:32:18 -0500 (CDT) Resent-Date: Tue, 14 Jul 1998 13:32:18 -0500 (CDT) Message-Id: <199807141832.NAA25881@Eng.Auburn.EDU> X-Mailer: Microsoft Outlook Express for Macintosh - 4.0c (197) Date: Tue, 14 Jul 1998 07:54:01 -0500 Subject: sp package in octave From: "A. Scottedward Hodel" To: help-octave@bevo.che.wisc.edu Mime-version: 1.0 X-Priority: 3 Content-type: text/plain; charset="KOI8-R" Content-transfer-encoding: 7bit Resent-Message-ID: <"MR37FzQu-UI.A.t2H.yQ6q1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Has anyone ported Vandenberghe and Boyd's semidefinite programming software (sp) to Octave? I'm starting to work on an .oct file but would like to avoid the work if I can. A S Hodel Dept Elect Eng, Auburn Univ,AL 36849-5201 On leave at NASA Marshall Space Flight Center (205) 544-1426 From help-octave-request@bevo.che.wisc.edu Wed Jul 15 13:12:23 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id NAA15203 for help-octave-dist@bevo.che.wisc.edu; Wed, 15 Jul 1998 13:12:23 -0500 (CDT) Resent-Date: Wed, 15 Jul 1998 13:12:23 -0500 (CDT) Date: Wed, 15 Jul 1998 14:12:20 -0400 (EDT) From: Michael J Ursiak To: help-octave@bevo.che.wisc.edu Subject: passing a matrix to a fortran subroutine Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"Zrdlu82addG.A.GnD.GEPr1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Can anybody help me with the following problem: I am attempting to pass a matrix to fortran subroutine from octave. I have written the following wrapper function based upon your svd.cc function: #include #include #include #include "defun-dld.h" #include "error.h" #include "gripes.h" #include "help.h" #include "mappers.h" #include "oct-obj.h" #include "pr-output.h" #include "utils.h" #include "f77-fcn.h" #include "dDiagMatrix.h" #include "dMatrix.h" #include "ndbleSVD.h" extern "C" { int F77_FCN (xstuff, XSTUFF) ( int&,int&,double*,int&); } DEFUN_DLD (newsvd, args, , "The `newsvd'.") { int z; octave_value arg = args(0); Matrix tmp = arg.matrix_value(); int m = tmp.rows (); int n = tmp.cols (); double *tmp_data = tmp.fortran_vec (); F77_XFCN (xstuff, XSTUFF, (m,n,tmp_data,z)); return octave_value(z); } I get the following error when I try to execute this newvd in octave: octave:2> y = newsvd(h) ld.so.1: octave: fatal: relocation error: symbol not found: __f_open_nv: referenced in /a/kenner/lime/homes/ursiak/octave/newstuff/newsvd.oct Killed I am running Octave, version 2.0.10 on SunOS 5.5.1 sun4m sparc SUNW,SPARCstation-20. From help-octave-request@bevo.che.wisc.edu Wed Jul 15 13:51:07 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id NAA15940 for help-octave-dist@bevo.che.wisc.edu; Wed, 15 Jul 1998 13:51:07 -0500 (CDT) Resent-Date: Wed, 15 Jul 1998 13:51:07 -0500 (CDT) Message-Id: <199807151850.NAA14474@Eng.Auburn.EDU> X-Mailer: Microsoft Outlook Express for Macintosh - 4.0c (197) Date: Wed, 15 Jul 1998 13:43:56 -0500 Subject: Re: passing a matrix to a fortran subroutine From: "A. Scottedward Hodel" To: Michael J Ursiak CC: help-octave@bevo.che.wisc.edu Mime-version: 1.0 X-Priority: 3 Content-type: text/plain; charset="US-ASCII" Content-transfer-encoding: 7bit Resent-Message-ID: <"gIA6-F_GmlE.A.H3D.boPr1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu ---------- >From: Michael J Ursiak >To: help-octave@bevo.che.wisc.edu >Subject: passing a matrix to a fortran subroutine >Date: Wed, Jul 15, 1998, 1:12 PM >extern "C" >{ > int F77_FCN (xstuff, XSTUFF) ( int&,int&,double*,int&); >} > > >I get the following error when I try to execute this newvd in octave: > >octave:2> y = newsvd(h) >ld.so.1: octave: fatal: relocation error: symbol not found: __f_open_nv: >referenced in /a/kenner/lime/homes/ursiak/octave/newstuff/newsvd.oct >Killed I haven't seen this particular error, but here's some things I had to do to get my .oct files working: (1) Make sure gcc was compiled with --enable-shared. (2) Run mkoctfile with sources, not object files, e.g. mkoctfile newsvd.cc xstuff.f not f77 -c xstuff.f mkoctfile newsvd.cc xstuff.o (3) Make sure your LD_LIBRARY_PATH points to your octave libs. If you configured octave with --prefix=prefix_dir, then LD_LIBRARY_PATH should include prefix_dir/lib/octave-2.x.x (your version here) Hope this helps. From help-octave-request@bevo.che.wisc.edu Wed Jul 15 14:51:20 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id OAA16767 for help-octave-dist@bevo.che.wisc.edu; Wed, 15 Jul 1998 14:51:20 -0500 (CDT) Resent-Date: Wed, 15 Jul 1998 14:51:20 -0500 (CDT) From: lash@tellabs.com Message-Id: <199807151950.OAA14144@sunl10.tellabs.com> Subject: Re: passing a matrix to a fortran subroutine To: help-octave@bevo.che.wisc.edu Date: Wed, 15 Jul 1998 14:50:39 -0500 (CDT) In-Reply-To: <199807151850.NAA14474@Eng.Auburn.EDU> from "A. Scottedward Hodel" at Jul 15, 98 01:43:56 pm X-Mailer: ELM [version 2.4 PL24] MIME-Version: 1.0 Content-Type: text/plain; charset=US-ASCII Content-Transfer-Encoding: 7bit Resent-Message-ID: <"sStsrEMPX1D.A.YwD.4gQr1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu > > ---------- > >From: Michael J Ursiak > >To: help-octave@bevo.che.wisc.edu > >Subject: passing a matrix to a fortran subroutine > >Date: Wed, Jul 15, 1998, 1:12 PM > >extern "C" > >{ > > int F77_FCN (xstuff, XSTUFF) ( int&,int&,double*,int&); > >} > > > > > > >I get the following error when I try to execute this newvd in octave: > > > >octave:2> y = newsvd(h) > >ld.so.1: octave: fatal: relocation error: symbol not found: __f_open_nv: > >referenced in /a/kenner/lime/homes/ursiak/octave/newstuff/newsvd.oct > >Killed > > I haven't seen this particular error, but here's some things I had to do > to get my .oct files working: > > (1) Make sure gcc was compiled with --enable-shared. > (2) Run mkoctfile with sources, not object files, e.g. > mkoctfile newsvd.cc xstuff.f > not > f77 -c xstuff.f > mkoctfile newsvd.cc xstuff.o > (3) Make sure your LD_LIBRARY_PATH points to your octave libs. If you > configured octave with --prefix=prefix_dir, then LD_LIBRARY_PATH > should include prefix_dir/lib/octave-2.x.x (your version here) > > Hope this helps. > > I am having a similar problem with octave 2.0.13 under solaris 2.5.1 (SunOS 5.5.1) I wrote a c++ routine to do the inner loop of a routine that calculates tvar (time variance). My first problem was that mkoctfile would dump core if I left it as it was, using g++ to do the linking. If I changed it to use gcc or ld, the .oct file would get generated, but when I tried to use it (or even do a which tvar_oct) octave would report: error: ld.so.1: octave: fatal: relocation error: symbol not found: _t12basic_string3ZcZt18string_char_traits1ZcZt24__default_alloc_template2b0i0.nilRep: referenced in /home/tellabx-5/lash/tdv/tvar_oct.oct: /home/tellabx-5/lash/tdv/tvar_oct.oct `FStvar_oct__Fv' I had assumed that this was due to our tools group not installing gcc correctly (they have recently gone to egcs). Is there any easy way to check if they compiled gcc with --enable-shared? I compiled octave with this compiler and doing an ldd on it shows: sunl10% ldd ~/progs/octave/bin/octave libm.so.1 => /usr/lib/libm.so.1 libdl.so.1 => /usr/lib/libdl.so.1 libsocket.so.1 => /usr/lib/libsocket.so.1 libc.so.1 => /usr/lib/libc.so.1 libnsl.so.1 => /usr/lib/libnsl.so.1 libintl.so.1 => /usr/lib/libintl.so.1 libmp.so.1 => /usr/lib/libmp.so.1 libw.so.1 => /usr/lib/libw.so.1 Which makes me believe that it can handle shared libraries. I do notice that we don't have a shared version of libstdc++ which is where _t12basic_string3ZcZt18string_char_traits1ZcZt24__default_alloc_template2b0i0 is defined. (Don't you just love that C++ name mangling!!!). This symbol does also appear to be in the octave executable file though. I also made sure that LD_LIBRARY_PATH included the octave library directory, and the same problems occur. (Just tried that, thanks for the tip.) Anyway, if anyone has any thoughts about what might be wrong, I would appreciate it. I've attatched the c++ file that I am using below, but the same problems occure with the example files (hello.cc, etc.) Bill Lash lash@tellabs.com =====tvar_oct.cc========================= #include DEFUN_DLD (tvar_oct,args,,"TVAR: tvar for particular observation interval n") { double retval; ColumnVector arg_samples; int arg_N; int arg_n; double inner; int j; int i; int nargin = args.length (); retval = 0.0; if (nargin !=3) { return retval; } arg_samples = args(0).vector_value (); arg_N = NINT(args(1).double_value ()); arg_n = NINT(args(2).double_value ()); for(j=0;j X-Sender: borosy@fabi.vub.ac.be X-Mailer: QUALCOMM Windows Eudora Light Version 3.0.2 (32) Date: Thu, 16 Jul 1998 17:37:09 +0200 To: help-octave@bevo.che.wisc.edu From: =?iso-8859-1?Q?Andr=E1s?= Borosy Subject: path Mime-Version: 1.0 Content-Type: text/plain; charset="us-ascii" Resent-Message-ID: <"vUv4RTk8__P.A.mEG.L5hr1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu What are the counterpart of the following Matlab commands in Ocatve: plot3, surf, peak Many thank, Andras Borosy From help-octave-request@bevo.che.wisc.edu Sun Jul 19 16:22:31 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id QAA25595 for help-octave-dist@bevo.che.wisc.edu; Sun, 19 Jul 1998 16:22:31 -0500 (CDT) Resent-Date: Sun, 19 Jul 1998 16:22:31 -0500 (CDT) Date: Sun, 19 Jul 1998 23:19:32 +0200 (MEST) From: Marco Bravi Reply-To: Marco.Bravi@pcainf1.ing.uniroma1.it To: help-octave@bevo.che.wisc.edu cc: Giuseppe Baffi Subject: FIFO read by Octave does not block as desired Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"IqS4B15sLaH.A.IW.WOms1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Dear Octave users, I would like to use Octave as the control environment for a simulation environment, and the two should communicate over two FIFOs. I am experiencing a misbehaviour of this type: once the input FIFO has been read once (blocking read, as expected), all subsequent reads return immediately, even if the other program is not blocked in the corresponding write. No data is actually input by this read. This is the infinite loop I use on the Octave part: fid1 = -1; fid2 = -1; times = 10; while (times) times--; if (fid1 < 0) fprintf(stderr, "Open gPROMS.pipe..."); fid1 = fopen("../gproms/gPROMS.pipe", "r"); fprintf(stderr, "done.\n"); end fprintf(stderr, "I am about to read..."); [a, b, c, d, e, f] = fscanf(fid1, "%g %g %g %g %g %g\n", "C") fprintf(stderr, "done.\n"); fprintf(stderr, "Going to sleep..."); sleep(15); fprintf(stderr, "awaken.\n"); if (fid2 < 0) fprintf(stderr, "Opening Matlab.pipe..."); fid2 = fopen("../gproms/Matlab.pipe", "w+"); fprintf(stderr, "done.\n"); end fprintf(stderr, "About to send Wf..."); fwrite(fid2, "4\n"); fflush(fid2); fprintf(stderr, "done.\n"); endwhile The other side has been implemented in several ways. For test purposes, a small C program using stdin/stdout and talking to the FIFOS is what I use for debugging all my programs uising FIFOS. What's most puzzling, is that once the first read has been performed correctly (i. e., with blocking behaviour) the fscanf call returns *even if the other program is terminated*, i. e., no other program is writing the FIFO. This appears rather strange to me since I did not open the FIFO for read/write (as it is often seen in order to avoid blocking on open). Could you please suggest an explanation for this and a workaround for correctly syncyng the two concurrent processes with this FIFO read/write system? Thanks a lot in advance. Marco Bravi -- Marco Bravi Marco.Bravi@ingchim.ing.uniroma1.it Dip. Ing. Chimica tel. +39-6-44585587 / 612 v. Eudossiana, 18 fax +39-6-4827453 I-00184 Roma (Italy) From help-octave-request@bevo.che.wisc.edu Mon Jul 20 11:53:05 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id LAA14345 for help-octave-dist@bevo.che.wisc.edu; Mon, 20 Jul 1998 11:53:05 -0500 (CDT) Resent-Date: Mon, 20 Jul 1998 11:53:05 -0500 (CDT) From: Smaranda_Paun@otpp.com X-Lotus-FromDomain: OTPPB To: help-octave@bevo.che.wisc.edu Message-ID: <85256647.005B730D.00@smtp.otpp.com> Date: Mon, 20 Jul 1998 12:51:31 -0400 Subject: Win NT installation problem Mime-Version: 1.0 Content-type: text/plain; charset=us-ascii Resent-Message-ID: <"n8W-Fv0Za2P.A.r0D.xX3s1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu I have tried to install Octave version 2.0.13 on Win NT 4.00 1381 using gnu-win32 b19 to install it. I installed everything on drive D (I don't know if it makes any difference but should I be using C instead of D?) and when I run on bash install-octave I am getting the following errors: installing octave-sh as /octave/bin/octave install-octave: /octave/bin/octave: No such file or directory chmod: not network device >D<. installing src/octave as /octave/bin/octave.bin cp: cannot create regular file `/octave/bin/octave.bin': No such file or directo ry chmod: not network device >D<. installing octave-bug as /octave/bin/octave-bug cp: cannot create regular file `/octave/bin/octave-bug': No such file or directo ry chmod: not network device >D<. /oct/i386-pc-cygwin32 in /octave/libexec/octave/2.0.13 cp: FIND:: No such file or directory chmod: not network device >D<. cp: Parameter: No such file or directory chmod: not network device >D<. cp: format: No such file or directory chmod: not network device >D<. cp: not: No such file or directory chmod: not network device >D<. : No such file or directory chmod: not network device >D<. /mstalling .m files in /octave/share/octave/2.0.13 /m: No such file or directoryhare/octave/2.0.13 chmod: not network device >D<. cp: cannot create regular file `/octave/info/octave.info-4': No such file or dir ectory chmod: not network device >D<. cp: cannot create regular file `/octave/info/octave.info-5': No such file or dir ectory chmod: not network device >D<. cp: cannot create regular file `/octave/info/octave.info-6': No such file or dir ectory chmod: not network device >D<. cp: cannot create regular file `/octave/info/octave.info-7': No such file or dir ectory chmod: not network device >D<. cp: cannot create regular file `/octave/info/octave.info-8': No such file or dir ectory chmod: not network device >D<. cp: cannot create regular file `/octave/info/octave.info-9': No such file or dir ectory chmod: not network device >D<. installing man page in /octave/man/man1 cp: cannot create regular file `/octave/man/man1/octave.1': No such file or dire ctory chmod: not network device >D<. There's more of them, but I just copied some. Does anyone have any idea what went wrong? Thanks a lot, Smaranda. From help-octave-request@bevo.che.wisc.edu Mon Jul 20 18:55:59 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id SAA19108 for help-octave-dist@bevo.che.wisc.edu; Mon, 20 Jul 1998 18:55:59 -0500 (CDT) Resent-Date: Mon, 20 Jul 1998 18:55:59 -0500 (CDT) Message-ID: <19980720185743.10663@lymatech.com> Date: Mon, 20 Jul 1998 18:57:43 -0500 From: Suttipan Limanond To: help-octave@bevo.che.wisc.edu Subject: [MAILER-DAEMON: Returned mail: Host unknown (Name server: lymatech.com: host not found)] Mime-Version: 1.0 Content-Type: text/plain; charset=us-ascii X-Mailer: Mutt 0.84e Resent-Message-ID: <"tiOodv88qbF.A.2rD.Pk9s1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu -----Forwarded message from Mail Delivery Subsystem ----- Received: from localhost (localhost) by dugdae.lymatech.com (8.8.7/8.8.7) with internal id RAA00489; Mon, 20 Jul 1998 17:30:34 -0500 (CDT) (envelope-from MAILER-DAEMON) Date: Mon, 20 Jul 1998 17:30:34 -0500 (CDT) From: Mail Delivery Subsystem Message-Id: <199807202230.RAA00489@dugdae.lymatech.com> To: suttipan MIME-Version: 1.0 Content-Type: multipart/report; report-type=delivery-status; boundary="RAA00489.900973834/dugdae.lymatech.com" Subject: Returned mail: Host unknown (Name server: lymatech.com: host not found) Auto-Submitted: auto-generated (failure) This is a MIME-encapsulated message --RAA00489.900973834/dugdae.lymatech.com The original message was received at Mon, 20 Jul 1998 15:24:19 -0500 (CDT) from suttipan@localhost ----- The following addresses had permanent fatal errors ----- octave@lymatech.com ----- Transcript of session follows ----- 550 octave@lymatech.com... Host unknown (Name server: lymatech.com: host not found) --RAA00489.900973834/dugdae.lymatech.com Content-Type: message/delivery-status Reporting-MTA: dns; dugdae.lymatech.com Arrival-Date: Mon, 20 Jul 1998 15:24:19 -0500 (CDT) Final-Recipient: RFC822; octave@lymatech.com Action: failed Status: 5.1.2 Remote-MTA: DNS; lymatech.com Last-Attempt-Date: Mon, 20 Jul 1998 17:30:34 -0500 (CDT) --RAA00489.900973834/dugdae.lymatech.com Content-Type: message/rfc822 Return-Path: Received: (from suttipan@localhost) by dugdae.lymatech.com (8.8.7/8.8.7) id PAA00280; Mon, 20 Jul 1998 15:24:19 -0500 (CDT) (envelope-from suttipan) Message-ID: <19980720152418.08961@lymatech.com> Date: Mon, 20 Jul 1998 15:24:18 -0500 From: Suttipan Limanond To: octave@lymatech.com Subject: ODE solver Mime-Version: 1.0 Content-Type: text/plain; charset=us-ascii X-Mailer: Mutt 0.84e I wonder if we can pass extra parameters to the user-supplied function (.oct file) associated with `lsode', other than putting everything inside that function. Thanks, Suttipan Limanond. --RAA00489.900973834/dugdae.lymatech.com-- -----End of forwarded message----- From help-octave-request@bevo.che.wisc.edu Thu Jul 23 10:42:58 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id KAA02327 for help-octave-dist@bevo.che.wisc.edu; Thu, 23 Jul 1998 10:42:58 -0500 (CDT) Resent-Date: Thu, 23 Jul 1998 10:42:58 -0500 (CDT) Date: Thu, 23 Jul 1998 23:41:09 +0800 (WST) From: Mike Alder To: help-octave@bevo.che.wisc.edu cc: Mike Alder Subject: RedHat Linux Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"Z3MHWyJNZ3F.A.mqF.Bo1t1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Hi, I have downloaded what seems to be specified as LINUX binaries. I found, after some misdirection, the install instructions. This failed. I logged in as root and it installed. Sort of. I moved to the /usr/local/bin directory and there was an executable called octave. I typed `octave' and got the message: error: Segmentation fault -- stopping myself... error: attempted clean up apparently failed -- aborting... Aborted So far all I seem to have accomplished is to clog up my hard disk with a lot of files that don't do anything. The download was from your homepages. Any suggestions? (Other than buying MATLAB). mike alder. From help-octave-request@bevo.che.wisc.edu Thu Jul 23 16:23:40 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id QAA13424 for help-octave-dist@bevo.che.wisc.edu; Thu, 23 Jul 1998 16:23:40 -0500 (CDT) Resent-Date: Thu, 23 Jul 1998 16:23:40 -0500 (CDT) Date: Thu, 23 Jul 1998 16:22:35 -0500 (CDT) From: Jesse Bennett Reply-To: Jesse Bennett To: help-octave@bevo.che.wisc.edu Subject: C language interface to octave (eng* function library) Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"bICd-FDyW-P.A.UOF.cn6t1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Hi, Is there any octave support (or ongoing development) for an Octave interface library? I find that Matlab is very effective as a user interface for much of the R&D that I am involved in, and I am interested in the possibility of using octave for this work. In other words, is there any compatibility with the libmat.a functions eng*? E.g., > nm /apps/matlab/extern/lib/sun4/libmat.a | grep eng fengintf.o_sun4: U _engClose U _engEvalString U _engGetFull U _engGetMatrix U _engOpen U _engOutputBuffer U _engPutFull U _engPutMatrix 000000cc T _engclose_ 0000002c T _engevalstring_ [ ... ] Thanks in advance, Jesse Bennett From help-octave-request@bevo.che.wisc.edu Thu Jul 23 17:22:58 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id RAA23871 for help-octave-dist@bevo.che.wisc.edu; Thu, 23 Jul 1998 17:22:58 -0500 (CDT) Resent-Date: Thu, 23 Jul 1998 17:22:58 -0500 (CDT) Message-Id: <199807232222.RAA17343@Eng.Auburn.EDU> X-Mailer: Microsoft Outlook Express for Macintosh - 4.0c (197) Date: Thu, 23 Jul 1998 17:21:25 -0500 Subject: Re: C language interface to octave (eng* function library) From: "A. Scottedward Hodel" To: Jesse Bennett , help-octave@bevo.che.wisc.edu Mime-version: 1.0 X-Priority: 3 Content-type: text/plain; charset="US-ASCII" Content-transfer-encoding: 7bit Resent-Message-ID: <"WZy_Zf02oiH.A.Qs.Cf7t1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Are you aware of the .oct file support for dynamically linked functions? There are several examples of these in the octave distribution in examples/oregonator.cc examples/hello.cc and in the src/DLD-FUNCTIONS directory (with some straightforward editing of the include file paths for mkoctfile). I've used this resource to develop an update of the qz function in Octave to solve generalized eigenvalues (to appear in a later release; requires a few extra fortran routines, etc.) I've also implemented some examples of .oct files working with structures in the octave controls toolbox (see ftp://ftp.eng.auburn.edu/pub/hodel), which is also undergoing some rewriting at this time (maintaining compatibility with previous releases!). I found it somewhat less painful to write .oct files than I did mex-files (although I had more experience working with this sort of interface, so it's not a fair comparison). The .oct files have access to the full Octave library, and so you'll already have access to much of BLAS, LAPACK, etc. when writing your C++ code (which can easily act as a wrapper to call your own C-code if you wish). Those specific functions you list in your note are proprietary to MATLAB; I don't think anyone has written direct interfaces to them. Backward compatibility and portability among versions of matlab (4,5) may be a problem to hinder such an implementation, since Octave ranges over a wide range of target machines. A S Hodel Dept Elect Eng, Auburn Univ,AL 36849-5201 On leave at NASA Marshall Space Flight Center (205) 544-1426 ---------- >From: Jesse Bennett >To: help-octave@bevo.che.wisc.edu >Subject: C language interface to octave (eng* function library) >Date: Thu, Jul 23, 1998, 4:22 PM [edited] >Is there any octave support (or ongoing development) for an Octave >interface library? I find that Matlab is very effective as a user >interface for much of the R&D that I am involved in, and I am interested >in the possibility of using octave for this work. In other words, is >there any compatibility with the libmat.a functions eng*? E.g., From help-octave-request@bevo.che.wisc.edu Thu Jul 23 17:30:46 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id RAA17847 for help-octave-dist@bevo.che.wisc.edu; Thu, 23 Jul 1998 17:30:46 -0500 (CDT) Resent-Date: Thu, 23 Jul 1998 17:30:46 -0500 (CDT) Date: Fri, 24 Jul 1998 10:34:52 +1200 From: Ben Skelton Subject: Re: RedHat Linux To: help-octave@bevo.che.wisc.edu Reply-to: bcs25@cantva.canterbury.ac.nz Message-id: <98072410430502.09962@avalon.elec.canterbury.ac.nz> MIME-version: 1.0 X-Mailer: KMail [version 0.7.9] Content-type: text/plain Content-transfer-encoding: 8bit References: Resent-Message-ID: <"2-FQB39v1LM.A.z1F.Wm7t1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Hi Mike There are some big problems with RH binaries if you install a binary which was compiled with a different compiler than the compiler system you have installed on your system...most likely a pre RH5.1 binary installed on a RH5.1 system. This is a gerneral RH problem and not I suspect unique to the octave binary that you downloaded---I may be wrong I dont know what you downloaded. See if you can get a hold of an SRPM and compile it yourself...I think that is the safest until RH get their act together and sort out their compilier woes. Hope this helps. --Ben On Fri, 24 Jul 1998, Mike Alder wrote: >Hi, >I have downloaded what seems to be specified as LINUX >binaries. I found, after some misdirection, the install >instructions. This failed. I logged in as root and it >installed. Sort of. I moved to the /usr/local/bin directory >and there was an executable called octave. I typed `octave' >and got the message: > >error: Segmentation fault -- stopping myself... >error: attempted clean up apparently failed -- aborting... >Aborted > >So far all I seem to have accomplished is to clog up my >hard disk with a lot of files that don't do anything. > >The download was from your homepages. > >Any suggestions? (Other than buying MATLAB). > >mike alder. From help-octave-request@bevo.che.wisc.edu Fri Jul 24 02:19:17 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id CAA31264 for help-octave-dist@bevo.che.wisc.edu; Fri, 24 Jul 1998 02:19:17 -0500 (CDT) Resent-Date: Fri, 24 Jul 1998 02:19:17 -0500 (CDT) Date: Fri, 24 Jul 1998 09:02:34 +0200 (MET DST) From: Jeroen Schaap To: help-octave@bevo.che.wisc.edu Subject: Segmentation fault. Libraries seem to be OK. What now? Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"fzGTcrNKXuF.A.9fH.1VDu1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Hi I've a small linux-box and really like to install octave. It's a 468-dx50 with 8 meg internal memory. A little bit too small to compile octave. Yes, I know that i could add a swapfile, but if the binaries could run that would be fine to me. I got me this file from che.wisc.edu: octave-2.0.13-i686-pc-linux-gnulibc1.tar.gz Ran sh ./install-octave It installed... But trying to run it results in this: error: Segmentation fault -- stopping myself... error: attempted clean up apparently failed -- aborting... IOT trap/Abort (core dumped) I asked in the local newsgroup... they said to check the libraries, since it was a point in the readme/install's... so: rulffh:/usr/local/src# ldd /usr/local/bin/octave liboctinterp.so => /usr/local/lib/octave-2.0.13/liboctinterp.so liboctave.so => /usr/local/lib/octave-2.0.13/liboctave.so libcruft.so => /usr/local/lib/octave-2.0.13/libcruft.so libm.so.5 => /lib/libm.so.5.0.9 libncurses.so.3.0 => /lib/libncurses.so.3.0 libdl.so.1 => /lib/libdl.so.1.8.2 libstdc++.so.27 => /usr/lib/libstdc++.so.27.1.4 libc.so.5 => /lib/libc.so.5.4.44 seems to be OK. A little bit strange that the libstdc++ points to .27 since I have .2.8.0: ( i guess .27 means 2.7..... rulffh:/usr/local/src# ls /usr/lib/libst* /usr/lib/libstdc++.a /usr/lib/libstdc++.so.2.8.0* /usr/lib/libstdc++.so@ /usr/lib/libstdc++.so.27@ /usr/lib/libstdc++.so.2.8@ /usr/lib/libstdc++.so.27.1.4* So is there a way to get octave linking to the 2.8 libs? Or could i try relinking so.27 to so.2.8.0? Thanks!! Jeroen Jeroen Schaap.............I was dreaming of guitarnotes that would irritate Homepage: ..| ^|^ |...an executive kind Keywords: Guitars, Linux, Mac and SCN...........\_`-'_/..............of guy Tel: (0)20-5667127................................| |...........Frank Zappa From help-octave-request@bevo.che.wisc.edu Sat Jul 25 20:07:12 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id UAA27884 for help-octave-dist@bevo.che.wisc.edu; Sat, 25 Jul 1998 20:07:12 -0500 (CDT) Resent-Date: Sat, 25 Jul 1998 20:07:12 -0500 (CDT) Message-Id: <199807260107.VAA002.71@wkim> MIME-Version: 1.0 Date: Sat, 25 Jul 98 20:59:22 -0400 From: "Wonkoo Kim" To: "Octave Users" Reply-To: "Wonkoo Kim" Subject: diary X-Mailer: Ultimedia Mail/2 Lite, IBM T. J. Watson Research Center Content-ID: <270_86_4_901414762> Content-Type: text/plain; charset="US-ASCII" Content-Transfer-Encoding: 7bit Resent-Message-ID: <"qLiRq6tN9GN.A.SJD.AFou1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu I'm using Octave 2.0.13 for OS/2. I noticed a problem that "diary" doesn't work well if PAGER was invoked (due to a long Octave output screen). Paged screen output were not saved in diary file. Is there any way to redirect all octave screen outputs to diary file regardless of the length of output? Thanks. //-------------------------------------------------------------------- // Wonkoo Kim (wkim+@pitt.edu) From help-octave-request@bevo.che.wisc.edu Sat Jul 25 20:19:09 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id UAA06836 for help-octave-dist@bevo.che.wisc.edu; Sat, 25 Jul 1998 20:19:09 -0500 (CDT) Resent-Date: Sat, 25 Jul 1998 20:19:09 -0500 (CDT) Message-Id: <199807260119.VAA002.86@wkim> MIME-Version: 1.0 Date: Sat, 25 Jul 98 21:13:10 -0400 From: "Wonkoo Kim" To: "Octave Users" Reply-To: "Wonkoo Kim" Subject: diary X-Mailer: Ultimedia Mail/2 Lite, IBM T. J. Watson Research Center Content-ID: <270_86_4_901415590> Content-Type: text/plain; charset="US-ASCII" Content-Transfer-Encoding: 7bit Resent-Message-ID: <"0z35F5MJt7N.A.9ZH.NQou1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu I think the problem is solved. I had put diary "file.log" at the beginning of a .m file, and diary off at the end of the .m file. However, it looks like I have to remove "diary off" from the .m file, and I have to issue it directly from Octave session. I don't know why, but it worked that way. Thanks. //-------------------------------------------------------------------- // Wonkoo Kim (wkim+@pitt.edu) From help-octave-request@bevo.che.wisc.edu Mon Jul 27 17:40:18 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id RAA31888 for help-octave-dist@bevo.che.wisc.edu; Mon, 27 Jul 1998 17:40:18 -0500 (CDT) Resent-Date: Mon, 27 Jul 1998 17:40:18 -0500 (CDT) Date: Mon, 27 Jul 1998 18:38:13 -0400 From: "Dimitrie O. Paun" X-Sender: dimi@qew.cs To: help-octave@bevo.che.wisc.edu Subject: binomial normal cdf Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"2Q2wvrUA-oF.A.srH.SHQv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Hello, I am looking for a Octave/Matlab function that can return the binomial normal cdf. Does anyone has such a beast? Thanks, Dimi -- Dimitrie O. Paun | Graduate Student, Computer Science | Friends don't let friends boot NT! University of Toronto | From help-octave-request@bevo.che.wisc.edu Mon Jul 27 18:32:58 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id SAA30973 for help-octave-dist@bevo.che.wisc.edu; Mon, 27 Jul 1998 18:32:58 -0500 (CDT) Resent-Date: Mon, 27 Jul 1998 18:32:58 -0500 (CDT) Message-ID: X-Mailer: XFMail 1.3-alpha-031298 [p0] on Linux X-Priority: 3 (Normal) Content-Type: text/plain; charset=iso-8859-1 Content-Transfer-Encoding: 8bit MIME-Version: 1.0 In-Reply-To: Date: Tue, 28 Jul 1998 00:31:03 +0100 (BST) Reply-To: Ted.Harding@nessie.mcc.ac.uk From: (Ted Harding) To: "Dimitrie O. Paun" Subject: RE: binomial normal cdf Cc: help-octave@bevo.che.wisc.edu Resent-Message-ID: <"P5DMIRPPHK.A.A1H.q4Qv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu On 27-Jul-98 Dimitrie O. Paun wrote: > I am looking for a Octave/Matlab function that can return the binomial > normal cdf. Does anyone has such a beast? For the binomial probability function try function P=binompdf(n,p) # For integer n, real p (0 <= p <= 1), # returns the entire binomial distribution on (0,1,2,...,n) # as a column vector. if (( p<0 ) | ( p>1) ) error ( "p is out of range" ); endif if ( p == 0 ) P = zeros(n+1,1); P(1) = 1; elseif ( p == 1 ) P = zeros(n+1,1); P(n+1) = 1; else r = (0:n)'; P = exp(lgamma(n+1)-lgamma(r+1)-lgamma(n-r+1)+r*log(p)+(n-r)*log(1-p)); P = P/sum(P); endif endfunction and for the binomial CDF try function P=binomCDF(n,p) # For integer n, real p (0 <= p <= 1), # returns the entire binomial cumulative distribution on (0,1,2,...,n) # as a column vector. P = cumsum(binompdf(n,p)); endfunction For the Normal CDF try function p = normCDF(x); # function p = normCDF(x) = int from (-inf) to x of normpdf(t) dt s=sqrt(2); p=0.5*(1 + erf(x/s)); endfunction I don't know what would be meant by "binomial normal CDF" unless one of the above. Best wishes, Ted. -------------------------------------------------------------------- E-Mail: (Ted Harding) Date: 28-Jul-98 Time: 00:31:03 -------------------------------------------------------------------- From help-octave-request@bevo.che.wisc.edu Mon Jul 27 18:43:26 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id SAA27736 for help-octave-dist@bevo.che.wisc.edu; Mon, 27 Jul 1998 18:43:26 -0500 (CDT) Resent-Date: Mon, 27 Jul 1998 18:43:26 -0500 (CDT) Date: Mon, 27 Jul 1998 19:35:51 -0400 From: "Dimitrie O. Paun" X-Sender: dimi@qew.cs Reply-To: "Dimitrie O. Paun" To: help-octave@bevo.che.wisc.edu Subject: bivariate normal cdf (sorry, NOT binomial;) In-Reply-To: <002a01bdb9b4$9ef21900$3ec8c58f@thuban.llumc.edu> Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"cBdNTvjIx_D.A.6MF.eCRv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu I am sorry for the confusion, I don't know what I was thinking...;( I was looking for the bivariate normal cdf. TIA, Dimi -- Dimitrie O. Paun | Graduate Student, Computer Science | Friends don't let friends boot NT! University of Toronto | From help-octave-request@bevo.che.wisc.edu Tue Jul 28 03:47:47 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id DAA18640 for help-octave-dist@bevo.che.wisc.edu; Tue, 28 Jul 1998 03:47:47 -0500 (CDT) Resent-Date: Tue, 28 Jul 1998 03:47:47 -0500 (CDT) >Received: from localhost by csnio.ren.nic.in with smtp id m0z11GF-0001YuC (Debian Smail-3.2 1996-Jul-4 #2); Tue, 28 Jul 1998 09:51:15 +0530 (IST) Date: Tue, 28 Jul 1998 09:51:15 +0530 (IST) From: "Mr.Jaydeep Pattanaik IODC NIO" To: help-octave@bevo.che.wisc.edu Subject: Re: netcdf (fwd) In-Reply-To: Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"HrdTEuWl9bN.A.7r.zAZv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Dear all, I am a first timer in Octave. Want to know whether 1. I can directly use netcdf data in octave. 2. Fit a second order polynomial to my data. Any hints will be greatly appreciated. Jaydeep ---------------------------------------------------------------------- > Jaydeep Pattanaik, |Phone - 091-0832-229900 Ext- 4241 > Indian Oceanographic Data Centre, |Fax - 091-0832-223340 > National Institute of Oceanography, | 091-0832-229102 > Dona Paula, Goa - 403 004, INDIA |e-mail: jaydeep@csnio.ren.nic.in > | jaydeep@darya.nio.org > --------------------------------------------------------------------- From help-octave-request@bevo.che.wisc.edu Tue Jul 28 10:17:11 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id KAA31612 for help-octave-dist@bevo.che.wisc.edu; Tue, 28 Jul 1998 10:17:11 -0500 (CDT) Resent-Date: Tue, 28 Jul 1998 10:17:11 -0500 (CDT) Message-ID: <35BDEB4B.5A8091E7@ffa.se> Date: Tue, 28 Jul 1998 17:16:27 +0200 From: Daniel Tourde Organization: The Aeronautical Research Institute of Sweden X-Mailer: Mozilla 3.04Gold (X11; I; Linux 2.0.33 i686) MIME-Version: 1.0 To: jwe@bevo.che.wisc.edu, help-octave@bevo.che.wisc.edu Subject: Eigenvalues computation Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Resent-Message-ID: <"HsNpI0rBtbL.A.H3B.3tev1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Hi ! I would like to trace the variation of the eigenvalues of a matrix D(k) in function of k. Until now I have considered that Octave kept the same order for my eigenvalues but in view of some strange results I have I'm wondering if at a certain step of my computation the order does not change (lambda8 <-> lambda9 for instance). So, here is my question : Does octave classify the eigenvalues of a matrix ? If yes, how ? If yes, how can I force Octave to keep this classification all along my computation of lambda(k) x = D(k) x ? If no, what can I do ? .... Thanks for any answer. Best regards Daniel Tourde -- ********************************************************************* Daniel TOURDE E-mail : ted@ffa.se The Aeronautical Research Institute of Sweden Tel : +46 8 634 13 44 P.O. Box 11021 S-161 11 BROMMA, Sweden Fax : +46 8 25 34 81 ********************************************************************* From help-octave-request@bevo.che.wisc.edu Tue Jul 28 10:42:36 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id KAA08084 for help-octave-dist@bevo.che.wisc.edu; Tue, 28 Jul 1998 10:42:36 -0500 (CDT) Resent-Date: Tue, 28 Jul 1998 10:42:36 -0500 (CDT) Message-Id: <35BDF151.F843703E@bmw.de> Date: Tue, 28 Jul 1998 17:42:09 +0200 From: Daniel Heiserer X-Mailer: Mozilla 4.04 [en] (X11; I; IRIX64 6.2 IP28) Mime-Version: 1.0 To: Daniel Tourde , "help-octave@bevo.che.wisc.edu" Subject: Re: Eigenvalues computation References: <35BDEB4B.5A8091E7@ffa.se> Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Resent-Message-ID: <"E3UYVae9-sG.A.noG.sFfv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Daniel Tourde wrote: > > Hi ! > > I would like to trace the variation of the eigenvalues of a matrix D(k) > in function of k. Until now I have considered that Octave kept the same > order for my eigenvalues but in view of some strange results I have I'm > wondering if at a certain step of my computation the order does not > change (lambda8 <-> lambda9 for instance). > > So, here is my question : > Does octave classify the eigenvalues of a matrix ? > If yes, how ? > If yes, how can I force Octave to keep this classification all along my > computation of lambda(k) x = D(k) x ? > If no, what can I do ? .... > > Thanks for any answer. Best regards Hi. What is an eigenvalue (in a physical sense)? It's not only a mathematical solution of solving the saekular-equation det(A-L1)=0!. In fact an eigenvalue describes some geometry ==> eigenvector. If you are looking on eigenvalues only you probably assume that the eigenvectors don't change or change only a little so you could map the old values to the new values. But anyway they are not equivalent, because their "shapes" change also (even slightly). So if you want to trace eigenvalues don't look at the values first, look at their shapes (eigenvectors). [value_1,shape_1=eig(mymatrix(x1)); [value_2,shape_2=eig(mymatrix(x2)); just make the scalar product of the new shapes with the old shapes: mac=shape_1*shape_2'; the mac (modal assurance criterion) then tells you if (if their index is close to 1) and how (look at the indices) you can correlate the old and the new shapes. Bye correlating the shapes you can look at their eigenvalues and you are done. Have a nice day. Bye Daniel ----- -------------------------------------------------------------- Daniel Heiserer, BMW AG, Knorrstrasse 147, 80788 Muenchen Abteilung EK-20 Tel.: 089-382-21187, Fax.: 089-382-42820 mailto:daniel.heiserer@bmw.de From help-octave-request@bevo.che.wisc.edu Tue Jul 28 11:13:23 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id LAA09691 for help-octave-dist@bevo.che.wisc.edu; Tue, 28 Jul 1998 11:13:23 -0500 (CDT) Resent-Date: Tue, 28 Jul 1998 11:13:23 -0500 (CDT) Message-Id: <199807281613.LAA05073@Eng.Auburn.EDU> X-Mailer: Microsoft Outlook Express for Macintosh - 4.0c (197) Date: Tue, 28 Jul 1998 11:12:15 -0500 Subject: Re: Eigenvalues computation From: "A. Scottedward Hodel" To: Daniel Tourde , jwe@bevo.che.wisc.edu, help-octave@bevo.che.wisc.edu Mime-version: 1.0 X-Priority: 3 Content-type: text/plain; charset="US-ASCII" Content-transfer-encoding: 7bit Resent-Message-ID: <"bxTWU8Xm2mG.A.IeB.jifv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Octave's eig routine does not sort eigenvalues; if you're trying to track eigenvalue variation (which I've also been doing for my work here), you may wish to use the sortcom routine that I wrote for the octave controls toolbox ftp://ftp.eng.auburn.edu/pub/hodel. There's an updated version that I need to put out that will take the same calling syntax as the sort routine; the current version only returns the sorted values and doesn't always put conjugate pairs in the same order. I'll have the new version out later today. A S Hodel Dept Elect Eng, Auburn Univ,AL 36849-5201 On leave at NASA Marshall Space Flight Center (205) 544-1426 ---------- >From: Daniel Tourde >To: jwe@bevo.che.wisc.edu, help-octave@bevo.che.wisc.edu >Subject: Eigenvalues computation >Date: Tue, Jul 28, 1998, 10:16 AM > >Hi ! > >I would like to trace the variation of the eigenvalues of a matrix D(k) >in function of k. Until now I have considered that Octave kept the same >order for my eigenvalues but in view of some strange results I have I'm >wondering if at a certain step of my computation the order does not >change (lambda8 <-> lambda9 for instance). > >So, here is my question : >Does octave classify the eigenvalues of a matrix ? >If yes, how ? >If yes, how can I force Octave to keep this classification all along my >computation of lambda(k) x = D(k) x ? >If no, what can I do ? .... > >Thanks for any answer. Best regards > > Daniel Tourde >-- >********************************************************************* >Daniel TOURDE E-mail : ted@ffa.se >The Aeronautical Research Institute of Sweden Tel : +46 8 634 13 44 >P.O. Box 11021 S-161 11 BROMMA, Sweden Fax : +46 8 25 34 81 >********************************************************************* > From help-octave-request@bevo.che.wisc.edu Tue Jul 28 12:53:56 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id MAA27412 for help-octave-dist@bevo.che.wisc.edu; Tue, 28 Jul 1998 12:53:56 -0500 (CDT) Resent-Date: Tue, 28 Jul 1998 12:53:56 -0500 (CDT) From: Dirk Laurie Message-Id: <199807281752.TAA31620@calvyn.puk.ac.za> Subject: Re: Eigenvalues computation To: help-octave@bevo.che.wisc.edu Date: Tue, 28 Jul 1998 19:52:43 +0200 (SAT) Reply-to: wskdpl@pukrs12.puk.ac.za Content-Type: text Resent-Message-ID: <"V3lFH2C7D0C.A.sSE.0Ahv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu In response to some discussion by: A. Scottedward Hodel, Daniel Heiserer, Daniel Tourde (proposer) on eigenvalues and eigenvectors as continuous functions of a matrix. The question is a comfortable one only when the matrix is normal (i.e. A'*A = A*A'). In that case eigenvalues are very well-behaved functions of the matrix. If a normal matrix A is close to another normal matrix B which has eigenvector matrix Q (in this case known to be unitary, i.e. inv(Q)=Q') then Q'*A*Q will be close to a diagonal matrix, and in fact its diagonal will be extra close to the eigenvalues of A. More precisely, if A and B differ by O(epsilon) then the off-diagonal part of diag(Q'*A*Q) is of magnitude O(epsilon) and diag(Q'*A*Q) differs from the eigenvalues of A by O(epsilon^2). When A is not normal, the eigenvalues are still continuous functions of A, but these functions may be extremely ill-behaved. Even in the case where A is normal, it can be numerically delicate to trace an eigenvector as a function of A when A has nearly equal eigenvalues. If you use the method [S,X]=eig(A); [y,p]=sort(diag(X)); S=S(:,p); X=X(p,p); you get the eigenvalues sorted but the corresponding eigenvectors of A and B may nevertheless differ drastically. These issues are explained clearly in the classical books by J H Wilkinson: "The Algebraic Eigenvalue Problem" and "Rounding Errors in Algebraic Processes", as well as later tomes like "Accuracy and Stability of Numerical Algorithms" by Nicholas J Higham. -- Dirk Laurie From help-octave-request@bevo.che.wisc.edu Tue Jul 28 13:38:44 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id NAA25245 for help-octave-dist@bevo.che.wisc.edu; Tue, 28 Jul 1998 13:38:44 -0500 (CDT) Resent-Date: Tue, 28 Jul 1998 13:38:44 -0500 (CDT) Message-Id: <199807281838.NAA14238@Eng.Auburn.EDU> X-Mailer: Microsoft Outlook Express for Macintosh - 4.0c (197) Date: Tue, 28 Jul 1998 13:37:53 -0500 Subject: memory exhausted From: "A. Scottedward Hodel" To: help-octave@bevo.che.wisc.edu Mime-version: 1.0 X-Priority: 3 Content-type: text/plain; charset="ISO-8859-1" Content-transfer-encoding: 7bit Resent-Message-ID: <"icUsBSyQ4sC.A.WtC.0qhv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu I ran the following code on 1 172x172 matrix: # truncate acd to 6 digits of precision... for ii=1:rows(acd) for jj = 1:columns(acd) eval(["acd(ii,jj) = ",sprintf("%12.6E",acd(ii,jj));]); endfor endfor After some time into the loop, Octave repeatedly displays an error message: "memory exhausted --- trying to return to prompt" I recall some discussion of this from a month or so ago, but I don't recall if anyone encountered a solution. It appears that what's happening is that free memory is getting chopped up into little pieces and the memory manager eventually can't allocate a new acd each time through the loop. What this means is that the total amount of computation I can do in a single octave run becomes quite limited; is there some way we can "help" the memory manager? [The import of the question goes far beyond the code above, which is being used to try to duplicate results from another package; I've had to modify the OCST code I'm writing so that the number of calls to "new" is kept as low as possible, otherwise some code dies after 40-50 iterations on systems of this size...] Details: Octave-2.1.7 on SunOS; dunamic loading enabled; gcc-2.8.1. uname -s -n -m -v -r SunOS eds2 5.6 Generic sun4u A S Hodel Dept Elect Eng, Auburn Univ,AL 36849-5201 On leave at NASA Marshall Space Flight Center (205) 544-1426 From help-octave-request@bevo.che.wisc.edu Tue Jul 28 14:44:12 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id OAA08839 for help-octave-dist@bevo.che.wisc.edu; Tue, 28 Jul 1998 14:44:12 -0500 (CDT) Resent-Date: Tue, 28 Jul 1998 14:44:12 -0500 (CDT) Message-ID: X-Mailer: XFMail 1.2 [p0] on HPUX X-Priority: 3 (Normal) Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 8bit MIME-Version: 1.0 In-Reply-To: <199807281838.NAA14238@Eng.Auburn.EDU> Date: Tue, 28 Jul 1998 21:43:38 +0200 (METDST) Reply-To: Dietz@AIA.RWTH-Aachen.De Organization: Aerodynamisches Institut der RWTH Aachen From: Guido Dietz To: "A. Scottedward Hodel" Subject: RE: memory exhausted Cc: help-octave@bevo.che.wisc.edu Resent-Message-ID: <"GwnTnRcSuhM.A.OdB.Loiv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu On 28-Jul-98 A. Scottedward Hodel wrote: > I ran the following code on 1 172x172 matrix: > ># truncate acd to 6 digits of precision... > for ii=1:rows(acd) > for jj = 1:columns(acd) > eval(["acd(ii,jj) = ",sprintf("%12.6E",acd(ii,jj));]); > endfor > endfor What about: acd=round(acd*1e6)/1e6; Should a) work and b) be faster ... G. From help-octave-request@bevo.che.wisc.edu Tue Jul 28 14:55:24 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id OAA11727 for help-octave-dist@bevo.che.wisc.edu; Tue, 28 Jul 1998 14:55:24 -0500 (CDT) Resent-Date: Tue, 28 Jul 1998 14:55:24 -0500 (CDT) Message-Id: <199807281954.OAA19134@Eng.Auburn.EDU> X-Mailer: Microsoft Outlook Express for Macintosh - 4.0c (197) Date: Tue, 28 Jul 1998 14:53:52 -0500 Subject: Re: memory exhausted From: "A. Scottedward Hodel" To: dietz@aia.rwth-aachen.de CC: help-octave@bevo.che.wisc.edu Mime-version: 1.0 X-Priority: 3 Content-type: text/plain; charset="US-ASCII" Content-transfer-encoding: 7bit Resent-Message-ID: <"5jcD_PTEd3.A.1X.syiv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Thanks, but it's not equivalent if you have widely varying coefficients... Which I do. (from 1e-5 to 1e+7, balancing inappropriate in this appl. for non-numerical reasons.). A S Hodel Dept Elect Eng, Auburn Univ,AL 36849-5201 On leave at NASA Marshall Space Flight Center (205) 544-1426 ---------- >From: Guido Dietz >To: "A. Scottedward Hodel" >Cc: help-octave@bevo.che.wisc.edu >Subject: RE: memory exhausted >Date: Tue, Jul 28, 1998, 2:43 PM > >On 28-Jul-98 A. Scottedward Hodel wrote: >> I ran the following code on 1 172x172 matrix: >> >># truncate acd to 6 digits of precision... >> for ii=1:rows(acd) >> for jj = 1:columns(acd) >> eval(["acd(ii,jj) = ",sprintf("%12.6E",acd(ii,jj));]); >> endfor >> endfor >What about: >acd=round(acd*1e6)/1e6; >Should a) work and b) be faster ... > >G. > From help-octave-request@bevo.che.wisc.edu Wed Jul 29 01:14:40 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id BAA08694 for help-octave-dist@bevo.che.wisc.edu; Wed, 29 Jul 1998 01:14:40 -0500 (CDT) Resent-Date: Wed, 29 Jul 1998 01:14:40 -0500 (CDT) From: Dirk Laurie Message-Id: <199807290613.IAA01800@calvyn.puk.ac.za> Subject: Re: memory exhausted To: help-octave@bevo.che.wisc.edu Date: Wed, 29 Jul 1998 08:13:23 +0200 (SAT) In-Reply-To: <199807281954.OAA19134@Eng.Auburn.EDU> from "A. Scottedward Hodel" at Jul 28, 98 02:53:52 pm Reply-to: wskdpl@pukrs12.puk.ac.za Content-Type: text Resent-Message-ID: <"PFcePFpIscI.A.rkE.Q3rv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu A. Scottedward Hodel wrote: > > Thanks, but it's not equivalent if you have widely varying coefficients... > Which I do. (from 1e-5 to 1e+7, balancing inappropriate in this appl. for > non-numerical reasons.). > > ---------- > >From: Guido Dietz > >Date: Tue, Jul 28, 1998, 2:43 PM > > >On 28-Jul-98 A. Scottedward Hodel wrote: > >> I ran the following code on 1 172x172 matrix: > >> > >># truncate acd to 6 digits of precision... > >> for ii=1:rows(acd) > >> for jj = 1:columns(acd) > >> eval(["acd(ii,jj) = ",sprintf("%12.6E",acd(ii,jj));]); > >> endfor > >> endfor > >What about: > >acd=round(acd*1e6)/1e6; > >Should a) work and b) be faster ... > > If the only purpose of that computation is to simulate floating-point with lower precision, then acd=acd.*(1+1e-6*(rand(size(acd)-0.5))) also does the trick. Dirk From help-octave-request@bevo.che.wisc.edu Wed Jul 29 01:30:32 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id BAA15128 for help-octave-dist@bevo.che.wisc.edu; Wed, 29 Jul 1998 01:30:32 -0500 (CDT) Resent-Date: Wed, 29 Jul 1998 01:30:32 -0500 (CDT) From: Dirk Laurie Message-Id: <199807290629.IAA01904@calvyn.puk.ac.za> Subject: Re: precision (was: memory exhausted) To: help-octave@bevo.che.wisc.edu Date: Wed, 29 Jul 1998 08:29:16 +0200 (SAT) In-Reply-To: <199807290613.IAA01800@calvyn.puk.ac.za> from "Dirk Laurie" at Jul 29, 98 08:13:23 am Reply-to: wskdpl@pukrs12.puk.ac.za Content-Type: text Resent-Message-ID: <"HpGBtlzv-xL.A.qzD.IGsv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Dirk Laurie wrote: > > acd=acd.*(1+1e-6*(rand(size(acd)-0.5))) > Ugh! I meant acd=acd.*(1+1e-6*(rand(size(acd))-0.5)) Dirk From help-octave-request@bevo.che.wisc.edu Wed Jul 29 03:32:01 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id DAA13960 for help-octave-dist@bevo.che.wisc.edu; Wed, 29 Jul 1998 03:32:01 -0500 (CDT) Resent-Date: Wed, 29 Jul 1998 03:32:01 -0500 (CDT) Message-Id: <35BEDDD9.1F73D247@bmw.de> Date: Wed, 29 Jul 1998 10:31:21 +0200 From: Daniel Heiserer X-Mailer: Mozilla 4.04 [en] (X11; I; IRIX64 6.2 IP28) Mime-Version: 1.0 To: "A. Scottedward Hodel" , "help-octave@bevo.che.wisc.edu" Subject: Re: memory exhausted, precision ok References: <199807281838.NAA14238@Eng.Auburn.EDU> Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Resent-Message-ID: <"TWp6kJnpxxG.A.GPD.B4tv1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu A. Scottedward Hodel wrote: > > I ran the following code on 1 172x172 matrix: > > # truncate acd to 6 digits of precision... > for ii=1:rows(acd) > for jj = 1:columns(acd) > eval(["acd(ii,jj) = ",sprintf("%12.6E",acd(ii,jj));]); > endfor > endfor > Hi, here is my solution of your precision stuff (no ideas about the memory): ------------------ chop --------------------------- function back=precision(in,prec) % PRECISION(X,PREC) % precision chops some digits. The accuracy of % the output equals to prec. Default for prec is 6. % WARNING: works for real only % mailto: daniel.heiserer@bmw.de if (nargin <2) prec=6; disp('Warning using default precision.'); else prec=floor(prec); end disp(['Precision is ',num2str(prec),' digits.']); prec=prec-1; % take care of the signature signature=sign(in); absin=abs(in); % split into exponent and base exponent=log(absin)./log(10); digits=absin./10.^floor(exponent); % well, no lets take the precision digits=round(digits*10^prec)/10^prec; % lets merge the stuff again back=signature.*digits.*10.^floor(exponent); ---------------------- chop -------------------- -- Mit freundlichen Gruessen Daniel Heiserer ----- -------------------------------------------------------------- Daniel Heiserer, BMW AG, Knorrstrasse 147, 80788 Muenchen Abteilung EK-20 Tel.: 089-382-21187, Fax.: 089-382-42820 mailto:daniel.heiserer@bmw.de From help-octave-request@bevo.che.wisc.edu Wed Jul 29 13:53:51 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id NAA18923 for help-octave-dist@bevo.che.wisc.edu; Wed, 29 Jul 1998 13:53:51 -0500 (CDT) Resent-Date: Wed, 29 Jul 1998 13:53:51 -0500 (CDT) From: "John W. Eaton" MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Message-ID: <13759.28602.453901.931488@tillamook-sharp.bogus.domain> Date: Wed, 29 Jul 1998 13:53:46 -0500 (CDT) To: "A. Scottedward Hodel" Cc: help-octave@bevo.che.wisc.edu Subject: memory exhausted In-Reply-To: <199807281838.NAA14238@Eng.Auburn.EDU> References: <199807281838.NAA14238@Eng.Auburn.EDU> X-Mailer: VM 6.46 under Emacs 19.34.1 Resent-Message-ID: <"qc-YY_OFzeK.A.YHH._-2v1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu On 28-Jul-1998, A. Scottedward Hodel wrote: | I ran the following code on 1 172x172 matrix: | | # truncate acd to 6 digits of precision... | for ii=1:rows(acd) | for jj = 1:columns(acd) | eval(["acd(ii,jj) = ",sprintf("%12.6E",acd(ii,jj));]); | endfor | endfor | | After some time into the loop, Octave repeatedly | displays an error message: | | "memory exhausted --- trying to return to prompt" As other people have noted, there may be more efficient ways to do the above operation. My comments are only about the memory problems. The trouble is that for each element of the matrix acd, the statement eval(["acd(ii,jj) = ",sprintf("%12.6E",acd(ii,jj));]); will cause Octave to evaluate something like acd(ii,jj) = 5.746209E-01 which also results in Octave printing the entire matrix. Unless paging is turned off or immediate paging is enabled, Octave will buffer all the output before displaying it. Using the default output format for a 172x172 matrix requires about 350 kilobytes of storage. Since the entire matrix is displayed for each assignment (172x172 of them), the total storage required for buffering is approximately 10354400 kilobytes, or about 9.87 gigabytes! There are a couple of ways to keep Octave's eval statement from printing the result. One is to end the string to be evaluated with a semicolon: eval(["acd(ii,jj) = ",sprintf("%12.6E;",acd(ii,jj));]); Another is to set the built-in variable `default_eval_print_flag' to 0. Also, I would write the eval statement this way: eval (sprintf ("acd(ii,jj) = %12.6E;", acd(ii,jj))); which should be slightly faster as it eliminates the relatively expensive string concatenation inside the square brackets. jwe From help-octave-request@bevo.che.wisc.edu Wed Jul 29 15:53:58 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id PAA30318 for help-octave-dist@bevo.che.wisc.edu; Wed, 29 Jul 1998 15:53:58 -0500 (CDT) Resent-Date: Wed, 29 Jul 1998 15:53:58 -0500 (CDT) Date: Wed, 29 Jul 1998 14:53:32 -0600 (MDT) From: Herng-Jeng Jou To: help-octave@bevo.che.wisc.edu Subject: Dynamical Linking on RedHat 5 Linux Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"SXgAEUT2sTJ.A.J0E.mv4v1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Hello, I use 2.0.13 (compiled from source gcc 2.7.2.3, g77 egcs-1.02) on Linux RedHat 5.0. When I tested the hello.cc example, it compiled with mkoctfile ok. But when loading it within octave, I got error : octave: error in loading shared libraries /export/home2/hjjou/code/expfm/hello.oct: undefined symbol: stream__19octave_pager_stream octave was compiled with -DWITH_DYNAMIC_LINKING=1 Is it a installation problem? Thanks, Jou From help-octave-request@bevo.che.wisc.edu Thu Jul 30 00:14:23 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id AAA18158 for help-octave-dist@bevo.che.wisc.edu; Thu, 30 Jul 1998 00:14:23 -0500 (CDT) Resent-Date: Thu, 30 Jul 1998 00:14:23 -0500 (CDT) From: "John W. Eaton" MIME-Version: 1.0 Content-Type: text/plain; charset=us-ascii Content-Transfer-Encoding: 7bit Message-ID: <13760.301.38389.229331@tillamook-sharp.bogus.domain> Date: Thu, 30 Jul 1998 00:14:21 -0500 (CDT) To: "Wonkoo Kim" Cc: "Octave Users" Subject: diary In-Reply-To: <199807260119.VAA002.86@wkim> References: <199807260119.VAA002.86@wkim> X-Mailer: VM 6.46 under Emacs 19.34.1 Resent-Message-ID: <"i0tQ5tuOC-N.A.qEG.vEAw1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu On 25-Jul-1998, Wonkoo Kim wrote: | I think the problem is solved. | | I had put | diary "file.log" | at the beginning of a .m file, and | diary off | at the end of the .m file. However, it looks like I have to | remove "diary off" from the .m file, and I have to issue it | directly from Octave session. I don't know why, but it worked | that way. I agree that this is a bug, but I see no easy way to fix it. The problem is that nothing is written to the diary file until output is flushed to the pager, and that normally only happens just before a prompt is printed. jwe From help-octave-request@bevo.che.wisc.edu Thu Jul 30 02:03:27 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id CAA00568 for help-octave-dist@bevo.che.wisc.edu; Thu, 30 Jul 1998 02:03:27 -0500 (CDT) Resent-Date: Thu, 30 Jul 1998 02:03:27 -0500 (CDT) Message-Id: <199807300703.DAA000.78@wkim> MIME-Version: 1.0 Date: Thu, 30 Jul 98 02:23:17 -0400 From: "Wonkoo Kim" To: "John W. Eaton" Reply-To: "Wonkoo Kim" Cc: "Octave Users" Subject: Re: diary X-Mailer: Ultimedia Mail/2 Lite, IBM T. J. Watson Research Center Content-ID: <76_74_4_901779798> Content-Type: text/plain; charset="US-ASCII" Content-Transfer-Encoding: 7bit Resent-Message-ID: <"DDnfiU9W0WB.A.hyB._qBw1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu > Date: Thu, 30 Jul 98 00:14:21 -0500 > From: "John W. Eaton" > Subject: diary > > On 25-Jul-1998, Wonkoo Kim wrote: > > | I think the problem is solved. > | > | I had put > | diary "file.log" > | at the beginning of a .m file, and > | diary off > | at the end of the .m file. However, it looks like I have to > | remove "diary off" from the .m file, and I have to issue it > | directly from Octave session. I don't know why, but it worked > | that way. > > I agree that this is a bug, but I see no easy way to fix it. > > The problem is that nothing is written to the diary file until output > is flushed to the pager, and that normally only happens just before a > prompt is printed. > > jwe Is ``diary off'' in .m file already executed when pager gets output? It looks like Octave assumes ``diary off'' at the end of .m file was issued at the beginning of paged output, and so I got nothing in diary file for the paged output. (Thus I had to remove ``diary off'' from .m file.) Perhaps Octave needs a special treat for ``diary off'' in .m file when output is paged. Thanks for your answer. //-------------------------------------------------------------------- // Wonkoo Kim (wkim+@pitt.edu) From help-octave-request@bevo.che.wisc.edu Thu Jul 30 08:56:35 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id IAA06110 for help-octave-dist@bevo.che.wisc.edu; Thu, 30 Jul 1998 08:56:35 -0500 (CDT) Resent-Date: Thu, 30 Jul 1998 08:56:35 -0500 (CDT) Date: Thu, 30 Jul 1998 07:56:13 -0600 (MDT) From: Herng-Jeng Jou To: help-octave@bevo.che.wisc.edu Subject: Dynamical Linking on RedHat 5.0 Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"xGYtV26gNyP.A.o1B.TuHw1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Hi, I have problem using dynamical linking on RedHat Linux 5.0. I tried hello.cc example (no problem compiled with mkoctfile). When loaded within octave, here is the error message I got: octave: error in loading shared libraries /export/home2/hjjou/code/expfm/hello.oct: undefined symbol: stream__19octave_pager_stream I used octave 2.0.13 compiled with egcs-2.90.27 980315 (egcs-1.0.2 release) on RedHat 5.0 Linux system. octave_config_info does show -DWITH_DYNAMIC_LINNKING=1 Is there installation problem. Thanks, Jou From help-octave-request@bevo.che.wisc.edu Thu Jul 30 15:15:24 1998 Received: (from slist@localhost) by bevo.che.wisc.edu (8.8.8/8.8.8) id PAA11733 for help-octave-dist@bevo.che.wisc.edu; Thu, 30 Jul 1998 15:15:24 -0500 (CDT) Resent-Date: Thu, 30 Jul 1998 15:15:24 -0500 (CDT) >Received: from localhost by csnio.ren.nic.in with smtp id m0z1vze-0001Z8C (Debian Smail-3.2 1996-Jul-4 #2); Thu, 30 Jul 1998 22:25:54 +0530 (IST) Date: Thu, 30 Jul 1998 22:25:54 +0530 (IST) From: "Mr.Jaydeep Pattanaik IODC NIO" To: Help Octave Subject: netcdf In-Reply-To: <13758.48604.974930.604960@tillamook-sharp.bogus.domain> Message-ID: MIME-Version: 1.0 Content-Type: TEXT/PLAIN; charset=US-ASCII Resent-Message-ID: <"dfblRjoMq1N.A.eTG.cRNw1"@bevo.che.wisc.edu> Resent-From: help-octave@bevo.che.wisc.edu X-Mailing-List: X-Loop: help-octave@bevo.che.wisc.edu Precedence: list Resent-Sender: help-octave-request@bevo.che.wisc.edu Sender: help-octave-request@bevo.che.wisc.edu Hello, Has anybody worked with NETCDF (NETwork Common Data Format) data in octave. If yes, please provide me some insights. Thanks in advance. Jaydeep ---------------------------------------------------------------------- On Wed, 29 Jul 1998, John W. Eaton wrote: > On 28-Jul-1998, Mr.Jaydeep Pattanaik IODC NIO wrote: > > | I am resending this mail as I did not receive any reply. > > Your first two messages were sent to help-octave-request, which is > only read by me. Your last message was sent to help-octave (the whole > list) but it was empty. Perhaps you should resend it to the list to > find out if anyone has done any work with netcdf (I haven't). As for > the second question, try the polyfit function. > > jwe > > | ---------- Forwarded message ---------- > | Date: Tue, 21 Jul 1998 10:23:53 +0530 (IST) > | From: "Mr.Jaydeep Pattanaik IODC NIO" > | To: help-octave-request@bevo.che.wisc.edu > | Subject: netcdf > | > | Dear all, > | > | I am a first timer in Octave. Want to know whether > | > | 1. I can directly use netcdf data in octave. > | 2. Fit a second order polynomial to my data. > | > | Any hints will be greately appriciated. > | > | > | Jaydeep > | > | ---------------------------------------------------------------------- > | Jaydeep Pattanaik, |Phone - 091-0832-229900 Ext-4241 > | Indian Oceanographic Data Centre, |Fax - 091-0832-223340 > | National Institute of Oceanography, | 091-0832-229102 > | Dona Paula, Goa - 403 004, INDIA |e-mail: jaydeep@csnio.ren.nic.in > | | jaydeep@darya.nio.org > | ---------------------------------------------------------------------- >