Re: [julia-users] poor performance of threads

Thanks a lot for the suggestions. As I have mentioned, it was really a toy problem, but I am not getting a significant speedup on a bigger problem, where threads are nicely separated either and the problem is very CPU bound either. I would be very interested to know about tool that would point o

[julia-users] Re: Sessions in Escher

Hi Leonardo, I'm not an Escher user, but I use SecureSessions.jl when building web apps with plain HttpServer.jl. If this is a little over the top then the source code should be helpful for implementing a simpler scheme. Hope this helps, Jock

[julia-users] Re: malformed expr concerning a module

I think i forgot to add toplevel esc(Expr(:toplevel, ex)) Thanks everyone On Thursday, March 3, 2016 at 6:31:05 PM UTC-8, Julia Tylors wrote: > > Hi everyone; > > I am writing a macro which uses module: > > @f module X > x= 4 > end > > macro f(ex) > isa(ex,Expr) && (ex.head == :module) || erro

[julia-users] Integer enums + some questions

I want to use enums as integers, especially for array indexing. So I started to implement an IntEnum type, so I could do something like the following: julia> @intenum Animal pig dog cat julia> pig pig=1 julia> dog dog=2 julia> A = zeros(length(Animal)) 3-element Array{Float64,1}: 0.0 0.0

[julia-users] malformed expr concerning a module

Hi everyone; I am writing a macro which uses module: @f module X x= 4 end macro f(ex) isa(ex,Expr) && (ex.head == :module) || error("@f : syntax error") module_name = ex.args[2] module_blk = ex.args[3].args import_stmts = quote import CSS.Leech end insert!(module_blk,length(module_blk),impo

Re: [julia-users] ERROR: LoadError: syntax: incomplete: premature end of input

Many thanks for the reply. I re-checked my code. It seems like that I have left out a ")". On Friday, 4 March 2016 02:37:58 UTC+11, Stefan Karpinski wrote: > > It means what entered in new.jl isn't a complete Julia expression. Perhaps > you're missing a closing `end` or something? > > On Thu,

[julia-users] Re: PyPlot ArgumentError: haskey of NULL PyObject error

It sounds like "using PyPlot" probably threw an error that you ignored? To make sure that Python and Matplotlib are properly installed, the easiest thing is to have PyCall install its own Miniconda Python distro: ENV["PYTHON"]="" Pkg.build("PyCall") You only need to do this once, then

[julia-users] Re: Plotting the Rosenbrock Function using PyPlot

This works, for example: x, y = -1.5:0.1:1.5, -1.5:0.1:1.5 z = Float64[rosenbrock([x,y]) for x in x, y in y] using PyPlot pcolor(x,y,z)

[julia-users] Why does using Compose not work behind a firewall with packages in LOAD_PATH

I'm trying to set up julia on lab computers, with the packages shared in a single directory, that is included in the default LOAD_PATH in Julia. For some reason Compose tries to initialize a local METADATA. This fails since git is not allowed. Other packages work fine, including all of Compo

Re: [julia-users] poor performance of threads

Your "one" function is not type-stable (f::Int->f::Float64). From the allocation, two() must also be boxing something, but I'm not sure why. (This might be a julia bug.) In any event, because of all this allocation, it's not telling you anything about what the actual performance will be when it'

[julia-users] subscribe() in Escher

Hi all, I've some doubt about correct positioning of subscribe() in Escher. I've slightly modified 'form' example (attached) to decouple subscribe() and map() call to handle events triggered. Why my example not works as expected, printing all content of form when button is pressed? Many thanks i

[julia-users] Re: subscribe() in Escher

Sorry, I've forgotten the attachment ... Leonardo Il giorno giovedì 3 marzo 2016 23:21:15 UTC+1, Leonardo ha scritto: > > Hi all, > I've some doubt about correct positioning of subscribe() in Escher. > I've slightly modified 'form' example (attached) to decouple subscribe() > and map() call to

[julia-users] Sessions in Escher

Hello, user sessions are supported in Escher? (e.g. with some API to get unique ID) and cookies are supported (or planned)? Many thanks in advance Leonardo

Re: [julia-users] Replacing multiple characters in a string

Here is one solution: str = "ACGT" d = Dict("A"=>"T", "C"=> "G", "G"=>"C", "T"=>"A") y = replace(str, r"[ACGT]{1}", x->d[x]) # "TGCA" reverse() is on the top of that On Thu, Mar 3, 2016 at 3:26 PM Carlos Guzman wrote: > I'm relatively new to programming. Started picking up Julia yesterday, a

Re: [julia-users] Replacing multiple characters in a string

Carlos You are trying to use pattern matching, which is mostly used to find patterns in strings. In your case, it seems the string consists exclusively of these characters, and you want to replace all of them. I would use `map` for this. function dual(c::Char) c=='A' && return 'T' c=='C'

[julia-users] Re: 3D scatter / meshed surface color map

In case you work with an IJulia notebook or Escher, you can use ThreeJS.jl . If you use Escher, you can also create UI's to control the variables and get the plots to update interactively. Please see the examples for more details. There isn't support

[julia-users] ZMQ poll

Hello, if anyone can help me: poll() API exists into ZeroMQ port in Julia? If so, how to call it? Many thanks in advance Leonardo

[julia-users] Re: Plotting the Rosenbrock Function using PyPlot

If you're willing to use Plots, here's the syntax (might be very close to Pyplot, not sure). It uses the pyplot backend using Plots function rosenbrock(x::Vector) return (1.0 - x[1])^2 + 100.0 * (x[2] - x[1]^2)^2 end default(size=(600,600), fc=:heat) x, y = -1.5:0.1:1.5, -1.5:0.1:1.5 z = Surf

[julia-users] Replacing multiple characters in a string

I'm relatively new to programming. Started picking up Julia yesterday, and am having trouble with something I suspect is very simple to do. I've spent some time searching, but haven't been able to get this to work. I am essentially trying to answer this question here: http://rosalind.info/probl

[julia-users] Plotting the Rosenbrock Function using PyPlot

Hi, I have defined Rosenbrock function as *function rosenbrock(x::Vector)* * return (1.0 - x[1])^2 + 100.0 * (x[2] - x[1]^2)^2* *end* *m = rand(51,2)* I want to check the unimodality of the function by plotting(surface as well as contour) over the points given by rows of matrix m. Can anyone h

[julia-users] Re: PyPlot ArgumentError: haskey of NULL PyObject error

this also works on JuliaBox, https://www.juliabox.org/, although with the 0.4 and 0.5 kernels, users will need to add the PyPlot package themselves ... result: In [9]: scatter(x, y, s=area, alpha=0.5) Out[9]: PyObject

Re: [julia-users] poor performance of threads

Tomas In your example, different threads access nearby elements in the array f. They are likely to be in the same cache line, leading to much unnecessary inter-CPU communication. This is also called "false sharing". Leaving a gap between the elements used by each thread might speed things up. You

[julia-users] My module MathToolkit (PSLQ algorithm, recurrence detection, etc.)

Hi, I changed my mind, deleted the previous repository and created a more general package where I want to keep the functions I generally write when using some mathematical software (I wrote the same recvec function in Sympy and Maxima; I also wrote a PSLQ version for Maxima). Right now, I have

[julia-users] poor performance of threads

Hi All, I would like to ask if someone has an experience with Threads as they are implemented at the moment in the master branch. After the successful compilation (put JULIA_THREADS=1 to Make.user) I have played with different levels of granularity, but usually the code was slower or more or les

Re: [julia-users] Using CUDArt on remote machines - 'illegal memory access'

Thanks for the update. On Thursday, March 03, 2016 10:52:26 AM Matthew Pearce wrote: > To get it straight if A is a matrix in main memory, the corresponding GPU > memory object is d_A = CudaArray(A) then: > > A[i, j] = d_A[ j * nrows + i] Yes, that should be right. Please do feel free to edit th

Re: [julia-users] Re: Set Environment Variable to a Number?

That makes sense given what I'm seeing since OMP_NUM_THREADS and other environmental variables work. I guess this isn't the issue then (still don't know why the Phi won't read the environmental variable, instead I have to pass it... but that's a whole separate issue). Thanks On Thursday, March

Re: [julia-users] Re: Set Environment Variable to a Number?

On Thu, Mar 3, 2016 at 1:36 PM, Chris Rackauckas wrote: > Oh, meant the line with the mic environment variable: > > printf("%s\n",getenv("MIC_OMP_NUM_THREADS")); > > It's clear from context though. Sorry for the mistake. Env's are strings and you can't store a number in it. Especially not by le

Re: [julia-users] Using CUDArt on remote machines - 'illegal memory access'

Thanks again. Think the problem may have been with my kernel and getting confused about the row major and column major ordering of the layout of the array. I thought I'd checked it was producing the correct norms yesterday, but I must have changed something... To get it straight if A is a matr

[julia-users] Re: Set Environment Variable to a Number?

Oh, meant the line with the mic environment variable: printf("%s\n",getenv("MIC_OMP_NUM_THREADS")); It's clear from context though. Sorry for the mistake. On Thursday, March 3, 2016 at 10:35:42 AM UTC-8, Chris Rackauckas wrote: > > Hey, > I found a solution for strings here >

[julia-users] Set Environment Variable to a Number?

Hey, I found a solution for strings here , but I really need it to be a number. If I do ENV["MIC_OMP_NUM_THREADS"]=240 Then call a C function, I get printf("%s\n",getenv("OMP_NUM_THREADS")); as 240 whereas printf("%d\n",g

Re: [julia-users] Using CUDArt on remote machines - 'illegal memory access'

Oh, drat, I misread it as if you were trying to create an n-dimensional array of size one in each dimension. I forgot you could move a CPU-array to the host with a construct like that. Sorry about that. More efficient, however, is to use fill or fill! for such things, since there is no memory

[julia-users] Current status of in-place ops

Hi Marcin, Here is an effort under Julia to conclude a final aproach based on actual user code. https://groups.google.com/forum/m/#!topic/julia-users/_LduslpljpE If you like you can contribute too.

[julia-users] Current status of in-place ops

Hello, One of the biggest downsides of Julia is its lack of support for in-place operations. E.g. in numpy I can write import numpy.random as nr A = nr.randn(1, 1) B = nr.randn(3000, 3000) A[0:3000, 0:3000] += B whereas in Julia to achieve the same speed and memory usage I need to writ

Re: [julia-users] Interested in GSoC idea "Random number generation"

Yup, Seems good :)

[julia-users] Re: 3D scatter / meshed surface color map

Ok, thanks. Joaquim quinta-feira, 3 de Março de 2016 às 16:45:20 UTC, Josef Heinen escreveu: > > The build problem is fixed (in the master branch). > peaks() is part of GR - also in the master branch. > > The next GR.jl release will be released this weekend. > > On Thursday, March 3, 2016 at 5:28

[julia-users] Re: 3D scatter / meshed surface color map

The build problem is fixed (in the master branch). peaks() is part of GR - also in the master branch. The next GR.jl release will be released this weekend. On Thursday, March 3, 2016 at 5:28:29 PM UTC+1, J Luis wrote: > > Hmm, is 'peaks' part of GR? Because > > julia> using GR > > julia> z = peak

[julia-users] Re: 3D scatter / meshed surface color map

Hmm, is 'peaks' part of GR? Because julia> using GR julia> z = peaks(); ERROR: UndefVarError: peaks not defined quinta-feira, 3 de Março de 2016 às 16:24:07 UTC, Josef Heinen escreveu: > > I could re-produce the problem on a fresh Windows installation. > Removing the downloaded archive leads t

[julia-users] Re: 3D scatter / meshed surface color map

Oops, don't know what happened to the message asking me what version was it, and my reply to it (both disappeared), but here it goes again. It eas trying to install v0.0.17 julia> Pkg.add("GR") INFO: Cloning cache of GR from git://github.com/jheinen/GR.jl.git INFO: Installing GR v0.9.17 INFO: Bu

[julia-users] Re: 3D scatter / meshed surface color map

I could re-produce the problem on a fresh Windows installation. Removing the downloaded archive leads to a chmod error, although I don't use this function in build.jl. Will have to use another method to delete the file. Nevertheless, the package should be installed On Thursday, March 3, 201

[julia-users] Re: 3D scatter / meshed surface color map

Which version are you using? => *Pkg.installed("GR")* On Thursday, March 3, 2016 at 3:13:33 PM UTC+1, J Luis wrote: > > But unfortunately on Windows > > INFO: Downloading pre-compiled GR 0.17.3 binary > =[ ERROR: GR > ]== > > LoadE

Re: [julia-users] Acessing raw binary data

Could you use BitVectors? They store boolean vectors using one bit per boolean. On Thu, Mar 3, 2016 at 8:29 AM, Christina Lee wrote: > > > I'm interested in looking at how to perform exact diagonalization of > various quantum states in Julia. > > This problem is often memory limited, but the sta

Re: [julia-users] ERROR: LoadError: syntax: incomplete: premature end of input

It means what entered in new.jl isn't a complete Julia expression. Perhaps you're missing a closing `end` or something? On Thu, Mar 3, 2016 at 6:58 AM, Fei Ma wrote: > Hello. Here I have a problem. > > I run a program, but stopped it when it is running, and then the terminal > output this : > >

[julia-users] Re: Creating docstrings using code generation

The `eval` is going to do some really odd things there. If you're on a very recent 0.5 build you should probably see the following output after doing that ``` julia> for name in (:year, :month, :day, :hour, :minute, :second, :millisecond) func = eval(name) @doc """

Re: [julia-users] Using CUDArt on remote machines - 'illegal memory access'

Thanks Tim. For me `elty=Float32' so if I use `CudaArray(elty, ones(10))' or `CudaArray( elty, ones(10)...)' I get a conversion error. [I am running Julia 0.5.0-dev+749] The result of my CudaArray creation above looks like: julia> to_host(CudaArray(map(elty, ones(10' 1x10 Array{Float32,2}:

[julia-users] ERROR: LoadError: syntax: incomplete: premature end of input

Hello. Here I have a problem. I run a program, but stopped it when it is running, and then the terminal output this : ERROR: LoadError: syntax: incomplete: premature end of input in include at ./boot.jl:261 in include_from_node1 at ./loading.jl:304 while loading /home/i/Documents/Julia/new.j

[julia-users] Acessing raw binary data

I'm interested in looking at how to perform exact diagonalization of various quantum states in Julia. This problem is often memory limited, but the states can often be represented by strings of 0 and 1 (binary format). Not only that, but various binary operations, and converting to other bas

Re: [julia-users] julia make error on amazon linux ec2

also gcc version is 4.8.3 20140911 (red hat 4.8.3-9) On 03/03/2016 14:07, Stefan Karpinski wrote: What is the content of /etc/issue – i.e. what is the Linux distro and what is the output of `gcc --version`? On Thu, Mar 3, 2016 at 8:52 AM, Chris Wade mailto:chris.w...@thecitysecret.com>> wrot

Re: [julia-users] julia make error on amazon linux ec2

all it says is amazon linux AMI release 2015.09 kernel \r on an \m On 03/03/2016 14:07, Stefan Karpinski wrote: What is the content of /etc/issue – i.e. what is the Linux distro and what is the output of `gcc --version`? On Thu, Mar 3, 2016 at 8:52 AM, Chris Wade mailto:chris.w...@thecitysecr

[julia-users] Re: Creating docstrings using code generation

Thanks Mike, I also managed to come up with an alternative solution: using Base.Dates for name in (:year, :month, :day, :hour, :minute, :second, :millisecond) func = eval(name) @doc """ $name(dt::TimeType) -> Int64 Return the $name of a `Date` or `DateTime` as an `Int64`.

[julia-users] Re: Creating docstrings using code generation

Just needs extra `$`s in the docstring, one for expression interpolation and one for string interpolation: julia> for name in (:year, :month, :day, :hour, :minute, :second, : millisecond) @eval begin @doc """ $($name)(dt::TimeType) -> Int64

Re: [julia-users] [Plots.jl] Sequential color of consecutive data series

Thank you, this works fine for me. I understand the motivation of generating the palette based on the background color and I don't expect this to be implemented. No worries.

[julia-users] Creating docstrings using code generation

I was hoping to generate several redundant docstrings using code generation. Unfortunately I have run into an issue where $name isn't being replaced in the docstring: julia> using Base.Dates julia> for name in (:year, :month, :day, :hour, :minute, :second, : millisecond) @eval begin

Re: [julia-users] julia make error on amazon linux ec2

On Thu, Mar 3, 2016 at 9:07 AM, Stefan Karpinski wrote: > What is the content of /etc/issue – i.e. what is the Linux distro and what > is the output of `gcc --version`? >From the output,the gcc version seems to be 4.8.3 > > On Thu, Mar 3, 2016 at 8:52 AM, Chris Wade > wrote: >> >> the linux ver

Re: [julia-users] Using CUDArt on remote machines - 'illegal memory access'

I can tell you I've gotten CUDArt-based operations working on remote machines. I don't really know what the issue is, but I did notice a couple of concerns: - what's up with `CudaArray(map(elty, ones(n1)))`? That doesn't look right at all. Don't you mean `CudaArray(elty, ones(n1)...)`? - don't y

[julia-users] Re: 3D scatter / meshed surface color map

But unfortunately on Windows INFO: Downloading pre-compiled GR 0.17.3 binary =[ ERROR: GR ]== LoadError: chmod: no such file or directory (ENOENT) while loading C:\j\.julia\v0.4\GR\deps\build.jl, in expression starting on line 25

Re: [julia-users] julia make error on amazon linux ec2

What is the content of /etc/issue – i.e. what is the Linux distro and what is the output of `gcc --version`? On Thu, Mar 3, 2016 at 8:52 AM, Chris Wade wrote: > the linux version is > 4.1.13-19.31.amzn1.x86_64 GNU/Linux > > > On 03/03/2016 13:42, Stefan Karpinski wrote: > > I was wondering about

Re: [julia-users] julia make error on amazon linux ec2

the linux version is 4.1.13-19.31.amzn1.x86_64 GNU/Linux On 03/03/2016 13:42, Stefan Karpinski wrote: I was wondering about the Linux version, not what AMI. If it's an older version of CentOS, for example, it may have a version of GCC that's too old. On Thu, Mar 3, 2016 at 4:08 AM, Chris

[julia-users] Re: Abstract typed array construction JIT performance

All right that makes sense. On Thursday, March 3, 2016 at 1:50:09 PM UTC, Lutfullah Tomak wrote: > > Hi Colin, > I think if getindex inlined in vect it is not compiled seperately for your > second run. Check code generated with @code_* for youreself.

[julia-users] Abstract typed array construction JIT performance

Hi Colin, I think if getindex inlined in vect it is not compiled seperately for your second run. Check code generated with @code_* for youreself.

Re: [julia-users] julia make error on amazon linux ec2

I was wondering about the Linux version, not what AMI. If it's an older version of CentOS, for example, it may have a version of GCC that's too old. On Thu, Mar 3, 2016 at 4:08 AM, Chris wrote: > Linuz VM is running: Amazon Linux AMI ElasticBeanstalk > 4.1.13-19.31.amzn1.x86_64 GNU/Linux > > On

[julia-users] Re: Using CUDArt on remote machines - 'illegal memory access'

I should add that I don't think the error lies in my `reap' function for remote calls. As I can correctly call the cudakernels.sqrownorms function on the host from the remote: julia> reap(3, :(reap(1, :(sum(cudakernels.sqrownorms(d_M[1] ))[3] 5.149127f6 (The above gets process three to ca

[julia-users] Using CUDArt on remote machines - 'illegal memory access'

Hello I've come across a baffling error. I have a custom CUDA kernel to calculate squared row norms of a matrix. It works fine on the host computer: julia> d_M = residual_shared(Y_init,A_init,S_init,k,sig) CUDArt.CudaArray{Float32,2}(CUDArt.CudaPtr{Float32}(Ptr{Float32} @ 0x000b037a),(40

Re: [julia-users] 3D scatter / meshed surface color map

These are all really cool packages. Don't pick just one! ;) https://github.com/tbreloff/Plots.jl One of these days I'll get GLVisualize working... On Wednesday, March 2, 2016, wrote: > Hello Julia users, > > I would like to know if anyone knows how we can have 3D scatter plots or > meshed surf

Re: [julia-users] Re: JuliaCon 2016: Program Committee Office Hours on the 5th of March at 15:00 UTC

On 3 March 2016 at 11:47, Bart Janssens wrote: > > Should we already fill in the proposal form for JuliaCon before this > session, to improve upon it later, or is it better to just submit after? Either way is fine. I imagine that it is easiest for us on the PC if you have a link to a draft (or d

[julia-users] Re: [ANN] GLVisualize

On Monday, 29 February 2016 16:03:10 UTC+1, Simon Danisch wrote: > > > If not MatplotLib, could this become the Processing (and by extension > OpenFrameworks, LibCinder) of Julia? > > That's definitely more the direction I'd like to take (although with a > very different approach). > I hope that

[julia-users] Re: JuliaCon 2016: Program Committee Office Hours on the 5th of March at 15:00 UTC

Hi, Should we already fill in the proposal form for JuliaCon before this session, to improve upon it later, or is it better to just submit after? Cheers, Bart On Wednesday, March 2, 2016 at 11:55:22 PM UTC+1, Pontus Stenetorp wrote: > > Everyone, > > With the JuliaCon 2016 proposal deadline a

[julia-users] Abstract typed array construction JIT performance

Hi all I have got some questions about array construction JIT speed. abstract Aimmutable B <: A endimmutable C <: A end b = B() c = C() @time getindex(B, b, b) @time Base.vect(b, b); 0.006129 seconds (6.03 k allocations: 276.191 KB) 0.005339 seconds (3.43 k allocations: 175.144 KB) Rest

[julia-users] [Plots.jl] Sequential color of consecutive data series

Hi I was making a plot using Plots.jl and the gadfly backend to draw different series of data on the same plot and I found that if I declare the colour for just a specific data series the default palette advances to the next colour even though the next available colour in the palette has not be

Re: [julia-users] julia make error on amazon linux ec2

Linuz VM is running: Amazon Linux AMI ElasticBeanstalk 4.1.13-19.31.amzn1.x86_64 GNU/Linux On Wednesday, March 2, 2016 at 6:17:25 PM UTC, Stefan Karpinski wrote: > > What version of Linux is the VM running? > > On Wed, Mar 2, 2016 at 12:06 PM, Chris > wrote: > >> Im trying to install julia on am

[julia-users] Re: julia make error on amazon linux ec2

Linuz VM is running: Amazon Linux AMI ElasticBeanstalk 4.1.13-19.31.amzn1.x86_64 GNU/Linux On Wednesday, March 2, 2016 at 5:06:07 PM UTC, Chris wrote: > > Im trying to install julia on amazon linux instance. > > i ran: > > sudo yum -y install clang m4 patch ncurses-devel python-devel > git clone

[julia-users] Re: PyPlot ArgumentError: haskey of NULL PyObject error

Works find on my machines. Does you example work in the Python environment you are using within Julia? On Thursday, March 3, 2016 at 8:18:52 AM UTC+1, Kai Xin Thia wrote: > > Using the sample code below, Version 0.4.3 on Mac, I keep getting the > error. Is there a way to fix it?: > > > using PyP

[julia-users] Re: 3D scatter / meshed surface color map

*Pkg.add("GR")* *using GR* *z = peaks()* *surface(z)* On Thursday, March 3, 2016 at 5:10:52 AM UTC+1, mauri...@gmail.com wrote: > > Hello Julia users, > > I would like to know if anyone know