[julia-users] CPU_CORES is deprecated

2016-10-28 Thread Nils Gudat
As the title says - getting a warning but the constant still exist. There's none of the ususal "use ... instead" in the warning, what's the replacement (if any) for this in 0.5?

[julia-users] Re: Segfaults in SharedArrays

2016-06-26 Thread Nils Gudat
Oh well that would explain why it had no effect then - is there some version that contains this fix then? Would building the current 0.4.7+ master from source do? (Sorry, I've stuck strictly to stable versions so far so am not well versed in how the master branches work)

[julia-users] Re: Segfaults in SharedArrays

2016-06-25 Thread Nils Gudat
No luck unfortunately on 0.4.6 either, so it seems the SharedArray fix (assuming it made it into 0.4.6) didn't help. Any other ideas? How does one figure out the source of a segfault?

[julia-users] Re: How to install 0.4.5 on Ubuntu?

2016-06-16 Thread Nils Gudat
Apologies for decreasing quality of questions, but how do I actually start julia after building from source on Linux? After having successfully (I believe) built the latest 0.4.6-pre master, I can't find an exectuable to start julia. I tried grepping through the folder (find | grep '\.exe$'),

Re: [julia-users] Segfaults in SharedArrays

2016-06-15 Thread Nils Gudat
That looks familiar, yes. Would this be available in 0.4.6 then? And if so, when is 0.4.6 expected to be released?

[julia-users] Segfaults in SharedArrays

2016-06-15 Thread Nils Gudat
I've had a longstanding problem with random segfaults in a program I've written which recently re-appeared (see original discussion here: https://groups.google.com/forum/#!topic/julia-users/FV9AmKsZNuk) It appears that currently these segfaults are caused by the usage of SharedArrays (or at

[julia-users] Re: How to install 0.4.5 on Ubuntu?

2016-06-14 Thread Nils Gudat
What (might have) helped in the end was to remove the nightlies PPA - I had added both PPAs for release version and nightly, so maybe that was my mistake.

[julia-users] (How?) Has parallel execution changed in 0.5?

2016-06-13 Thread Nils Gudat
I feel like I've posted on here with the same title when we went from 0.3 to 0.4, but I'm finding myself in a similar situation. When running my code on 0.5, it simply stops at a function definition - there's no error, the program just doesn't move on. The function definition is this:

[julia-users] Re: How to install 0.4.5 on Ubuntu?

2016-06-13 Thread Nils Gudat
Sorry to come back to this, but I've now tried the instructions given above (which are also on the julialang website as I've now seen), but I keep getting v0.5.0-2975-ubuntu14.04.1. When I do sudo apt-get install julia=0.4.5 as before, I'm being told that version 0.4.5 is not found. What could

Re: [julia-users] Re: Double free or corruption (out)

2016-06-12 Thread Nils Gudat
So it looks like I'm having the same issue - have been running the code without parallelization (defining my SharedArrays as regular ones), and it has now been going for about 3 days without any segfaults. Is this a known issue? If so, do we know whether there's a Julia version one can revert

[julia-users] Re: ANN: DC.jl - automagical linked plots in your IJulia Notebook

2016-06-08 Thread Nils Gudat
Looks great, and after some fights with different csv readers in Julia (is there still none that can compete with pandas/fread?) I got it to work. Two questions after a couple of minutes of trying it out - note that I have zero experience with DC.js or javascript in general! - The package

Re: [julia-users] Re: Can I somehow get Julia without standard library?

2016-06-07 Thread Nils Gudat
Don't know whether this is something you'd be interested in, but Scott Jones has worked on building a Julia "lite" version with fewer components: https://groups.google.com/forum/?fromgroups#!searchin/julia-users/raspberry$20pi/julia-users/WStpLtrKiFA/JhiAbc-vAwAJ

Re: [julia-users] Re: Double free or corruption (out)

2016-06-02 Thread Nils Gudat
Hm, interesting observation... I suppose the issue in my case is that the code as it is takes about 3-4 days to complete, so running it on 1 instead of 15 cores means I'm unlikely to ever get my PhD! I will at least try to run a shorter version that might be solvable in a day or two without

Re: [julia-users] Re: Double free or corruption (out)

2016-06-02 Thread Nils Gudat
Fair enough. Does anyone have any clues as to how I would go about investigating this? As has been said before, the stacktraces aren't very helpful for segfaults, so how do I figure out what's going wrong here?

[julia-users] Re: julia equivalent of python [] (Part II)

2016-06-01 Thread Nils Gudat
A roundabout way: lol = ["a",1,"b",(2,3)] lol[4] = collect(lol[4]) Output: 4-element Array(Any,1): "a" 1 "b" [2,3]

Re: [julia-users] Re: Double free or corruption (out)

2016-05-31 Thread Nils Gudat
machine with Julia 0.4.5 installed an all packages on the latest tagged versions. On Tuesday, September 29, 2015 at 1:43:31 PM UTC+1, Nils Gudat wrote: > > The code usually segfaults after 2-5 hours, and is available at > http://github.com/nilshg/LearningModels, however I haven't writ

[julia-users] Re: How to install 0.4.5 on Ubuntu?

2016-05-22 Thread Nils Gudat
Found it, sudo-apt get install julia=0.4.5-trusty did the trick!

[julia-users] How to install 0.4.5 on Ubuntu?

2016-05-22 Thread Nils Gudat
Slightly stupid question from someone who is forced to use Linux on a server: how can I get the latest stable version as opposed ot the dev release? I followed the instructions from here (http://julialang.org/downloads/platform.html), but always end up with a 44-day old 0.5 master version

[julia-users] Re: Parallel computing

2016-05-22 Thread Nils Gudat
This is a rather broad question, and I would suggest that you have (another?) close look at the parallel documentation. The errors about functions not defined on worker processes e.g. suggest that you aren't using the

[julia-users] Re: Help! Problem with repeated calling of a function that creates SharedArray.

2016-05-13 Thread Nils Gudat
This code runs for me without problems (Julia 0.4.5 64bit on an 8GB Win8.1 laptop)

[julia-users] Process medium-sized files in JuliaBox (or alternative)

2016-05-13 Thread Nils Gudat
Slightly weird scenario: Have to do a small bit of analysis on a heavily locked down machine, meaning the only way for me to use any type of decent software is to use it online. First port of call was obviously JuliaBox, but unfortunately I need to process a couple of csv files on the order of

[julia-users] Re: Executing a Julia script from python

2016-02-04 Thread Nils Gudat
If you have a few minutes watch this video, in which Stefan explains how to use Julia from within Python: https://www.youtube.com/watch?v=PsjANO10KgM (you'll be especially interested in the bit starting around minute 25)

[julia-users] Re: how to get back the objective value, the number of explored nodes, number of simplex iterations from gurobi solver

2016-02-02 Thread Nils Gudat
Are you using Gurobi.jl?

[julia-users] Re: What happens when I do x[:] when x::SharedArray{Float64,1}?

2016-01-11 Thread Nils Gudat
I actually remember this changing around v0.3.8(?), as I had a few similar(SharedArray(..)) calls in a buried in a larger model and spend some time trying to figure out why the results suddenly made no sense at all... But I think my original question still stands: what does the [:] call

[julia-users] Re: What happens when I do x[:] when x::SharedArray{Float64,1}?

2016-01-11 Thread Nils Gudat
Oh sorry, I misinterpreted you saying [:] was "equivalent" to the tmp/similar construction. Thanks!

[julia-users] What happens when I do x[:] when x::SharedArray{Float64,1}?

2016-01-09 Thread Nils Gudat
Usually when doing calculations on SharedArrays, I would extract their contents using sdata() to be able to pass them to functions only defined for regular arrays, and then write the return value back into the SharedArray. However, I just realized that instead of doing say: shared[:] =

[julia-users] Re: Pkg.add() though a corporate network proxy won't work, how to fix? - Windows Environment

2015-12-29 Thread Nils Gudat
This thread got me halfway to configuring Julia to update packages through my work network. The missing piece was adding my workplace's security certificate to the list of certificates used by my git installation, which is very well explained in this post

[julia-users] Re: Nettle has build errors - no provider can provide libnettle on Win 10

2015-12-20 Thread Nils Gudat
>>> >>> Please post the entire output of what happens when you try >>> Pkg.build("Nettle"). >>> >>> >>> On Saturday, December 19, 2015 at 2:17:35 PM UTC-8, Nils Gudat wrote: >>>> >>>> Getting this error when try

[julia-users] Nettle has build errors - no provider can provide libnettle on Win 10

2015-12-19 Thread Nils Gudat
Getting this error when trying to add IJulia to my 0.4.2 installation. I found this old IJulia issue which in turn references this BinDeps issue as a possible solution, but both of these

[julia-users] Re: Nicer syntax collect(linspace(0,1,n))?

2015-10-01 Thread Nils Gudat
Given that I had been using convert(Array{Float64,1}, linspace(0,1,n)) so far, I think the collect() syntax is great :)

[julia-users] Re: How to import all symbols from a module with @everywhere?

2015-09-26 Thread Nils Gudat
The way to do this is addprocs(CPU_CORES-1) import Images, Colors @everywhere using Images, Colors I have to admit that I still don't understand the difference between import and using, but I think import doesn't make the exported functions available in the namespace of each worker, for that

Re: [julia-users] Re: Double free or corruption (out)

2015-09-26 Thread Nils Gudat
The minimization itself is NLopt, the problem is to solve an economic model (which takes around 2 minutes to solve on 16 cores) and compare its output (a 100x4 Float64 Array) to some data moments. The model results depend on two parameters. The model itself is mostly minimization (via Optim)

[julia-users] Re: Double free or corruption (out)

2015-09-26 Thread Nils Gudat
That's the problem I alluded to in my question: This happened in the middle of a very lengthy minimization problem, which had been running for a couple of hours. On a previous run, a very similar version of the code finished successfully after about 10 hours. I was hoping that someone could at

[julia-users] Blosc has build errors

2015-09-23 Thread Nils Gudat
When trying to precompile HDF5, I'm getting: This is 0.4.0-pre+7340 on Ubuntu 14.04.3 Any ideas on how to solve this?

[julia-users] Re: Blosc has build errors

2015-09-23 Thread Nils Gudat
Yep, that did the trick! That's what happens when Windows users have to use Linux on servers... Thanks!

[julia-users] (Why?) Has the behaviour of similar() changed in 0.4.0-rc1?

2015-09-18 Thread Nils Gudat
When working with multiple SharedArrays, I used to be able to do: julia> x = SharedArray(Float64,1) 1-element SharedArray{Float64,1}: 0.0 julia> similar(x) 1-element SharedArray{Float64,1}: 0.0 When I got to the office this morning, I was close to a heart attack, as all my model outputs that

[julia-users] Re: Question about parallel computing

2015-09-18 Thread Nils Gudat
Your question seems to be a bit generic - there's lot of different ways in Julia to do parallelized computing, depending on your problem and the hardware setup you're trying to use. If you just search this forum for "parallel", you'll find lots of threads, many of them containing code snippets

[julia-users] Re: (Why?) Has the behaviour of similar() changed in 0.4.0-rc1?

2015-09-18 Thread Nils Gudat
Ah, sorry, indeed! As the issue mentions, I wanted to highlight the fact that similar now returns an Array, not a SharedArray. I guess this makes sense, I had just never used SharedArrays for anything other than results collection (i.e. I do all other calculations on sdata(SharedArray)), so I

[julia-users] Re: Juno stopped working - error message

2015-09-18 Thread Nils Gudat
I'd say just switch to JunoAtom, it's much nicer :) But in all seriousness, these "utf8sizeof" related bugs can often be solved by doing Pkg.checkout("JuliaParser")

Re: [julia-users] Re: IDE for Julia

2015-09-17 Thread Nils Gudat
I've been using the Atom client for a couple of weeks now an I think it's brilliant - there's a lot of active development going on currently with many quick bugfixes almost daily and lots of great mid- to long-term feature additions planned. It'd be great if the community could focus their

Re: [julia-users] Debugging help with NLopt

2015-09-17 Thread Nils Gudat
Argh, I can't believe I actually made this mistake - this is the umptiest minimum distance estimator I'm writing using NLopt, and I always ensure that a scalar instead of a vector is returned by doing exactly what you're suggesting, but for whatever reason I simply missed it this time around.

[julia-users] Debugging help with NLopt

2015-09-17 Thread Nils Gudat
Just putting this out here in the hope that someone has an idea: I'm trying to fit a model to some targets using a minimum distance estimator with NLopt. The objective function basically computes the distance between model output given two parameters and the targets I set. The structure is:

[julia-users] Re: Juno + Julia 0.4

2015-09-15 Thread Nils Gudat
Yes, it does (that's the point of having the bundle(!) available on the website). To connect to a different instance of Julia, follow the steps here: http://junolab.org/docs/install-manual.html However, I don't think JunoLT is working with 0.4 (at least it wasn't when I last checked two weeks

[julia-users] Re: How to append data to a data frame?

2015-09-14 Thread Nils Gudat
Works the same way it would in pandas: df = DataFrame(col1 = rand(10), col2 = rand(10)) add a column by doing: df[:col3] = rand(10)

[julia-users] Re: How to append data to a data frame?

2015-09-14 Thread Nils Gudat
Sorry for not actually paying attention to your question, you want to add a row! If you can get the data into another dataframe with the same column names (e.g. if you have another CSV with the same header), you can do: append!(df, newdata)

[julia-users] Re: META: What are the chances of moving this forum?

2015-09-12 Thread Nils Gudat
I'll have two weeks off at the beginning of October, during which I'd be happy to read up on how to get a Discourse forum up and running.

[julia-users] Re: How to evaluate variable inside Formula type? (Linear Mixed Models)

2015-09-11 Thread Nils Gudat
Do you mean something like this? f(x) = x^2 y = 5 f(eval(parse("y"))) # = 25 this converts the string to a symbol, which is then evaluated (giving the 5 previously assigned to it). Of course the standard disclaimer applies; you should probably put all those "X1045" variables in an Array and

[julia-users] Re: product of vectors and matrix elements

2015-09-10 Thread Nils Gudat
You've (almost) already written out a list comprehension syntax to do this. Generally I think you wouldn't want to have a couple of vectors and a couple of matrices, but the vectors stored as (e.g.) rows of a matrix, and the matrices stored as slices of a three-dimensional array. E.g.: vectors

Re: [julia-users] Re: META: What are the chances of moving this forum?

2015-09-10 Thread Nils Gudat
The Rust example probably explains this: https://youtu.be/pYYmnBma9qk?list=PLP8iPy9hna6Sdx4soiGrSefrmOPdUWixM=337 :)

Re: [julia-users] Re: META: What are the chances of moving this forum?

2015-09-10 Thread Nils Gudat
Juno is already on (I think) discourse: http://discuss.junolab.org/about That's what already brought me to post this, as the forum over there is infinitely better than Google Groups. Maybe Mike has some comments about time & monetary effort involved in setting this up? On Thursday, September

[julia-users] Re: cleaning up objects in parallel processes?

2015-09-09 Thread Nils Gudat
I think whos() should help you, it'll give a list of all objects defined, as well as their size in memory. If you run it as remotecall_fetch(2, whos, ), where 2 is the number of the worker process (of course you could pick any number returned by procs()), you should be able to figure out what's

[julia-users] Re: a floating point test for other's environments

2015-09-09 Thread Nils Gudat
Not sure that this is what you meant by "(shh, it's new hardware)", but on a Surface Pro 3 running Windows 10 the test passes.

[julia-users] META: What are the chances of moving this forum?

2015-09-09 Thread Nils Gudat
Was just thinking about this as first I had to try opening a thread here a couple of times before any posts were actually displayed, and then read through this thread on Mike's Juno Discuss forum,

[julia-users] Re: META: What are the chances of moving this forum?

2015-09-09 Thread Nils Gudat
Not sure about this (maybe Mike can chime in), but I believe all the categorization is from the users themselves, as the forum software asks you to put your post in a category when you ask it. Also, I wasn't implying that the auto tips alone would be a reason to move, the slowness of Google

Re: [julia-users] Re: META: What are the chances of moving this forum?

2015-09-09 Thread Nils Gudat
That's a good point, I hadn't thought of that. Discourse (which I think is what's behind the Juno forum, and the Atom discussion group) seems to do: Discourse Mailing list mode Discourse reply via email

[julia-users] Re: cleaning up objects in parallel processes?

2015-09-07 Thread Nils Gudat
Aren't you locally creating state on each of the worker processes?

[julia-users] Re: cleaning up objects in parallel processes?

2015-09-06 Thread Nils Gudat
Not entirely sure about this, but wouldn't you have to first re-allocate those large arrays before gc() can free up the memory? That's how I tend to do it, based on what it says in the manual here : "For

[julia-users] Why does 0.4 error silently in parallel?

2015-09-04 Thread Nils Gudat
In 0.3.11, when I would (as usual) run some idiotic code with lots of typos in it, I'd be flooded with error messages from all workers. Now, on 0.4, I am often puzzled by weird results, which I can trace down after a while to the fact that some arrays that should contain results are actually

[julia-users] Re: @parallel and HDF5 is posible to work together?

2015-09-01 Thread Nils Gudat
As the warning says, the module is not defined on the workers. Do addprocs() import HDF5 @everywhere using HDF5

Re: [julia-users] IDE for Julia

2015-08-31 Thread Nils Gudat
Just check the contributor list here: https://github.com/JunoLab/atom-julia-client/graphs/contributors Let's say it's fairly strongly tilted towards Mike at the moment (who also originally built most of the Juno LightTable IDE) :)

Re: [julia-users] IDE for Julia

2015-08-31 Thread Nils Gudat
I've been using JunoLT for about a year and switched to the new Atom client about 2 weeks ago - I think it's already really good, despite a long to-do list and some rough edges that need ironing out.

[julia-users] Multidimensional linear extrapolation on irregular grid

2015-08-31 Thread Nils Gudat
Is there a package that can do a multidimensional (3D and higher) linear interpolation on an irregular grid, and extrapolate values beyond gridpoints? >From what I can see: - Grid.jl can do linear extrapolation, but only on regular grids - The same is true for Interpolations.jl - ApproXD

[julia-users] Re: Multidimensional linear extrapolation on irregular grid

2015-08-31 Thread Nils Gudat
Sounds interesting, I'd definitely give it a go! I'd be doing 3 to 6 dimensional interpolation, but I'm sure a 3D example would scale up easily!?

[julia-users] Re: Multidimensional linear extrapolation on irregular grid

2015-08-31 Thread Nils Gudat
I thought that Interpolations.jl could do this (it's ticked off the feature wishlist here: https://github.com/tlycken/Interpolations.jl/issues/5), but from what I understand this only works for regularly spaced data?

[julia-users] Re: Multidimensional linear extrapolation on irregular grid

2015-08-31 Thread Nils Gudat
I wish I wouldn't have to submit my PhD thesis in 5 weeks, between this and Mike's new work on the julia atom client there's never been so much stuff that I would happily try to contribute too... I'll see how the QuantEcon bit works, but otherwise I might just have a go at the gradient

[julia-users] Re: help with parallelization

2015-08-29 Thread Nils Gudat
There could be different problems with your code: 1. dijkstra_shortest_paths might not be defined on the remote workers (your snippet doesn't show where you define this); you can check this by doing remotecall_fetch(3,isdefined,:dijkstra_shortest_paths) 2. d has to be a SharedArray if you want

Re: [julia-users] Extra '\n' char while creating a string spaning over two lines

2015-08-28 Thread Nils Gudat
I must admit that I don't entirely understand your original example, so I'm not sure this will help, but given that - as you pointed out - you can add linebreaks into code, why don't you create your string as this is a* multiline string

[julia-users] Re: (How?) Has parallel computing changed in 0.4?

2015-08-28 Thread Nils Gudat
Thanks! So this means the issue is fixed already? I just installed the latest 0.4 nightly (Windows), but performance hasn't changed.

[julia-users] Re: (How?) Has parallel computing changed in 0.4?

2015-08-28 Thread Nils Gudat
0.4.0-dev+7053, commit ff77f73 (2015-08-28 04:25 UTC) Windows (x86_64-w64-mingw32)

[julia-users] Re: (How?) Has parallel computing changed in 0.4?

2015-08-27 Thread Nils Gudat
Thanks for the clarification :) Any pointers to where I could read up on the reasons for this, and (maybe?) ways around it?

[julia-users] Re: (How?) Has parallel computing changed in 0.4?

2015-08-27 Thread Nils Gudat
Great, this does the trick! So as a general rule, for parallel processing the way to go would be 1. addprocs() 2. import all modules used 3. @everywhere include auxiliary files such as function definitions, which should include a using statement for all modules so that they are in the

[julia-users] Re: (How?) Has parallel computing changed in 0.4?

2015-08-27 Thread Nils Gudat
Now that my code is running, in case anyone is still following this thread: is it possible that there has been a regression in parallel speed in 0.4? I recall that running my code in 0.3.11 resulted in a more or less linear speedup on my local machine (i.e. adding seven processors on my machine

[julia-users] (How?) Has parallel computing changed in 0.4?

2015-08-27 Thread Nils Gudat
Might be in relation to this thread https://groups.google.com/forum/#!msg/julia-users/FjGXSTzvfmc/j0ZDG629IwAJ - I've ported some code from 0.3.11 to 0.4, but when trying to run it in parallel I encounter all kinds of issues. My main file starts with nprocs()==CPU_CORES ||

[julia-users] Multiple dispatch confusion

2015-08-26 Thread Nils Gudat
I'm defining a function with two methods as follows: a = 2. b = 0.8 function u(c::Float64, h::Float64, b=b, a=a) ((c/(h^b))^(1-a))/(1-a) end function u(c::Float64, a=a) c^(1-a)/(1-a) end However, when calling u(2.), the result is -0.87..., which should be the result of the function

Re: [julia-users] Re: PyPlot: How to change the aspect ratio of a plot?

2015-08-26 Thread Nils Gudat
https://github.com/stevengj/PyPlot.jl#exported-functions

Re: [julia-users] Re: PyPlot: How to change the aspect ratio of a plot?

2015-08-26 Thread Nils Gudat
FWIW, I always think the fig/ax syntax aligns best with Python, e.g. fig, ax = PyPlot.subplots(2,2, ...) ax[1,1][:plot](..., linestyle = :, label = ...) ax[1,1][:legend](loc=best) ax[1,2][:hist](..., alpha = 0.5, label = ...) fig[:title](Some title) with this syntax, you can more or less convert

[julia-users] Re: parallel threads broken, replacing module

2015-08-25 Thread Nils Gudat
Can I just briefly revive this thread, as I don't want to spam the issue on GitHub? I just switched over from 0.3.11, and I'm having some issues with the parallel parts of my code. I have a couple of functions written in .jl files which I include at the beginning of my code, with each of these

[julia-users] Weird error in @everywhere begin...end block

2015-08-21 Thread Nils Gudat
I'm reading in an h5 file using the HDF5 package and need to invert it over the first two dimensions after reading it in, and I'm doing so in a @everywhere begin...end block, as the data has to be available on all workers. For some reason, doing both steps within one block fails. This works as

[julia-users] Re: Weird error in @everywhere begin...end block

2015-08-21 Thread Nils Gudat
NB: I just saw a post on the Juno forum http://discuss.junolab.org/t/compile-error-when-using-statement-and-for-loop-are-inside-begin-block/295 about this, so it might just be a Juno issue.

Re: [julia-users] Embarrassingly parallel workload

2015-08-20 Thread Nils Gudat
Sebastian, I'm not sure I understand you correctly, but point (1) in your list can usually be taken care of by wrapping all the necessary usings/requires/includes and definitions in a @everywhere begin ... end block. Julio, as for your original problem, I think Tim's advice about SharedArrays

[julia-users] LightXML has build errors

2015-08-20 Thread Nils Gudat
Has anyone seen this? - INFO: Building LightXML - INFO: Packages to update: zlib1 INFO: Packages to install: libxml2-2 - INFO: Downloading: zlib1 - INFO: Extracting: zlib1 - 7-Zip [64] 9.20 Copyright (c) 1999-2010 Igor

[julia-users] Why does enumerate fail in parallel?

2015-08-19 Thread Nils Gudat
I just rewrote one of my parallel loops and was surprised by this: list = [a, b, c] for i in enumerate(list) println(i) end (1,a) (2,b) (3,c) addprocs(2) @sync @parallel for i in enumerate(list) println(i) end ERROR: `getindex` has no method matching

[julia-users] What am I doing wrong when using QuantEcon's ARMA type?

2015-08-17 Thread Nils Gudat
I'm cleaning up a bit of code I wrote and wanted to replace a naive loop that simulates an AR(1) process with QuantEcon's built in ARMA type. However, I'm not getting the results I expected: ρ = 0.8 0.8 var_ɛ = 0.03 0.03 theoretical_variance = var_ɛ/(1-ρ^2) 0.08336 ar1 =

[julia-users] Re: What am I doing wrong when using QuantEcon's ARMA type?

2015-08-17 Thread Nils Gudat
Ah, I should have checked the docs more closely. I was confused because I went by the fields listed in the type defintion of ARMA (e.g. here in the repo) https://github.com/QuantEcon/QuantEcon.jl/blob/master/src/arma.jl#L48, and there it says variance of the white noise process (I think lines

[julia-users] Re: How has the behavior of include/require etc. changed in 0.4?

2015-08-15 Thread Nils Gudat
It seems to solve the issue regarding functions from base (upon checking, mean now has 70 methods after using QuantEcon), but the LoadError that's displayed in my original post persists. I guess that was expected though, given that it's unrelated to the problem of extending functions, and I

Re: [julia-users] Function for checking monotonicity of higher-dimensional arrays

2015-08-12 Thread Nils Gudat
Thanks for all the suggestions! Since I'm still on 0.3.11, and only have to check a couple of fairly small arrays, I simply went with function checkmonotonicity(v::Array{Float64}) counter = Array(Int64, ndims(v)) for i = 1:ndims(v) counter[i] = sum(mapslices(issorted, v, i)) end

[julia-users] Function for checking monotonicity of higher-dimensional arrays

2015-08-11 Thread Nils Gudat
This is a re-post of sorts of a question I already asked a year ago in the user group. At the time, Stefan Karpinski proposed that base.issorted might be expanded to accomodate checking in multiple dimensions, but looking at the docs it seems this hasn't happened. This post now has two

[julia-users] Re: How can I change R code to Julia code. It is always error. Thank you!

2015-08-03 Thread Nils Gudat
I think you really might want to take some time to go through the documentation and learn some of the basics of Julia, a good starting point might be http://julia.readthedocs.org/en/latest/manual/noteworthy-differences/ The first error is simply telling you that the result of

[julia-users] Re: How can I change the R code to Julia code. Thank you.

2015-07-31 Thread Nils Gudat
As your code is not a reproducible example, there's no way to say what's going wrong exactly. The error message is telling you that you are trying to access an array out of bounds, so it would be necessary to know the sizes of the arrays in your code to find the error. Two other things: -

[julia-users] Re: Irregular Interpolation

2015-07-29 Thread Nils Gudat
Besided GridInterpolations, you might want to look at ApproXD.jl, which works up to four dimensions. I once wrote a little script that compares different one- and two-dimensional interpolation schemes in Julia, it can be found here

[julia-users] Re: How Do You Filter Rows in a DataFrame Based on a Set of Interest or Specific Pattern?

2015-07-23 Thread Nils Gudat
I'm not a pro when it comes to DataFrames, but for your first application you might want to look at findin? I believe you could achieve your goal by doing: data_frame_value_in_set = data_frame[findin(data_frame[:quality], [7,8]), :] This should return all rows where quality is either 7 or 8,

Re: [julia-users] Extracting values from Array{ANY}

2015-07-19 Thread Nils Gudat
I think this is cross-posted from here: http://stackoverflow.com/questions/31498869/extract-values-from-arrayany-in-julia The text file is not in a standard format, so readdlm inserts a bunch of #undef elements into the resulting matrix. It's not really a problem with findfirst. On Sunday,

[julia-users] Re: Parallel execution, again

2015-07-08 Thread Nils Gudat
Please disregard the first part of my post, the @sync @parallel approach does work (after fixing a simple typo that actually had nothing to do with Julia's parallelism design). I'd still be interested in the second question re: format of functions that are passed to pmap.

[julia-users] Parallel execution, again

2015-07-08 Thread Nils Gudat
After succesfully parallelizing a larger number of small tasks using the @sync @parallel approach and SharedArrays, I'm currently stumped trying to parallelize a small number of large tasks. I'm basically trying to compare the results of minimization procedures using different NLopt

[julia-users] Problems with NLopt

2015-06-23 Thread Nils Gudat
After using the Optim package for quite a while, I'm now dabbling with NLopt.jl for the first time to solve a bounded multivariate optimization problem. For some reason, I can't seem to get NLopt to do anything though, and while the problem is hard to recreate as a small working example, I

[julia-users] Multivariate optimization with bounds

2015-06-18 Thread Nils Gudat
I'm trying to minimize a function of multiple variables using the Optim package. In my original Matlab code, I'm supplying two arrays to fmincon to set upper and lower bounds on each of the variables, but it seems the optimize function in Optim only allows for bounds in the univariate case.

[julia-users] Re: Need help writing parallel code with @sync and @everywhere

2015-06-17 Thread Nils Gudat
I haven't used @everywhere in combination with begin..end blocks, I usually pair @sync with @parallel - see an example here https://github.com/nilshg/LearningModels/blob/master/NHL/NHL_6_Bellman.jl, where I've parallelized the entire nested loop ranging from lines 25 to 47.

[julia-users] Re: Please sdvise how to plot the matrix (10x1000) with different color

2015-06-17 Thread Nils Gudat
It's not clear to me from your post what you are trying to do - do you want to plot each of the 10 rows of vec as a line with a different color? You might want to read this discussion https://github.com/dcjones/Gadfly.jl/issues/526 suggesting that it is much easier to to the kind of plot you

Re: [julia-users] Please sdvise how to plot the matrix (10x1000) with different color

2015-06-17 Thread Nils Gudat
^ True, I assumed in my post above that the rows where the object of interest here, just because it seems more natural to plot 10 1000 points data series than 1000 data series with 10 observations each... On Wednesday, June 17, 2015 at 3:35:18 PM UTC+1, Miguel Bazdresch wrote: Note that in

[julia-users] Re: Some DataFrames questions

2015-05-20 Thread Nils Gudat
I think I have to give up and grudgingly revert to pandas/R - I just tried to do this within a loop, dropping observations based on comparisons of a number of columns numbered by years with some transformations of other columns in the corresponding year. This is my (failed) attempt: for i =

  1   2   >