Sorry I forgot to add:
JuliaIO/Images.jl relies on having ImageMagick installed, whereas
RawArray.jl is a pure Julia solution without any dependencies.
On Monday, September 26, 2016 at 8:59:15 AM UTC-5, David Smith wrote:
>
> Hi, Isaiah. This is a valid question.
>
> 0. As a
rge.net/nrrd/descformat.html
>
> On Sunday, September 25, 2016, David Smith <david...@gmail.com
> > wrote:
>
>> Hi, all:
>>
>> I finally pushed this out, and it might satisfy some of your needs for a
>> simple way to store N-d arrays to disk. Hope you enj
Hi, all:
I finally pushed this out, and it might satisfy some of your needs for a
simple way to store N-d arrays to disk. Hope you enjoy it.
RawArray (.ra) is a simple file format for storing n-dimensional arrays.
RawArray was designed to be portable, fast, storage efficient, and future
t; @everywhere using MyModule
>
> It seems odd to call using twice but it seems to work.
>
> Nonetheless, I also would be interested in alternative methods that do not
> require modifying LOAD_PATH.
>
>
>
> On Friday, October 16, 2015 at 4:07:49 PM UTC-4, David Smith w
I'm also wondering how to replace 'require' without modifying the
LOAD_PATH. I can't expect my users to know how to make any path changes.
'using' and 'import' doesn't work for my parallel processes even though I'm
loading a registered package ("DCEMRI"). I get this error.
WARNING: Module
:836
in anonymous at task.jl:63
Worker 2 terminated.
It doesn't seem to be quite the same problem as #9245, but rather the
entire module is not able to be found.
On Friday, October 16, 2015 at 3:43:20 PM UTC-5, Tim Holy wrote:
>
> On Friday, October 16, 2015 01:07:49 PM David Smith
Thank you, this looks very nice.
On Sunday, September 13, 2015 at 4:29:15 AM UTC-5, Diego Javier Zea wrote:
>
> I was trying to save a matrix of *88000 x 88000* distances on a *Symmetric
> *but I got an *OutOfMemoryError*. Trying to use a sparse matrix solves
> the memory problem, but takes
I would buy stickers. I need to replace my poetically faded and secretly
embarrassing Matlab bumper sticker.
On Wednesday, March 4, 2015 at 8:32:15 PM UTC-6, Sorami Hisamoto wrote:
It would be nice if Julia T-shirt is on sale to fund JuliaCon 2015. I
guess there are a certain number of
I made this helper function a while back to display images using Gadfly:
function imshow(x, title::String=, units::String=, args...)
is, js, values = findnz(x)
m, n = size(x)
df = DataFrames.DataFrame(i=is, j=js, value=values)
plot(df, x=j, y=i, color=value,
...@gmail.com
javascript: wrote:
Have you tried Winston?
https://github.com/nolta/Winston.jl
On Tuesday, December 16, 2014 11:15:42 PM UTC+2, David Smith wrote:
Thank you.
I feel like Julia has matured enough finally to start migrating all of
my MRI research over to it. So far I have found
:
On Tuesday, December 16, 2014 4:15:42 PM UTC-5, David Smith wrote:
My biggest wish-list item (as a medical imager) would be native Julia
plotting that is similar to Matplotlib. I'd rather not have to require
that people have Python alongside Julia. Makes Julia sound less mature.
Julia
I'm very interested in this too.
On Wednesday, December 17, 2014 4:21:44 PM UTC-6, Seth wrote:
I'm wondering whether folks are actually basing their repos in ~/.julia
and doing their development and commits from there, or whether there's some
other process that allows development to
A few of us around here do medical imaging research, so I'm announcing the
release of DCEMRI.jl, a Julia module for processing dynamic contrast
enhanced magnetic resonance imaging (MRI) data.
http://github.com/davidssmith/DCEMRI.jl
To install,
julia Pkg.add(DCEMRI)
To run a quick demo,
.
On Tue, Dec 16, 2014 at 1:50 PM, David Smith david...@gmail.com
javascript: wrote:
A few of us around here do medical imaging research, so I'm announcing
the release of DCEMRI.jl, a Julia module for processing dynamic contrast
enhanced magnetic resonance imaging (MRI) data.
http://github.com
+1 for a runnable code widget.
I think the typical scientific user will want to see immediately how to
1. create some variables.
2. do some math, probably with a matrix or vector.
3. plot something.
4. make a function.
That is a basic fooling-around paradigm. Having a plot in the browser
Hear hear!
I hope the sheer number of contributors back to the language is evidence of
how appreciative we all are to have Julia. And besides the language being
great, the community is really extraordinary.
Beers are on me if you guys find yourselves in Nashville.
On Monday, December 8,
This is the code:
using PyPlot
x = rand(8,8)
x[1] = NaN
pcolor(x)
On Wednesday, November 26, 2014 8:05:32 AM UTC-6, Steven G. Johnson wrote:
On Wednesday, November 26, 2014 12:45:09 AM UTC-5, David Smith wrote:
My up-to-date version of Matplotlib won't plot NaNs either.
That's weird; it works
Sorry, forgot the version. That was
matplotlib1.4.0np19py27_0
Also consider using for loops to eliminate the need for storage of
intermediates. For loops are much faster in Julia than in Matlab.
On Tuesday, November 25, 2014 2:26:39 PM UTC-6, Peter Simon wrote:
You can find examples of both meshgrid and ndgrid implemented in
My up-to-date version of Matplotlib won't plot NaNs either.
On Tuesday, November 25, 2014 10:22:22 PM UTC-6, Steven G. Johnson wrote:
(Note that the example code for the Mandelbrot set above produces an array
containing NaN values. Maybe that doesn't work if you have an old version
of
Some ideas:
Is there a way to return an error for accesses before at least one
assignment in bits types? I.e. when the object is created uninitialized it
is marked dirty and only after assignment of some user values can it be
cleanly accessed?
Can Julia provide a thin memory management layer
To add some data to this conversation, I just timed allocating a billion
Int64s on my macbook, and I got this (I ran these multiple times before
this and got similar timings):
julia N=1_000_000_000
10
julia @time x = Array(Int64,N);
elapsed time: 0.022577671 seconds (800128
(814320 bytes allocated, 0.22% gc
time)
On Mon Nov 24 2014 at 3:18:39 PM Erik Schnetter schn...@cct.lsu.edu
javascript: wrote:
On Mon, Nov 24, 2014 at 3:01 PM, David Smith david...@gmail.com
javascript: wrote:
To add some data to this conversation, I just timed allocating a billion
Int64s
zero twice) to the cost of not
initializing the memory before writing to it. I could alternatively have
done: `@time x=zeros(); @time fill(x, 0)` to measure the same information.
On Mon Nov 24 2014 at 5:57:29 PM David Smith david...@gmail.com
javascript: wrote:
Did you mean to call zeros
measured the speed of 1+2+3 vs 1+2+3+3.
On Mon Nov 24 2014 at 6:59:50 PM David Smith david...@gmail.com
javascript: wrote:
But you initialized it in both cases. Is there a compiler optimization
going on here that combines the zeros() and fill()?
On Monday, November 24, 2014 5:12:56 PM UTC-6
The following code works as posted, but when I increase N to 1000, I
get a segfault. I don't see why that would break anything, as the size_t
type should be large enough to handle that.
Can someone point me in the right direction?
Thanks!
#include julia.h
#include stdio.h
int main()
nothing about gc.)
On Friday, November 21, 2014 12:08:44 PM UTC-7, Jameson wrote:
You are missing a gc root for x when you allocate y. I can't say if that
is the only issue, since you seem to be unboxing a jl_array_t* using
jl_unbox_voidpointer
On Fri, Nov 21, 2014 at 1:51 PM David Smith david
Thanks, Mauro. That's exactly what I needed.
On Wednesday, November 19, 2014 2:17:45 AM UTC-7, Mauro wrote:
immutable Foo{T,M,N} # probably want an immutable here, so nobody swaps
# out the arrays behind your back.
x::Array{T,M}
y::Array{T,N}
I would like to have a type such as this:
type foo{T,M,N}
x::Array{T,M}
y::Array{T,N}
end
Where T and M can be anything, but N = M + 1. I know I can't write
type foo{T,M}
x::Array{T,M}
y::Array{T,M+1}
end
So I wrote the following naive code, but I get an error:
julia
In case this helps, I wrote my own imshow() function a while back that
copied spy() but took varargs and had more customizations explicitly
included. It might fix the label problems or at least give you inspiration.
function imshow(x, title::String=, units::String=, args...)
is, js, values
Time is a flat circle.
On Monday, November 10, 2014 2:24:09 AM UTC-6, John Myles White wrote:
Yes, the use of zero is an anachronism from a design in which zero was
used to have a default value for arbitrary types.
-- John
On Nov 10, 2014, at 8:22 AM, Ivar Nesje iva...@gmail.com
+1 That would be slick.
On Wednesday, October 15, 2014 8:01:17 AM UTC-5, Erik Schnetter wrote:
Ooh. ∀ i ∈ R... a new way to write for loops!
/performance-tips/#avoid-global-variables
On Friday, October 10, 2014 10:56:36 PM UTC-5, David Smith wrote:
This code runs in 0.5 sec in v0.3.2, but takes 0.64 s in v0.4:
tic()
N = 256
n = 80
x = rand(N,N) + 1im*randn(N,N)
f = zeros(Complex128, N, N, n)
for t = 1:n
f[:,:,t
This code runs in 0.5 sec in v0.3.2, but takes 0.64 s in v0.4:
tic()
N = 256
n = 80
x = rand(N,N) + 1im*randn(N,N)
f = zeros(Complex128, N, N, n)
for t = 1:n
f[:,:,t] = fft(x + float(t))
end
toc()
If I enclose it in a function, as in the following, and then run it by
calling the
Great article! I used some of your tips on a complicated inner loop just
now and immediately got a 3.6x speedup.
On Monday, September 15, 2014 7:57:17 PM UTC-5, Arch Robison wrote:
Thanks to all for taking time to report my errors. I missed @code_llvm's
introduction into Julia. It let
What would be the best way to enclose a function inside of a module that
fires up additional worker processes, where each worker process also sees
the module?
Normally, if starting julia with `julia -p N`, I'd use an `@everywhere
using MyModule`, but this doesn't work here because the other
Is there a preferred idiom in Julia to match the equivalent Python:
if __name__ == '__main__':
run_like_a_command()
else:
be_a_module()
I tried searching, but apparently am not hitting the right keywords. Maybe
my terminology is incorrect.
Thanks, and sorry if this is obvious.
Ha ha ha.
Thank you. Clearly I need more coffee today.
On Wednesday, September 3, 2014 12:32:56 PM UTC-5, Ivar Nesje wrote:
May I suggest searching for __main__
https://groups.google.com/forum/?hl=en#!searchin/julia-users/__main__?
On Wednesday, September 3, 2014 7:20:46 PM UTC+2, David
Add me to the list of people excited about slicing changes.
Is there any way '|' could be distinguished from the boolean operation? As
in, a|1,:,:| ?
What about a[k,.,.]? The dots would remind me that it handles singleton
dimensions in a special way.
Ooh, what about a[k,\cdot,\cdot]?
.
On Wednesday, September 3, 2014 1:21:56 PM UTC-5, Tim Holy wrote:
On Wednesday, September 03, 2014 10:57:16 AM David Smith wrote:
Also, why don't we have a squeeze(A::AbstractArray{T,N}) method that
eliminates all singleton dimensions? If we had that, indexing returning
singleton dimensions
That does not do the same thing as squeeze on a 3+-dimensional array.
On Wednesday, September 3, 2014 1:40:53 PM UTC-5, John Myles White wrote:
Use vec
On Sep 3, 2014, at 11:35 AM, David Smith david...@gmail.com javascript:
wrote:
Ok, so you can continue using the old squeeze. Us
(x, 1) all the time.
On Wednesday, September 3, 2014 1:21:56 PM UTC-5, Tim Holy wrote:
On Wednesday, September 03, 2014 10:57:16 AM David Smith wrote:
Also, why don't we have a squeeze(A::AbstractArray{T,N}) method that
eliminates all singleton dimensions? If we had
Thanks! I'm a very happy user of Juno+LT.
I really enjoy the lesser-known key combos of C-d and C-m, which pop up
help() and methods() for functions.
On Wednesday, September 3, 2014 3:58:54 PM UTC-5, Mike Innes wrote:
Hey all,
Following IPython's rename to Jupyter I've updated the name
Ok, the Ctrl-d doc shortcut is broken now. I get this error. Do I still
need the Jewel julia package?
I get a similar error with Ctrl-m as well.
-
WARNING: LightTable.jl: `getthing` has no method matching getthing(::Nothing)
in anonymous at
Pkg.update()?
On 3 September 2014 17:21, David Smith david...@gmail.com javascript:
wrote:
Ok, the Ctrl-d doc shortcut is broken now. I get this error. Do I still
need the Jewel julia package?
I get a similar error with Ctrl-m as well.
-
WARNING: LightTable.jl: `getthing` has
Can anyone tell my why the code below fails?
Why is the error reporting an incorrect size for the right-hand side?
Is this a gotcha of the broadcast rules?
Is the only solution to boost S to a 1 x n x n array?
I prefer to avoid looping, if possible.
Thanks!
function f(n)
A =
Probably a dumb question, but I can't find what I'm wanting.
How do I turn this signature into one that is generic for any x that is an
Array of Float or Complex of any dimension?
func(x::Array{Complex{Float64},1}, q::Float64)
Thanks!
47 matches
Mail list logo