[julia-users] Re: brfft or irfft help

That does it! I never would have guessed that dims would refer to the side length of the matrix. Thanks a ton. On Wednesday, June 18, 2014 11:47:47 PM UTC-7, Alex wrote: > > Hi Ethan, > > irfft and brfft take an additional argument d, which roughly speaking > denotes the size of the transformed

[julia-users] Re: brfft or irfft help

Hi Ethan, irfft and brfft take an additional argument d, which roughly speaking denotes the size of the transformed matrix. See also the docs on brfft and irfft The following should work for you julia> a = rand(64

Re: [julia-users] Using Image.jl for false-color maps

How do you extract the colors from the colormap? So if I have some colormap: a = ColorMap([RGB(1,0,0),RGB(0,0,1)],100,1.) How do I get the vector of colors (should contain 100 colors from red to blue)? On Sunday, September 22, 2013 2:46:48 AM UTC+10, Steven G. Johnson wrote: > > See also https:/

[julia-users] brfft or irfft help

Does anyone have experience with brfft or irfft? I’m trying to optimize some code and I noticed a huge performance gain if I use rfft over fft for 2-d real matrices. However I need to filter in the Fourier domain, then use irfft (or brfft) to return to spatial corridinates but I can’t seem to

Re: [julia-users] Another Julia/Matlab timing comparison

Yeah, I could try a PR once I get more familiar with the behavior and get acquainted with GitHub. Cheers!

Re: [julia-users] Another Julia/Matlab timing comparison

On Wednesday, June 18, 2014 06:42:41 PM Ethan Anderes wrote: > Ah, ok. That's what I needed. Thanks a ton Tim! If you wouldn't mind distilling what you found missing and submitting a pull request for the docs, that would be awesome. --Tim

Re: [julia-users] Another Julia/Matlab timing comparison

Ah, ok. That's what I needed. Thanks a ton Tim!

[julia-users] Re: uncurried functions are ok... if you know what you are doing

See here (https://github.com/JuliaLang/julia/issues/554) for discussion on currying in Julia. On Wednesday, June 18, 2014 6:01:38 PM UTC-5, Peter Simon wrote: > > How about > > pairs = [e for e in enumerate(sqr)] > filter(p -> +(p...)%4 != 0, pairs) > > ? > > --Peter > > On Wednesday, June 18, 2

Re: [julia-users] Another Julia/Matlab timing comparison

Suppose you're going to be working with a 864 x 729 array of Float32s, and you want the rfft. Then just do this: A = rand(Float32, 864, 729) f = plan_rfft(A, (1,2), FFTW.PATIENT, 20) This will devote up to ~20s testing different algorithms to compute the rfft. Now say B has some data you actual

Re: [julia-users] type mismatch weirdness

Mystery solved thanks to Alireza clever guess (imImage was deffined in two included files) https://groups.google.com/forum/?fromgroups=#!topic/julia-dev/GojOx4nI-xo Terça-feira, 17 de Junho de 2014 17:33:28 UTC+1, J Luis escreveu: > > *imType* is a composite type > > Keno, is right. But now I do

[julia-users] Re: help with debugging julia code

…also you’ve got a spurious : in else: on line 77. On Wednesday, June 18, 2014 4:59:34 PM UTC-7, Ethan Anderes wrote: I also noticed that when D is an array, statements like D(k,i) = d should > be written D[k,i] = d instead. > > On Wednesday, June 18, 2014 4:40:11 PM UTC-7, shashank shekhar wr

[julia-users] Re: help with debugging julia code

I also noticed that when D is an array, statements like D(k,i) = d should be written D[k,i] = d instead. On Wednesday, June 18, 2014 4:40:11 PM UTC-7, shashank shekhar wrote: Hi > > I have written an algorithm for the K Center problem in 1 dimensional > using dp. Here : > https://github.com/

Re: [julia-users] help with debugging julia code

You have a colon after function KCenter1D(y::Array{Float64,1},K::Int64), which gets the parser confused. On Wed, Jun 18, 2014 at 7:40 PM, shashank shekhar wrote: > Hi > > I have written an algorithm for the K Center problem in 1 dimensional > using dp. Here : > https://github.com/shkr/Algorithm

[julia-users] help with debugging julia code

Hi I have written an algorithm for the K Center problem in 1 dimensional using dp. Here : https://github.com/shkr/Algorithms/blob/master/KCenter1D.jl But upon including the file *julia> **include("KCenter1D.jl")* *ERROR: syntax: malformed function arguments 1* * in include at boot.jl:238* *

Re: [julia-users] Another Julia/Matlab timing comparison

Thanks Patrick. Yep, I did see that but I couldn't parse it. I guess they are global variables but I'm not sure how to set them or use them. Sorry to be dense here... On Wednesday, June 18, 2014 2:54:31 PM UTC-7, Patrick O'Leary wrote: > > The flags are defined in the FFTW documentation: > http

[julia-users] Re: uncurried functions are ok... if you know what you are doing

How about pairs = [e for e in enumerate(sqr)] filter(p -> +(p...)%4 != 0, pairs) ? --Peter On Wednesday, June 18, 2014 3:34:24 PM UTC-7, gentlebeldin wrote: > > When you come across Julia, and know (a bit of) Haskell, comparisons are > inevitable. Haskell has only functions with one argument,

Re: [julia-users] Pass Keyword Argument Dictionary

Leah, this works great, thank you for the answer! On Wed, Jun 18, 2014 at 5:46 PM, Stefan Karpinski wrote: > This solution is what I've used in these situations, but it's not really > satisfactory – it's hard to execute efficiently and it obscures issues with > incorrect keyword usage. It wou

[julia-users] uncurried functions are ok... if you know what you are doing

When you come across Julia, and know (a bit of) Haskell, comparisons are inevitable. Haskell has only functions with one argument, but the function value may be another function. It takes some time to wrap your brain around that idea, but hey, it works. Another type of functions with just one ar

Re: [julia-users] When does inlining occur?

For easy linking, that's this: https://github.com/JuliaLang/julia/blob/5c01387135ba9376bce9c86d115655ec9342e1b8/base/inference.jl#L2455-L2469 It's small but a little hard to interpret without a bit more context (hint, hint, Jameson...). On Wed, Jun 18, 2014 at 6:06 PM, Jameson Nash wrote: > T

Re: [julia-users] When does inlining occur?

The actual set of heuristics is pretty small, and is implemented in inline_worthy On Wednesday, June 18, 2014, Stefan Karpinski wrote: > This is a fairly long, involved chunk of code that applies a lot of > heuristics: > > > https://github.com/JuliaLang/julia/blob/5c01387135ba9376bce9c86d115655e

Re: [julia-users] Another Julia/Matlab timing comparison

The flags are defined in the FFTW documentation: http://www.fftw.org/doc/Planner-Flags.html On Wednesday, June 18, 2014 4:30:30 PM UTC-5, Ethan Anderes wrote: > > Hi Tim: > > Do you happen to know where I can find some documentation or code snippets > on using FFTW.MEASURE or FFTW.PATIENT? I als

Re: [julia-users] Self-referential type array

Right, if that's what you want it to do. Although at that point the name is weird. On Wed, Jun 18, 2014 at 5:32 PM, Luke Stagner wrote: > I think an even better definition would be > julia> type SelfRef >objs::Array{SelfRef,1} >function SelfRef() >x = new

Re: [julia-users] Pass Keyword Argument Dictionary

This solution is what I've used in these situations, but it's not really satisfactory – it's hard to execute efficiently and it obscures issues with incorrect keyword usage. It would be good to come up with a better way to support this. On Wed, Jun 18, 2014 at 5:43 PM, Leah Hanson wrote: > Usin

Re: [julia-users] Pass Keyword Argument Dictionary

Using the same print_args, here's a revised arg_collector: julia> function arg_collector(;kwargs...) # do some operation on args print_args(;kwargs...) #always print args end arg_collector (generic function with 1 method) julia> arg_collector(x=5, y=6, z=7) x 5

Re: [julia-users] When does inlining occur?

This is a fairly long, involved chunk of code that applies a lot of heuristics: https://github.com/JuliaLang/julia/blob/5c01387135ba9376bce9c86d115655ec9342e1b8/base/inference.jl#L1978-L2663 There's a fairly old issue open about adding an inline keyword: https://github.com/JuliaLang/julia/issues

[julia-users] Pass Keyword Argument Dictionary

I have the following code that takes some keywords and does some operations. This is for a wrapper for Gcode (talking to 3D printers), so I'd like to specify an arbitrary axis to control. Below is what I would like to do. function print_args(;args...) > # do some formatting > str = "" >

Re: [julia-users] Self-referential type array

I think an even better definition would be julia> type SelfRef objs::Array{SelfRef,1} function SelfRef() x = new() x.objs = SelfRef[] return x end end julia> x=SelfRef() SelfRef([]) julia> push!(x.objs,SelfRef())

Re: [julia-users] Another Julia/Matlab timing comparison

Hi Tim: Do you happen to know where I can find some documentation or code snippets on using FFTW.MEASURE or FFTW.PATIENT? I also have some performance critical code that I've failed to speed up with plan_fft and plan_ifft. After reading your suggesting I tried to search for info on FFTW.MEASURE

Re: [julia-users] Self-referential type array

Right, this definition is better: julia> type SelfRef objs::Array{SelfRef,1} function SelfRef() x = new() x.objs = [x] return x end end julia> x = SelfRef() SelfRef([SelfRef(#= circular reference =#)]) We don't

Re: [julia-users] Self-referential type array

Thanks Stefan, although I feel like julia> x = SelfRef() 1-element Array{SelfRef,1}: SelfRef([SelfRef(#= circular reference =#)]) Should not return an array of type SelfRef. It should return just an instance of itself As it stands, to get the behavior I want I need to call it like julia> x=SelfR

Re: [julia-users] When does inlining occur?

I think the code that controls inlining is mostly here: https://github.com/JuliaLang/julia/blob/master/base/inference.jl But you'll want to get a response from someone who really understands it. -- John On Jun 18, 2014, at 1:30 PM, Sam L wrote: > Apologies if this is in the manual. I must ha

Re: [julia-users] Keeping cloned Packages up to date

Thanks, I was confused by that but now I have checked and it works! Nice! The reason I need this is to get all the OpenGL packages up to date that Simon is working on (they are not yet registered). Am Mittwoch, 18. Juni 2014 22:32:25 UTC+2 schrieb Stefan Karpinski: > > One caveat is that if a pac

Re: [julia-users] Keeping cloned Packages up to date

One caveat is that if a package cannot be fast-forwarded then it won't be updated. On Tue, Jun 17, 2014 at 5:01 PM, Tobias Knopp wrote: > Ok I thought this would be not the case. But one issue might be that I had > modified files in some packages and this cannot work then of course. > > Am Dien

[julia-users] When does inlining occur?

Apologies if this is in the manual. I must have missed it or it was added since I read through it thoroughly. Under what conditions do functions get inlined by the compiler? I would guess that this does not have a simple answer. Are there any rules of thumb or ways to reason about when inlining

Re: [julia-users] MethodError(getindex,(sqrt,837))

On Wednesday, June 18, 2014 03:26:07 PM Jacob Quinn wrote: > This however will throw a DomainError, since the argument is negative. You > may have to clarify what you're actually after here. But see http://docs.julialang.org/en/latest/manual/faq/#why-does-julia-give-a-domainerror-for-perfectly-se

Re: [julia-users] Re: Project organization and CLI

It's not uncommon in Julia to have stateless "empty" types just for dispatch. The real question here is whether the "report" verb really means the same thing or not for n-grams and words. If they're different variations on the same meaning, then the two methods should belong to the same generic fun

Re: [julia-users] MethodError(getindex,(sqrt,837))

To call a function, you use parenthesis, not square brackets, so you want sqrt(-3) This however will throw a DomainError, since the argument is negative. You may have to clarify what you're actually after here. -Jacob On Wed, Jun 18, 2014 at 1:14 PM, alex ortin wrote: > sqrt[-3] > >

Re: [julia-users] Self-referential type array

This works: julia> type SelfRef objs::Array{SelfRef,1} function SelfRef() x = new() x.objs = [x] end end julia> x = SelfRef() 1-element Array{SelfRef,1}: SelfRef([SelfRef(#= circular reference =#)]) A couple of points: The

Re: [julia-users] Self-referential type array

See: https://groups.google.com/forum/#!msg/julia-users/sv5UpxA79zQ/ot-9goKR554J On Wed, Jun 18, 2014 at 4:32 AM, Luke Stagner wrote: > Hello all, > > I know it is possible to have a self-referential type like the following > > type SelfRef >obj::SelfRef >SelfRef() = (x=new(); x.obj=x) >

Re: [julia-users] Self-referential type array

obj::Array{SelfRef,0} # The 2nd parameter to Array is the dimension, so you're declaring obj to be a zero-dimensional array; you probably want obj::Array{SelfRef,1} Similarly with: x.obj =Array(x,0) the 1st argument to Array is the *type* of Array you want to create, so this should be x.obj=Ar

[julia-users] Self-referential type array

Hello all, I know it is possible to have a self-referential type like the following type SelfRef obj::SelfRef SelfRef() = (x=new(); x.obj=x) end what I want to do is have obj be an array of type SelfRef so I can make a tree-like structure. I tried the following type SelfRef obj::Array

[julia-users] MethodError(getindex,(sqrt,837))

I've just started to program with julia and i'm trying use the sequence of cubic root and i got this message but i don't find documentation , how solve it. Left the copy of the program: thanks if (length(p)==4) println ("Polinomio Grado TRES\n") a=p[1]

Re: [julia-users] Re: default values and default constructor for Julia variables

Hi Steve, sorry the delay, I was out... I think the dictionary questions were answered by Kevin. > With regard to the inner and outer constructor both needed, can you give me > a reference? I have read the manual section you mentioned; the constructor > for the Point example in the manual take

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

Great lesson thanks, I'll try :) Paul W dniu 2014-06-18 16:51, Tim Holy pisze: On Wednesday, June 18, 2014 06:06:24 AM paul analyst wrote: >Thanks, but I'm an analyst market, not a programmer;/ Still, I really like >Julia. Historically, only SPSS. I can so deeply programmed? Any tips? Yes, wit

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

On Wednesday, June 18, 2014 06:06:24 AM paul analyst wrote: > Thanks, but I'm an analyst market, not a programmer;/ Still, I really like > Julia. Historically, only SPSS. I can so deeply programmed? Any tips? Yes, with Julia it's easy to reach in "deep." First you might want to look at existing m

Re: [julia-users] default values and default constructor for Julia variables

Hi Steve, I saw your original post, too, and started to respond, but got distracted and it fell off my radar. Sorry about that. Mayor gave good answers to your questions. To clarify one point, I'm in the process of adding functions for sorting OrderedDicts, which are normally in insertion order,

[julia-users] Re: Type variables acting differently when accessed through an array

Issue filed: https://github.com/JuliaLang/julia/issues/7302 On Tuesday, June 17, 2014 2:37:59 PM UTC-4, Jack Holland wrote: > > Does anyone have an explanation for why a type variable (e.g. Uint64) acts > differently when accessed in an array (e.g. [Uint64]) than by itself? > Here's some sample

[julia-users] Re: Type variables acting differently when accessed through an array

Notice also that convert(t, "4") just returns a SubString{ASCIIString} and not a Uint64. On Tuesday, June 17, 2014 7:37:59 PM UTC+1, Jack Holland wrote: > > Does anyone have an explanation for why a type variable (e.g. Uint64) acts > differently when accessed in an array (e.g. [Uint64]) than by

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

Thanks, but I'm an analyst market, not a programmer;/ Still, I really like Julia. Historically, only SPSS. I can so deeply programmed? Any tips? Paul W dniu środa, 18 czerwca 2014 14:51:03 UTC+2 użytkownik Tim Holy napisał: > > If you left them as ranges rather than converting them to vectors (if

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

If you left them as ranges rather than converting them to vectors (if that fits your usual usage case), it would be possible to have a special version of setdiff operating on ranges that would be much faster. Would be great if you could write & submit that. --Tim On Wednesday, June 18, 2014 02

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

Ok, txh, A have lot of long vectors. Is it good time or look for something faster? julia> a=[1:1:10^6]; julia> b=[1:10:10^6]; julia> @time setdiff(a,b) elapsed time: 1.339431633 seconds (252659008 bytes allocated, 14.36% gc time) 90-element Array{Int64,1}: Paul W dniu 2014-06-18 13:38,

[julia-users] Re: default values and default constructor for Julia variables

Mauro, Thanks for taking the time to respond. First, with regard to the existing Dict in Julia, just to be sure I understand, the order of the keys is not the order given by the '<' operation acting on the keys but is rather the order in which they were inserted, is that correct? With regard

Re: [julia-users] Re: `for i in ${SCALAR}`

On 18 June 2014 19:52, Tim Holy wrote: > > Pontus, if your function call is sometimes returning a range and sometimes > returning an integer, then you may have a type-stability issue anyway. To make > it type-stable, you may want to have it return 17:17. I think I formulated myself rather badly,

Re: [julia-users] Re: Array "chunks"

On Tuesday, June 17, 2014 10:35:18 PM UTC-4, Kevin Squire wrote: > > (It's not clear to me how your answer matches the OP's request...?) > > The partition function in https://github.com/JuliaLang/Iterators.jl will > do this: > Ah, partition. Thank you, Mr. Squire!

Re: [julia-users] Another Julia/Matlab timing comparison

Thanks Tim, Yes, profiling suggest that ffts are the main contributor. I remember the discussion (https://groups.google.com/d/msg/julia-users/kV-sKVFR2n0/X9F86aVKWDIJ), and I suppose we could come to the same conclusions, however I could be that I've overlooked something elementary. I was post

Re: [julia-users] default values and default constructor for Julia variables

Kevin is working on a sorted dict as we speak: https://github.com/JuliaLang/DataStructures.jl/pull/43 On Wed, 2014-06-18 at 13:03, mauro...@runbox.com wrote: > I re-read your original post. There are iterators available for > Associative KeyIterator{T<:Associative}, ValueIterator{T<:Associative}

Re: [julia-users] default values and default constructor for Julia variables

I re-read your original post. There are iterators available for Associative KeyIterator{T<:Associative}, ValueIterator{T<:Associative} and the usual iterator which returns key and value as a tuple. These are all also available for the dicts in DataStructures.jl. So you can do for (key, val) in

[julia-users] default values and default constructor for Julia variables

Dear Julia colleagues, Earlier I asked whether there is a datatype similar to "map" from C++, and nobody answered, so I decided as my first nontrivial Julia program to write such a type. This is an associative array (i.e., a dictionary), with the added property that the keys have a sort-order,

Re: [julia-users] Another Julia/Matlab timing comparison

Have you profiled it? Is most of the time spent in the ffts? Assuming it is the ffts that dominate, try preplanning with FFTW.MEASURE or FFTW.PATIENT. I seem to recall that some months ago I posted some Matlab/Julia FFT benchmarks to one of the mailing lists or the github repository. I suspect

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

setdiff --Tim On Wednesday, June 18, 2014 01:29:05 PM Paul Analyst wrote: > number removal line is dynamic. This time 4 then maybe 2 maybe 3 etc > > Question now : > How to construct fromvector [1:10] new vector without [2,4] like this > [1,3,5:10] > > Paul > > W dniu 2014-06-18 12:53, Tim Hol

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

number removal line is dynamic. This time 4 then maybe 2 maybe 3 etc Question now : How to construct fromvector [1:10] new vector without [2,4] like this [1,3,5:10] Paul W dniu 2014-06-18 12:53, Tim Holy pisze: A = rand(5,5) A = A[[1:2,3:5], :] --Tim On Wednesday, June 18, 2014 11:26:36 A

[julia-users] Re: Benchmarking study: C++ < Fortran < Numba < Julia < Java < Matlab < the rest

Good to know. Looks like that we should consider improve the performance of openlibm, especially for commonly used functions such as exp, log, etc. On Tuesday, June 17, 2014 7:23:37 PM UTC-5, Alireza Nejati wrote: > > Dahua: On my setup, most of the time is spent in the log function. > > On Tues

Re: [julia-users] Type variables acting differently when accessed through an array

Seems like a bug to me. Can you file as an issue? --Tim On Tuesday, June 17, 2014 11:37:59 AM Jack Holland wrote: > Does anyone have an explanation for why a type variable (e.g. Uint64) acts > differently when accessed in an array (e.g. [Uint64]) than by itself? > Here's some sample code to show

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

A = rand(5,5) A = A[[1:2,3:5], :] --Tim On Wednesday, June 18, 2014 11:26:36 AM Paul Analyst wrote: > Is possible or not ? Delete rows in array ? > Paul > > W dniu 2014-06-17 07:49, Paul Analyst pisze: > > Thx, for this info > > But not about items I need delete some rows ... > > > > How fast d

Re: [julia-users] Re: `for i in ${SCALAR}`

Pontus, if your function call is sometimes returning a range and sometimes returning an integer, then you may have a type-stability issue anyway. To make it type-stable, you may want to have it return 17:17. Regarding integers-as-iterables, I think the issue is that it makes it much easier to

[julia-users] Re: Another Julia/Matlab timing comparison

And versioninfo(): Julia Version 0.3.0-prerelease+3144 Commit 6ae8670* (2014-05-21 21:24 UTC) Platform Info: System: Darwin (x86_64-apple-darwin12.5.0) CPU: Intel(R) Core(TM)2 Duo CPU P8600 @ 2.40GHz WORD_SIZE: 64 BLAS: libgfortblas LAPACK: liblapack LIBM: libopenlibm Den onsdag d

[julia-users] Another Julia/Matlab timing comparison

Hello all, I'm computing a 2D convolution with fft and seeing a timing difference between Matlab and Julia, with Matlab being about 2x faster than Julia. I've simplified my code down to an example (see gist here: https://gist.github.com/1oly/51f39a831cef3931e7b8). The bottleneck of my code is

Re: [julia-users] Re: How Delete row or col of array? like deleteat! in vectors

Is possible or not ? Delete rows in array ? Paul W dniu 2014-06-17 07:49, Paul Analyst pisze: Thx, for this info But not about items I need delete some rows ... How fast delete rows in arrray, one [3,:]or more [[2,4],:] ? Paul W dniu 2014-06-16 21:22, Stefan Schwarz pisze: deleting in place a

[julia-users] Re: `for i in ${SCALAR}`

I don't think this is a parser issue. It would be strange to not allow integers specifically in a context where we accept anything else. The problem is how integers are defined as iterable in base/number.jl#L36-l39 , I'm su

[julia-users] `for i in ${SCALAR}`

Everyone, I am sure that there is a perfectly good reason for the parser to allow a construct like: for i in 17; println(i); end But, I am unable to see how this can be useful/helpful since a loop like this will only ever execute once. I have now had two bugs stemming from this behaviour, f