Dev March 2009

dev@dynare.org

6 participants
15 discussions

model derivatives with respect to parameters
by Michel Juillard 19 Aug '09

19 Aug '09

Dear Nikolay, we have started to work on producing the derivatives with respect to parameters. We are planing to create a new function, that will return a 3-dimensional array, by equation, variable and parameter. You could then select one parameter and squeeze the array in the parameter dimension to obtain a matrix of Jacobian derivatives with respect to a given parameter. Do you see any problem with this implementation? We should be able to provide in about two weeks. All the best, Michel

4 6

dlls for Matlab 7.8
by Michel Juillard 30 Mar '09

30 Mar '09

I fixed build_matlab.m and mjdgges.c for compatibility issue with Matlab 7.8 When we start releasing binaries for Windows 64bit, we will need to rethink the architecture of directories and distribution mode of the DLLs Best Michel

1 0

forum dynare v4
by Stéphane Adjemian 28 Mar '09

28 Mar '09

Hi Michel, How can we access to the (closed) forum of Dynare version 4? à+ Stéphane.

2 2

Matlab r2009a mexfiles BLAS LAPACK
by Michel Juillard 28 Mar '09

28 Mar '09

There is indeed a change in arguments to calls to Blas or Lapack functions on 64 bit platforms. See http://www.mathworks.com/access/helpdesk/help/techdoc/index.html?/access/he… "MEX-Files Calling BLAS or LAPACK Functions Must Be Updated On 64-Bit Platforms" I will change the mex sources expanding the mechanism we put in place for the pointers across different Matlab versions. Hopefully over the week-end. Best Michel

1 0

r2009a
by Michel Juillard 28 Mar '09

28 Mar '09

I have installed Matlab r2009a on my Linux Debian and the old DLLs work with it. I ran build_matlab.m and I don't have the link problem reported on the forum, but there are warnings on at the compilation of mjdgges.c and the DLL doesn't work. I will investigate the issue and let you know. Best Michel

1 0

Re: [DynareDev] AIM in Matlab (and DYNARE question)
by Michel Juillard 27 Mar '09

27 Mar '09

1) I can't either create pages on DynareWiki. The software was recently updated and this must have created the problem. 2) It makes good sense that problems with the solution method will trash estimation results. However, I prefer to test with stoch_simul, because the result is simpler and is supposed to be reproducible rather than with estimation that is a much longer chain of operations besides solution and may depend of random draws in the optimization step. Thanks again for your help on that. Best Michel G. Perendia wrote: > Dear Michel > > 1) I do not seem to be able to edit/create new Wiki page: when I add a new > title to my page (e.g MatlabEigsBug) and then click link to it after saving > my page, I get an unfamiliar page saying: > > This page does not exist yet. You can create a new empty page, or use one of > the page templates. > Action(edit,Create new empty page) > > Or choose a page template to create an already-formatted page: > Before creating the page, please check if a similar page already > exists. Here are some existing pages with similar names: > > EditTemplates > LikePages((none)) > > > > > but clicking on any of the given options makes no effect. > > > 2) Running estimation was useful after all: The eigs() problem was causing > very many AIM errors of type 4 - too few big roots.- to be reported during > the ML and MCMC estimation (see the ..AIM_truncated_run log I sent this > morning) but none of type 4 when I rerun the same estimation after the > change - only 2 of type 3 - too many big roots. Also, it affected the > results of the estimation : the likelihood was rather lower before than > that after the modification and the estimated parameters different. See new > log enclosed. > > Best regards > > George > Mob. +44(0)7951415480 > artilogica(a)btconnect.com > > ----- Original Message ----- > From: "Michel Juillard" <michel.juillard(a)ens.fr> > To: "G. Perendia" <george(a)perendia.orangehome.co.uk> > Sent: Thursday, March 26, 2009 8:03 PM > Subject: Re: [DynareDev] AIM in Matlab (and DYNARE question) > > > >> Thank you very much George. I can't access the bug report, but it's OK. >> >> Could you please add a page to DynareWiki describing the problem a) in >> Matlab b) its consequence for AIM >> >> For the other test cases that you prepare, could you please remove the >> estimation part and do the same comparison as I did with and without >> AIM? I'm curious whether the bug is frequent or if we found a crazy case >> by chance. >> >> All the best, >> >> Michel >> >> >> G. Perendia wrote: >> >>> Hi Michel >>> >>> FYI: Mathworks eigs() bug report record (though I'm not sure if you can >>> access it this way): >>> >>> >>> > http://www.mathworks.com/support/service_requests/Service_Request_Detail.do… > >>> Best regards >>> >>> George >>> ----- Original Message ----- >>> From: "Michel Juillard" <michel.juillard(a)ens.fr> >>> To: <Gary.Anderson(a)frb.gov>; "List for Dynare developers" >>> > <dev(a)dynare.org> > >>> Sent: Thursday, March 26, 2009 5:27 PM >>> Subject: Re: [DynareDev] AIM in Matlab (and DYNARE question) >>> >>> >>> >>> >>>> Hi Gary, >>>> >>>> I will keep you posted on Mathworks answer concerning eigs(). >>>> >>>> Cocnerning the tests, I have been accumulating a series of tests over >>>> the years, that I'm planing to commit to SVN in the coming weeks. >>>> Why don't you start by seeing what I already have and then we discuss a >>>> strategy to be more systematic. >>>> >>>> All the best, >>>> >>>> Michel >>>> >>>> Gary.Anderson(a)frb.gov wrote: >>>> >>>> >>>>> michel, >>>>> that is odd. when i run eigs, i get one of two sets of 10 >>>>> > eigenvalues. > >>> i >>> >>> >>>>> didn't investigate the invariant space vectors. >>>>> >>>>> >>>>> sometime ago i was supposed to develop a suite of tests for dynare. >>>>> >>>>> >>> since >>> >>> >>>>> it has been so long i thought i should touch base with you >>>>> to see if you think it would still be useful, and if so, what areas i >>>>> >>>>> >>> might >>> >>> >>>>> focus my energies on. >>>>> >>>>> >>>>> gary >>>>> >>>>> >>>>> >>>>> >>>>> >>>>> Michel Juillard >>>>> <michel.juillard@ >>>>> ens.fr> >>>>> >>>>> >>> To >>> >>> >>>>> Gary.Anderson(a)frb.gov, "G. >>>>> 03/26/2009 12:29 Perendia" >>>>> PM >>>>> >>>>> >>> <george(a)perendia.orangehome.co.uk>, >>> >>> >>>>> giovanni.lombardo(a)ecb.int, >>>>> Daniela.Bragoli(a)ecb.int >>>>> >>>>> >>>>> >>> cc >>> >>> >>> Subject >>> >>> >>>>> AIM in Matlab >>>>> >>>>> >>>>> >>>>> >>>>> >>>>> >>>>> >>>>> >>>>> >>>>> >>>>> Hi Gary, >>>>> >>>>> in experimenting with AIM inside Dynare, George came accross this very >>>>> weir behaviour of Matlab eigs() function. >>>>> >>>>> For a rather large matrix (101x101) >>>>> >>>>> [w,d]=eigs(a',10,'LM',opts);sort(diag(10)) >>>>> with opts containing disp:0, >>>>> >>>>> returns different eigenvalues when called repeatedly without changing >>>>> anything. >>>>> >>>>> I'm attaching a.mat that contains the a matrix (and opts, among other >>>>> things) >>>>> >>>>> For now, we are disabling the using of eigs() in Dynare and forcing >>>>> > the > >>>>> use of eig() >>>>> >>>>> We are going also to report the problem to Matlab. >>>>> >>>>> All the best, >>>>> >>>>> Michel >>>>> >>>>> (See attached file: a.mat) >>>>> >>>>> >>>> _______________________________________________ >>>> Dev mailing list >>>> Dev(a)dynare.org >>>> http://www.dynare.org/cgi-bin/mailman/listinfo/dev >>>> >>>> >>> >>> >>> >> > > > >

2 1

Re: [DynareDev] AIM in Matlab (and DYNARE question)
by Michel Juillard 26 Mar '09

26 Mar '09

Hi Gary, I will keep you posted on Mathworks answer concerning eigs(). Cocnerning the tests, I have been accumulating a series of tests over the years, that I'm planing to commit to SVN in the coming weeks. Why don't you start by seeing what I already have and then we discuss a strategy to be more systematic. All the best, Michel Gary.Anderson(a)frb.gov wrote: > michel, > that is odd. when i run eigs, i get one of two sets of 10 eigenvalues. i > didn't investigate the invariant space vectors. > > > sometime ago i was supposed to develop a suite of tests for dynare. since > it has been so long i thought i should touch base with you > to see if you think it would still be useful, and if so, what areas i might > focus my energies on. > > > gary > > > > > > Michel Juillard > <michel.juillard@ > ens.fr> To > Gary.Anderson(a)frb.gov, "G. > 03/26/2009 12:29 Perendia" > PM <george(a)perendia.orangehome.co.uk>, > giovanni.lombardo(a)ecb.int, > Daniela.Bragoli(a)ecb.int > cc > > Subject > AIM in Matlab > > > > > > > > > > > Hi Gary, > > in experimenting with AIM inside Dynare, George came accross this very > weir behaviour of Matlab eigs() function. > > For a rather large matrix (101x101) > > [w,d]=eigs(a',10,'LM',opts);sort(diag(10)) > with opts containing disp:0, > > returns different eigenvalues when called repeatedly without changing > anything. > > I'm attaching a.mat that contains the a matrix (and opts, among other > things) > > For now, we are disabling the using of eigs() in Dynare and forcing the > use of eig() > > We are going also to report the problem to Matlab. > > All the best, > > Michel > > (See attached file: a.mat)

1 0

Dynare 4.0.3
by Sébastien Villemot 25 Mar '09

25 Mar '09

Yesterday we released Dynare 4.0.3, which is a bug fixing release. The Windows and Linux packages are available as usual. Michel: it seems that you forgot to upload the source package. Best, Sébastien

1 0

Kirikou
by Sébastien Villemot 25 Mar '09

25 Mar '09

Hi all, Today I migrated our server "kirikou" to the latest Debian stable (codename "Lenny"). Everything went smoothly. The trickiest part was for the wikis: I had to upgrade the data for every wiki from version 1.5 to version 1.7 of Moinmoin. Please let me know if something appears to be broken. Best, Sébastien

1 0

linux mex files for dynare version 3
by Stéphane Adjemian 25 Mar '09

25 Mar '09

Hi Michel, I need to debug something on dynare version 3 (a recent post by glucke on the forum). Unfortunatly the mexglx files are not working on my linux... Where can I find the sources of the mex files? There is nothing on subversion (branches/3.xx). Stéphane.

2 5

2024

2023

2022

2021

2020

2019

2018

2017

2016

2015

2014

2013

2012

2011

2010

2009

2008

Dev March 2009