Commit Graph

453 Commits (e19026d130c384437e9d565505da0b11fb08ca3c)

Author SHA1 Message Date
Houtan Bastani e19026d130 SWZ: remove gsl from build system, to be reintstated if gsl remains a requirement in the ms new code 2011-02-02 18:58:36 +01:00
Ferhat Mihoubi 6f184371ef - Correction of temporary terms when bytecode is used and solve_algo <= 4 2011-02-02 00:31:37 +01:00
Ferhat Mihoubi 66a7884a15 - Corrects the "if (~eval)" displayed with print option in bytecode 2011-01-31 17:44:22 +01:00
Ferhat Mihoubi 3393eebf71 When a model is evaluated with bytecode the residual has to be returned in the equations order (not in the variables order) 2011-01-31 12:30:16 +01:00
Sébastien Villemot 80b5c23346 Build system for MEX files:
- remove the dirty vpath hack
- factorize between the two estimation DLLs
- factorize within SWZ makefile
2011-01-25 17:13:00 -05:00
Sébastien Villemot afeaf1ead0 bump min boost to 1.36 to accord with 72008502c1 2011-01-24 17:52:33 -05:00
Houtan Bastani 0f12abbe9f bump min boost to 1.36 to accord with 72008502c1 2011-01-24 17:35:12 -05:00
Ferhat Mihoubi ce07223628 - The Temporary terms management with bytecode when the model is block decomposed and solve_algo<5 is now compatible with octave
- Gets rid of warning message during the compilation of bytecode
2011-01-14 19:24:18 +01:00
Ferhat Mihoubi 206fff0e71 - Corrects the following problem:
Octave BiCGStab algorithm involves a 0 division in case of a preconditioner equal to the LU decomposition of the A matrix (in a linear system of the form A.x = b).
- The solution:
Checks if the linear system is solved simply using: x_new = x_old + U \ (L \ x_old)
Ticket #11
2011-01-14 10:39:33 +01:00
Sébastien Villemot ffc53ca910 Fix integer types in Kronecker DLLs 2011-01-13 19:34:49 +01:00
Sébastien Villemot e39cb94954 Build system: remove warning about parentheses (we know the priority rules of C++, no need to be reminded) 2011-01-13 19:05:58 +01:00
Ferhat Mihoubi 20f414e0d2 - Extends the error messages in bytecode
- Implements the print option in bytecode
- Minor modifications to manage the global temporary terms
2010-12-31 16:41:50 +01:00
Ferhat Mihoubi 313f64e153 - Adds option 'print' to bytecode
- Manages global temporary terms when the model is evaluated block by block
- Stores the result of a first order derivative only in case of numerical approximation. Do nothing if an external function is called to compute the first order derivatives (it has already been done during the function call)
- Cleans the code
2010-12-31 16:37:34 +01:00
Ferhat Mihoubi 816dbefa28 - Corrections in the debugging conditional preprocessor expressions 2010-12-31 16:29:39 +01:00
Sébastien Villemot 9aa79c6140 Another build system fix 2010-12-28 10:50:51 +01:00
Sébastien Villemot 3680d71fe6 Adapted build system for the previous commit 2010-12-27 10:45:42 +01:00
George Perendia 887209208f Estimation C++ DLL: Adding the new Proposal class with common, adjustable seed for both normal and uniform rng and the associated changes including the removal of now obsolete RandSampler class. 2010-12-24 09:49:59 +00:00
Houtan Bastani a11817cfa4 k-order: added support for m-files, added tests and modified manual 2010-12-22 09:48:05 +01:00
Houtan Bastani d0853e1658 k-order: removed unnecessary check 2010-12-22 09:48:05 +01:00
Ferhat Mihoubi 7a58637932 - Displays an error message and interrupts the execution of bytecode in case of unknown function 2010-12-17 18:54:03 +01:00
Ferhat Mihoubi 129d404eb8 - The binary operator oDerivPower is implemented in bytecode 2010-12-17 18:54:03 +01:00
Houtan Bastani d69d225b82 k-order: removed unused variable from class definition 2010-12-17 12:16:32 +01:00
Sébastien Villemot aa58f7bd37 Estimation DLL, QRDecomposition: removed unused data members 2010-12-17 10:52:00 +01:00
George Perendia 424bd8257b Estimation C++ DLL: Bug fixing in KalmanFilter.cc - re-setting a_init and some minor refactoring and tidying up 2010-12-17 09:02:24 +00:00
Houtan Bastani b63590712c fixed bug in argument number check and removed unused variable 2010-12-14 18:42:29 +01:00
Houtan Bastani e7843e349a Modify configure script to show Kalman Steady State seperately from other mex files 2010-12-14 10:27:39 +01:00
Sébastien Villemot e826fb8239 Build system: kalman_steady_state depends on libslicot and must be built after 2010-12-14 09:57:15 +01:00
Houtan Bastani 41bda1eb07 Mac libslicot compilation optional (because a Fortran compiler is not standard on OS X) 2010-12-13 18:15:11 +01:00
Sébastien Villemot 1329b57ff2 Estimation DLL: fix DecisionRules::compute() (was not reentrant, since GeneralizedSchurDecomposition::compute() modifies D and E) 2010-12-13 17:15:08 +01:00
George Perendia 915bb4c451 Estimation C++ DLL: Bug-fixing Sampler 2010-12-11 19:22:27 +00:00
Michel Juillard d38c4de498 estimation_dll: adding test for logposterior 2010-12-10 21:34:40 +01:00
Ferhat Mihoubi 97664607b2 - external functions are compatible with block decomposition and/or bytecode 2010-12-10 12:06:50 +01:00
Sébastien Villemot befa1b966b Estimation DLL: removed illegal casts from pointer to integer 2010-12-09 10:49:44 +01:00
Sébastien Villemot 965ea7280a Build system: add detection and linking of MatIO library for Estimation DLL under Octave 2010-12-08 10:22:44 +01:00
George Perendia fa3e19fdfb C++ Estimation DLL: Update of logMHMCMCposterior.cc with a draft octave MAT draws file save and adding a test random_walk_metropolis_hastings_core.m: Octave version crashes at start of DLL and Matlab version finishes with low acceptance due to frequent B&K and reports error within debugger too - needs more debugging! 2010-12-07 23:03:11 +00:00
Sébastien Villemot 2d84d8dc49 K-order DLL: code simplification 2010-12-06 17:01:06 +01:00
Sébastien Villemot a9b4331a95 Bugfix at order 3 (hopefully was not affecting final results, thanks to Pablo Winant for spotting this):
- in the 3rd derivatives matrix, among symmetric elements, the first one had
the right values but the following ones were set to zero
- moreover, the k-order DLL was trying to add all the symmetric elements in the
folded tensor, instead of only keeping one value among all the symmetric ones
- hopefully, Ondra's tensor library was (silently) refusing to add symmetric
elements after the first (and right) value had been added
- so the final result was correct
2010-12-06 17:00:46 +01:00
Ferhat Mihoubi bdfa881a96 Reorganization of error handling in bytecode 2010-11-25 16:11:05 +01:00
Stéphane Adjemian (Charybdis) 010c4c0c20 Fixed bug in kalman_steady_state.cc related to long integer and (mkl) lapack calls. Changed building system by adding -fdefault-integer-8 flag to gfortran when slicot is compiled on a 64-bits platform. 2010-11-25 14:35:46 +01:00
Sébastien Villemot 5a9615951c Kronecker DLL: fix minor issue (was creating a lot of warnings under Octave) 2010-11-24 19:12:03 +01:00
Sébastien Villemot 43f69d27e8 Kalman steady state DLL: minor fixes 2010-11-22 12:22:45 +01:00
Sébastien Villemot b8aae7bbf9 Fix underscoring problems with Fortran code 2010-11-22 12:22:22 +01:00
Sébastien Villemot 81342073c0 Build system: added kalman_steady_state DLL sources to dist file 2010-11-22 10:39:19 +01:00
Stéphane Adjemian (Charybdis) dfbbc415cd Merge remote branch 'origin/master' 2010-11-19 14:49:26 +01:00
Stéphane Adjemian (Charybdis) 160922a189 Added mex file for computing the fixed point of the riccati equation in the kalman filter.
The mex works with octave but not with matlab (lapack issue again ;-).
2010-11-19 14:39:54 +01:00
Sébastien Villemot b3ee886391 Build system: added libslicot to dist file 2010-11-19 09:22:00 +01:00
Stéphane Adjemian (Charybdis) d5b032c99f Moved test files for kronecker product mex files from ./mex/sources/kronecker/tests to ./tests/kronecker. 2010-11-18 16:01:21 +01:00
Sébastien Villemot 479b2be602 Build system: added missing PIC flag for Fortran code in Octave DLLs 2010-11-18 11:10:49 +01:00
Stéphane Adjemian (Charybdis) 382ab96cde Added the Slicot library (distributed under the GNU General Public Licence at www.slicot.org). This library provides routines
for the design and analysis of control systems. Adapted the building system to handle fortran 77 source files.
2010-11-18 10:11:27 +01:00
Stéphane Adjemian (Charybdis) 5c709e4733 Changed set_dynare_threads.m. The number of threads in parallelized mex files to be used if dynare is built with the openmp flag
(--with-openmp) is not passed by an environment variable anymore. The function set_dynare_threads changes the default value of the
number of threads (default is 1) in the options_.threads structure. Changed calls to sparse_hessian_times_B_kronecker_C and
A_times_B_kronecker_C dlls accordingly.
2010-11-15 13:06:39 +01:00