gmane.comp.gnu.octave.general

Optimising the Keen-Minsky Model

richard — 2012-05-25T19:18:35

Hi c. 
please find attached the Keen-Minsky Model ver02 written in
Octave. I trust you will find that the earlier cut down 
code simplifies the task of optimising the program without 
the complication of the full model 

HTH

Richard A Lough

#! ~/octave/RD/Keen-Minsky/KeenMinsky02.m
# An Octave program to verify model of firms and bank accounts
# including use of the the Phillips curve for non linearities, 
# as published by Prof Keen, his Model, MonetaryMinsky03  
# Octave Program author Richard A Lough. Published under the GPL License
# ver .02 published 20120525
######################################################################## 
# Set up the X matrix, initialise with ones
X= ones(12001,16); # Change this to ones(7202,16) to replicate Keen's model 
Y=zeros(12001,10); # Added for graphic output
for j = 1: 10      # this loop added to Keen's model for development work
Minsky = 12/(2+j); # Keen's model had a fixed value of 1.0
X(1,1)= 100; X(1,2)=12; X(1,3)=70; X(1,4)=5; X(1,5)= 13; X(1,6)= 1; X(1,7)=9

Re: mkoctfile and mex broken on Mac 10.6 Snow Leopard (with patchof a sort)

Sciencegal — 2012-05-25T19:16:23

RickG wrote

I know it's been awhile since you posted this, but maybe someone will pick
up on it. I included these statements in the mkoctfile_3.2.3 as you
suggested, but then when starting Octave I get the message:

/Applications/Octave.app/Contents/Resources/bin/mkoctfile-3.2.3: line 29:
i386: command not found

Any idea what's wrong here?

--
View this message in context: http://octave.1599824.n4.nabble.com/mkoctfile-and-mex-broken-on-Mac-10-6-Snow-Leopard-with-patch-of-a-sort-tp1636977p4630161.html
Sent from the Octave - General mailing list archive at Nabble.com.
_______________________________________________
Help-octave mailing list
Help-octave< at >octave.org
https://mailman.cae.wisc.edu/listinfo/help-octave

Re: Lomb-Scargle variant

Benjamin Lewis — 2012-05-25T17:24:56

Hi Mark,
This looks great! I'll definitely add it to the wishlist, probably near the top. I want to at least get one variant written, but this looks like a useful article overall. Thanks!
Ben
Le 2012-05-24 à 09:06, Mark Lytle a écrit :


_______________________________________________
Help-octave mailing list
Help-octave< at >octave.org
https://mailman.cae.wisc.edu/listinfo/help-octave

Re: Lomb-Scargle variant

Rhys Ulerich — 2012-05-25T17:19:36

Hi Mark,


Any chance you could provide a few uncommon phrases from that print
out?  I'm having trouble Googling for it from your hints and would
appreciate "some curious words" on which to search.

Thanks,
Rhys

Re: Are package folder supported?

Michael Goffioul — 2012-05-25T14:26:29

On Fri, May 25, 2012 at 9:39 AM, Jakob Lombacher
<lombacher< at >googlemail.com>wrote:


No, they're not supported.

Michael.
_______________________________________________
Help-octave mailing list
Help-octave< at >octave.org
https://mailman.cae.wisc.edu/listinfo/help-octave

Are package folder supported?

Jakob Lombacher — 2012-05-25T08:39:53

Hi,

I'd like to run some Matlab code with package folders (folders like
'+packageName')

http://www.mathworks.de/help/techdoc/matlab_oop/brfynt_-1.html

It seems that this doesn't work in octave 3.2.4.

Is there a workarround or description about this issue, I couldn't
find anything about this via google?

Thanks

Jakob

Re: Help in porting a toolbox

Juan Pablo Carbajal — 2012-05-25T07:15:01

On Fri, May 25, 2012 at 7:56 AM, Mathieu Dubois
<mathieu.dubois< at >ensta-paristech.fr> wrote:

Hi,

It looks like the problem is that those packages where built long ago
and now you have newer versions of the libraries. Since those aren't
good packages they can't check. I suggest you give it a try compiling
it yourself.
follow the instructions here
http://www.gnu.org/software/octave/doc/interpreter/Installation.html#Installation

If you find any errors, check here some Ubuntu specific explanations.
ther eis no hack involved just calling 'sudo apt-get install', quite
easy.
http://octave.org/wiki/index.php?title=Octave_for_GNU_Linux:_Binary_Octave_packages_for_GNU_Linux#Ubuntu

Cheers,

Re: Help in porting a toolbox

Mathieu Dubois — 2012-05-25T05:56:49


----- Mail original -----
De: "Jordi Gutiérrez Hermoso" <jordigh< at >octave.org>
À: "Octave-help" <help-octave< at >octave.org>
Envoyé: Jeudi 24 Mai 2012 23:17:57
Objet: Re: Help in porting a toolbox

On 24 May 2012 16:35, Mathieu Dubois <mathieu.dubois< at >ensta-paristech.fr> wrote:




I have installed octave 3.6.1 for Ubuntu 11.04 (in fact I'm under 11.04 - I'm in the process of upgrading to 12.04 :) but it fails to run:
$ octave
octave: /usr/lib/x86_64-linux-gnu/libstdc++.so.6: version `GLIBCXX_3.4.15' not found (required by /usr/local/lib/octave/3.6.1/liboctinterp.so.1)
octave: /usr/lib/x86_64-linux-gnu/libstdc++.so.6: version `GLIBCXX_3.4.15' not found (required by /usr/local/lib/octave/3.6.1/liboctave.so.1)
octave: /usr/lib/x86_64-linux-gnu/libgfortran.so.3: version `GFORTRAN_1.4' not found (required by /usr/local/lib/octave/3.6.1/libcruft.so.1)

Any solution?



Maybe I have not understood your message. You think it's a bug because we are supposed to have the same behavior as Matlba, right?
In that case, an

Re: Help in porting a toolbox

Carnë Draug — 2012-05-25T00:40:23

As Jordi has already mentioned, you should really report bugs when you
find them. The community needs to help if octave is to become even
better. That said, this type of code is necessary while full
compatibility does not exist. What you want to do is explained on the
Octave wiki at http://octave.org/wiki/index.php?title=Compatibility

Carnë

Re: GSoC: Student blogs

richard — 2012-05-24T23:50:18

On Thu, May 24, 2012 at 01:41:33PM +0200, c. wrote:  

I suspected we might arrive here :-)

I will ask and I will think about a stripped down version. 

For information, the original development was in MathCad, and published here: 
http://www.debtdeflation.com/blogs/qed/

Note that the comments state the code is not Open Source.   

There are some other, related, sources on the internet. 

If you are content with a _very_ stripped down version add these : 
initial conditions: 
x4=0.15; y4=0.33; s4=5; m4=0;
Algebra: 
h=i;
g = (y4-m4)*exp((X(i,7)-x4)*s4/(y4-m4)) +m4; # phillips curve calculation
g=g+Variable; 

# This is perhaps one twentieth of the calculations carried out in the inner loop. 
# I'd suppose if you did this 20 times it might approximate to the 
# CPU load for the model. 


Sometimes tools such as, quadratic optimisation, or least squared errors can 
be used to quickly maximise (or minimise) some feature of your model without 
much human intervention. For example, minimising the time taken for

Re: Help in porting a toolbox

Jordi Gutiérrez Hermoso — 2012-05-24T21:17:57


Can you try a newer Octave version?

   http://www.octave.org/wiki/index.php?title=Octave_for_GNU_Linux:_Binary_Octave_packages_for_GNU_Linux#Unofficial_binaries


No, yes, *please* report bugs if you find them. We can't make Octave
better if we don't even know what's wrong with it.

- Jordi G. H.
_______________________________________________
Help-octave mailing list
Help-octave< at >octave.org
https://mailman.cae.wisc.edu/listinfo/help-octave

Re: Extracting matrix rows that meet a condition

Ismael Núñez-Riboni — 2012-05-24T20:55:43

Maybe some values are identically zero or identically one? You should 
try >= or <= in one of your conditions...

Re: Help in porting a toolbox

Mathieu Dubois — 2012-05-24T20:35:01


----- Mail original -----
De: "Jordi Gutiérrez Hermoso" <jordigh< at >octave.org>
À: "Mathieu Dubois" <mathieu.dubois< at >ensta-paristech.fr>
Cc: help-octave< at >octave.org
Envoyé: Jeudi 24 Mai 2012 22:10:09
Objet: Re: Help in porting a toolbox

On 24 May 2012 15:43, Mathieu Dubois <mathieu.dubois< at >ensta-paristech.fr> wrote:

[snip]


I use Octave 3.2.4 (under Ubuntu 11.10).




I don't do that! I was wondering if it would be a way to work aroud the bug.

Mathieu

Re: Help in porting a toolbox

Jordi Gutiérrez Hermoso — 2012-05-24T20:10:09

[snip]

What Octave version? I just tried this example but got n = 6, not 0.



If you need to do this, then you have almost surely found an Octave
bug. Please report it:

    http://www.gnu.org/software/octave/bugs.html

However, if you need to work around the bug, you can check

    if (exist OCTAVE_VERSION)

- Jordi G. H.
_______________________________________________
Help-octave mailing list
Help-octave< at >octave.org
https://mailman.cae.wisc.edu/listinfo/help-octave

Help in porting a toolbox

Mathieu Dubois — 2012-05-24T19:43:57

Hello,

I'm a new Octave user (I have some experience using Matlab).

I wanted to port the SOM toolbox (see: http://www.cis.hut.fi/somtoolbox/) to Octave.
The code is a bit old but it works under Matlab 2008a (I'm stick to this version which is one of the reason why I want to leave Matlab).

I have some problems with code involving (s)scanf. I have read that there are a lot of differences between Octave and Matlab in this area.
The problematic line can be found in som_read_data.m, line 225:
  [s, c, e, n] = sscanf(li, '%s%[^ \t]');
The goal of the line is to read the header of a file which contains the data label.
The variable li (the line) contains "SepalL SepalW PetalL PetalW".
Under Matlab, s contains "SepalL" (the first label) and the variable n (which causes the error) contains 8.
Under Octave, s contains the correct value but n is 0.
This value creates an error in line 226 (because n is used to acces elements of li).

I have used the "C" form which is closer to Matlab if I understand correctly:
[s, c,

Re: Extracting matrix rows that meet a condition

Knof — 2012-05-24T20:01:47

Something is still broken somewhere as 426 rows match the first condition and 619 rows match the second condition.
426 + 619 = 1045 and my matrix only has 1008 rows.

--
View this message in context: http://octave.1599824.n4.nabble.com/Extracting-matrix-rows-that-meet-a-condition-tp4630112p4630135.html
Sent from the Octave - General mailing list archive at Nabble.com._______________________________________________
Help-octave mailing list
Help-octave< at >octave.org
https://mailman.cae.wisc.edu/listinfo/help-octave

Re: Problem installing packages

Sergei Steshenko — 2012-05-24T17:45:21




----- Original Message -----
[snip]


Oh, really ? What a long way to say "buggy".

And what about people who build Octave and packages from source ?
And what about GPL requiring both the source and build instructions from that source to be available ?
AFAIR, GPL explicitly forbids obscuring ways of building from source.

...

Was the OP trying to install packages globally ?

I judge from '/usr/local' in:

"
couldn't create installation directory
/usr/share/octave/packages/combinatorics-1.0.9 : Permission denied
error: called from `pkg>copy_files' in file
/usr/share/octave/3.6.1/m/pkg/pkg.m near line 1569, column 7
error: called from:
error:   /usr/share/octave/3.6.1/m/pkg/pkg.m at line 876, column 5
error:   /usr/share/octave/3.6.1/m/pkg/pkg.m at line 383, column 9
".

Regards,
  Sergei.

Re: Extracting matrix rows that meet a condition

Ismael Diego Nunez-Riboni — 2012-05-24T16:01:45

Then forget about the function "find", what you want to do is:

 ii = (x > 0 & x < 1);

The rows you want to save are:

x((sum(ii') == 3)',:)

They are only 619...

Re: Extracting matrix rows that meet a condition

Knof — 2012-05-24T15:36:13

The problem is I want to spit the matrix into two files. One file where all the rows that contain at least one value that meets the condition
(x < 0 | x > 1) are written to and the other file where all the rows where all three values match the condition (x > 0 & x < 1) are saved to it.
So the second file must contain all the rows that the first file does not.
I am sorry if I have not made that more clear before.

Knof



--
View this message in context: http://octave.1599824.n4.nabble.com/Extracting-matrix-rows-that-meet-a-condition-tp4630112p4630131.html
Sent from the Octave - General mailing list archive at Nabble.com._______________________________________________
Help-octave mailing list
Help-octave< at >octave.org
https://mailman.cae.wisc.edu/listinfo/help-octave

Re: Extracting matrix rows that meet a condition

Ismael Diego Nunez-Riboni — 2012-05-24T15:24:05

-------- Original-Nachricht --------


It works in my case, I get 1005 out of 1008 elements, i.e., apparently almost all rows have one element meeting the condition x > 0 & x < 1... I see no contradiction... Here it is how I see it (using of course the matrix you attached):

[i,j] = find(x < 0 | x > 1);

octave:41> size(x)
ans =

   1008      3

octave:42> size(i)
ans =

   426     1

octave:43> size(j)
ans =

   426     1


octave:45> [k,l] = find(x > 0 & x < 1);

octave:46> size(k)
ans =

   2592      1  % Repeated indices pointing to same row

octave:48> size(unique(k))
ans =

   1005      1  % Unique number of indices < total row no. of x

I fail to see the problem... Almost every single row in x have a value meeting
x > 0 & x < 1...

Re: Problem installing packages

Jordi Gutiérrez Hermoso — 2012-05-24T15:11:36

There is no GUI yet. If you mean QtOctave, it's abandoned and
scheduled to be deleted from the repositories.


You should not be using pkg install to install OF packages in Debian.
They are already packaged for Debian with apt. So e.g. instead of
doing "pkg install signal" from Octave, you should do "aptitude
install octave-signal" from bash.

In fact, I wish we could stop telling people to use pkg install,
because it's a very brittle way of getting packages to people.

HTH,
- Jordi G. H.