Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge -> triangular example

2011-11-12 Thread Michele Martone
On 2012@18:43, Carlo de Falco wrote: > On 12 Nov 2011, at 14:17, Michele Martone wrote: > > On 2008@15:41, Michele Martone wrote: > >> On 2008@15:12, Carlo de Falco wrote: > >>> ... > >>> There is an example of the correct way to create a new

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge

2011-11-12 Thread Michele Martone
On 2012@18:25, Jordi GutiƩrrez Hermoso wrote: > On 12 November 2011 18:18, Carlo de Falco wrote: > > On 12 Nov 2011, at 16:04, c. wrote: > >> On 12 Nov 2011, at 12:17, Alexander Barth wrote: > > > > It looks like it's working: > > > >>

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge

2011-11-12 Thread Jordi GutiƩrrez Hermoso
On 12 November 2011 18:18, Carlo de Falco wrote: > > On 12 Nov 2011, at 16:04, c. wrote: > >> >> On 12 Nov 2011, at 12:17, Alexander Barth wrote: >> >>> Dear Filippone, >>> >>> Also an iterative parallel solver would be very useful for me. >>> Currently, I'm using either a direct solver (the \ ope

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge

2011-11-12 Thread Carlo de Falco
On 12 Nov 2011, at 16:04, c. wrote: > > On 12 Nov 2011, at 12:17, Alexander Barth wrote: > >> Dear Filippone, >> >> Also an iterative parallel solver would be very useful for me. >> Currently, I'm using either a direct solver (the \ operator) or the >> conjugate gradient alogorithm. >> I tried

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge -> triangular example

2011-11-12 Thread c.
On 12 Nov 2011, at 18:43, Carlo de Falco wrote: > 7) with a few changes I can get the sparsersb to install cleanly, but I get a > bus-error when I try to use it: > > a = > > -1.214470 1.014121 -0.565035 > 0.372437 0.935208 -0.849173 > -0.260013 0.085202 -0.362010 > >>> ss = spars

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge -> triangular example

2011-11-12 Thread Carlo de Falco
On 12 Nov 2011, at 14:17, Michele Martone wrote: > On 2008@15:41, Michele Martone wrote: >> On 2008@15:12, Carlo de Falco wrote: >>> ... >>> There is an example of the correct way to create a new octave matrix class >>> in the package "triangular": >>> http://octave.svn.sourceforge.net/

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge

2011-11-12 Thread Salvatore Filippone
Il giorno sab, 12/11/2011 alle 12.17 +0100, Alexander Barth ha scritto: > Dear Filippone, > > Also an iterative parallel solver would be very useful for me. > Currently, I'm using either a direct solver (the \ operator) or the > conjugate gradient alogorithm. > I tried to find our more information

Re: [OctDev] Stable poles handed to "filter" can result in unstable filter

2011-11-12 Thread Alberto Cuadrado
I think the transfer function you provide is right (see help on "filter" function). In the vector "a" returned by butter.m, a(1) corresponds to 1 and a(N+1) corresponds to a_N*z^-N, being "N" the filter order. However, the poles of the filter are roots(a), not 1./roots(a), because if you multiply b

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge -> triangular example

2011-11-12 Thread Carlo de Falco
On 12 Nov 2011, at 14:17, Michele Martone wrote: > On 2008@15:41, Michele Martone wrote: >> On 2008@15:12, Carlo de Falco wrote: >>> ... >>> There is an example of the correct way to create a new octave matrix class >>> in the package "triangular": >>> http://octave.svn.sourceforge.net/

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge

2011-11-12 Thread c.
On 12 Nov 2011, at 12:17, Alexander Barth wrote: > Dear Filippone, > > Also an iterative parallel solver would be very useful for me. > Currently, I'm using either a direct solver (the \ operator) or the > conjugate gradient alogorithm. > I tried to find our more information about PSBLAS but it

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge

2011-11-12 Thread Salvatore Filippone
As the originator of PSBLAS, I should perhaps point out that it contains iterative solvers, not the direct ones that are implied by the \ operator. Our intention, as far as Octave is concerned, is to make available the Krylov solvers we have, coupled with the multilevel preconditioners of MLD2P4; b

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge -> triangular example

2011-11-12 Thread Michele Martone
On 2008@15:41, Michele Martone wrote: > On 2008@15:12, Carlo de Falco wrote: > > ... > > There is an example of the correct way to create a new octave matrix class > > in the package "triangular": > > http://octave.svn.sourceforge.net/viewvc/octave/trunk/octave-forge/extra/triangular/ > O

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge

2011-11-12 Thread Michele Martone
On 2012@11:44, Alexander Barth wrote: > On Fri, Nov 11, 2011 at 11:48 PM, Michele Martone > wrote: > > On 2011@22:29, Alexander Barth wrote: > ... > I have had some issues/doubts when I compiled librsb et sparsersb. > When I configured librsb (version 0.0.1445)

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge

2011-11-12 Thread Alexander Barth
Dear Filippone, Also an iterative parallel solver would be very useful for me. Currently, I'm using either a direct solver (the \ operator) or the conjugate gradient alogorithm. I tried to find our more information about PSBLAS but it seems that the web-page has currently some problems (time-out).

Re: [OctDev] "librsb"+"sparsersb" packages proposal for octave-forge

2011-11-12 Thread Alexander Barth
On Fri, Nov 11, 2011 at 11:48 PM, Michele Martone wrote: > On 2011@22:29, Alexander Barth wrote: >> > >> > hello, octave-forge :) >> > created main/sparsersb, in order to host files for the "sparsersb" >> > oct-file package; that is, an Octave function to have an additional Octave >> > type