Re: [julia-users] Handling non-convergence in roots.jl

Thanks! On Thu, Nov 10, 2016 at 6:03 PM, Mauro wrote: > > On Thu, 2016-11-10 at 20:07, Mauro wrote: > > On Thu, 2016-11-10 at 18:54, Joaquim Masset Lacombe Dias Garcia < > joaquimdgar...@gmail.com> wrote: > >> I am running a sequence of root finding functions in a

[julia-users] Handling non-convergence in roots.jl

I am running a sequence of root finding functions in a loop. For my problem the basic fzero always work. However, using Newtons method is very useful because it is very fast when it converges. The problem is that when it does not converge it throws an error instead of a non-convergence flag. Th

[julia-users] Re: jump ahead in Mersenne Twister

Thanks Dan, I will try using the dSFMTjump lib! best,

[julia-users] jump ahead in Mersenne Twister

Hi, Have anybody implemented an arbitrary jump ahead for Mersenne Twister in julia? I saw that the randjump documentation that goes as following: randjump(*r::MersenneTwister*, *jumps*[, *jumppoly*]) → Vector{MersenneTwister} Create an array of the size jumps of initialized MersenneTwister RN

[julia-users] Re: ANN: ParallelDataTransfer.jl

Very nice! Em segunda-feira, 3 de outubro de 2016 11:44:04 UTC-3, Andre Bieler escreveu: > > Thank you Sir! :) > > On Sunday, October 2, 2016 at 5:30:34 PM UTC+2, Chris Rackauckas wrote: >> >> ParallelDataTransfer.jl >> is a >> librar

[julia-users] MPI.jl and composite types

Hi all, I have the following type: type bar a::Int b::Vector{Matrix{Float64}} c::Vector{Float64} end I would like to send instances of that type via MPI. Can I do it? Is there any serialization/deserialization procedure i could use for that? thanks!

[julia-users] Re: Int Matrix multiplication versus Float Matrix multiplication

Thanks again! I just found this: http://stackoverflow.com/questions/2550281/floating-point-vs-integer-calculations-on-modern-hardware According to the information presented there, generic multiplication is slower for ints than for floats. Em quarta-feira, 21 de setembro de 2016 20:14:34 UTC-3,

[julia-users] Int Matrix multiplication versus Float Matrix multiplication

Does anyone knows why the float matrix multiplication is much faster in the following code: function foo() n = 2000 a = rand(1:10, n, n) b = rand(n, n) @time a*a @time b*b nothing end I get: julia> test() 6.715335 seconds (9 allocations: 30.518 MB, 0.00% gc time) 0.120801

[julia-users] Precompilation (using usrimg.jl) of pkg using binary lib

Hi all! I have a package which depends a external binary library(I will just write dll from now on), I want to figure out a way to call the dll in a way that if some user decides to precompile a system image with this library it can still work. My package is Xpress.jl (I know it wraps commerci

[julia-users] Julia with fewer dependencies?

The basic question is "Can I compile a smaller version of julia?" For instance, I want to ship some program as an executable, something similar to what BuildExecutable.jl does, btw its very nice! The only problem is that with BuildExecutable.jl or even with the pure usrimg.jl trick create a hug

Re: [julia-users] Release memory GC

So the recommendation for clearing memory is calling GC after some function so that no references exist anymore? Em quinta-feira, 7 de julho de 2016 17:17:03 UTC-3, Yichao Yu escreveu: > > On Thu, Jul 7, 2016 at 4:10 PM, Thuener Silva > wrote: > > I have a code that uses a lot of memory and I'm

[julia-users] Re: A question about allocations

I noticed that using "const temp = zeros(1)" and then for i in 1:3 @time temp[1]+=foo(a[1],a[2],3)*2+3.0 end results in zero allocations, but is this the best way? Em segunda-feira, 29 de fevereiro de 2016 12:12:20 UTC-3, Joaquim Masset Lacombe Dias Garcia escreveu: > > I ha

[julia-users] A question about allocations

I have the following code, and I was wondering if I could run it with zero allocations. If yes, how shoud I modify it, if no, is there a simple exaplanation? thanks in advance. a = ones(2) function foo(a,b,c) return a*1+c*b*1+3 end temp = 0.0 for i in 1:3 @time temp+=foo(a[1],a[2],3)*2+3.0 en

[julia-users] CUDArt: loop inside device do

Can I have a standard julia "for loop" inside a "device do" of CUDArt? I tried the following example: using CUDArt, MyCudaModule nrow = 10 ncol = 3000 mat = ones(Float64,nrow,ncol) out1 = zeros(Float64,nrow) vec = Float64[1:nrow;] out2 = zeros(Float64,nrow) d_mat = CudaArray(mat) d_out1 = Cud

[julia-users] Why is the order of these loops barely affects the time?

Why is the order of these loops barely affects the time? N=100 x=rand(N,N,N) y=rand(N,N,N) function abc(x::Array{Float64,3},y::Array{Float64,3}) for a=1:N for b=1:N for c=1:N y[a,b,c]=x[a,b,c] end end end end function acb(x::Array{Fl

[julia-users] DataFrames and special characters

I have noticed that DataFrames´ readtable cannot read by default special characters such ac $ % - & and it replaces everything by underscore. Is there any way to make it read those characters just like the R function fread does?

[julia-users] Re: Is build_executable.jl working?

|I am trying to build executable on windows, without success, my last error is: ERROR: chmod: no such file or directory (ENOENT) in chmod at fs.jl:184 does anybody know whats this means? Em terça-feira, 1 de dezembro de 2015 08:40:33 UTC-2, Ján Adamčák escreveu: > > I edit this line > jl_ate

Re: [julia-users] Re: CUDART and CURAND problem on running the same "do" loop twice

And it also works on windows if CUDA version is 7.5! On Wed, Nov 25, 2015 at 7:58 PM, Joaquim Masset Lacombe Dias Garcia < joaquimdgar...@gmail.com> wrote: > Wow! it seems that something is wrong in CUDA 7.0 ! > > I updated my MAC to cuda 7.5 and bith the C code and julia cod

[julia-users] Re: Is build_executable.jl working?

Have you tried this one: https://github.com/dhoegh/BuildExecutable.jl I did NOT have success but if you do you might even help me! this might help https://github.com/JuliaLang/julia/issues/11816#issuecomment-148784543 and this: https://groups.google.com/forum/#!searchin/julia-users/buil_executa

Re: [julia-users] Re: CUDART and CURAND problem on running the same "do" loop twice

; >> It works fine here, on Ubuntu 14.04 and CUDA 6.5. >> >> Cheers, >> [S]. >> >> >> >> >> On Tuesday, November 24, 2015 at 9:30:32 PM UTC-2, Joaquim Masset Lacombe >> Dias Garcia wrote: >>> >>> Interesting, both m

Re: [julia-users] Re: CUDART and CURAND problem on running the same "do" loop twice

Interesting, both my machines (windows and mac) have CUDA 7.0, possibly thats the issue, since the C code in https://github.com/JuliaGPU/CURAND.jl/issues/3#issuecomment-159319580 fails in both. If some linux user could test this version, our statistics would be more complete. I will try 6.5 an

Re: [julia-users] Re: CUDART and CURAND problem on running the same "do" loop twice

if I load: using CUDArt using CURAND and then I run either: d_a = curand(Float64, 1000); a = to_host(d_a); OR: result = devices(dev->capability(dev)[1]>=2) do devlist end I can repeat the block as many time I want. However if I run both combined (or alternatedly in any order): d_a = curand(F

[julia-users] CUDART and CURAND problem on running the same "do" loop twice

I was playing with the GPU (in both windows and mac) libraries and I came up with the following errors: The first time I execute the do loop, every thing goes well. But If i try a second time in a row I get this invalid pointer error. I got the same error in all my do loops, any ideas? best, Joa

Re: [julia-users] CUDArt main example does not run smoothly

> possible for me to fix myself. So if you can figure out what's needed on > your > system, a pull request would be most appreciated. If possible, try to keep > it > as general as possible with respect to CUDA version. > > Best, > --Tim > > On Thursday, Octobe

Re: [julia-users] CUDArt main example does not run smoothly

EDIT of last post, ignore my last 3 lines, the files are in the right place. Em quinta-feira, 15 de outubro de 2015 18:32:04 UTC-3, Joaquim Masset Lacombe Dias Garcia escreveu: > > Thats Great that you could update the README! > > I am now testing the pkg on mac (later I was on win

Re: [julia-users] CUDArt main example does not run smoothly

've updated the README, so you might want to try > again. I also tagged a new version of the package, so you'll have access > to > the newer features described in the README; you'll probably want to do > Pkg.update(). > > Best, > --Tim > > On Wednesday

[julia-users] DataFrames only with header line and no other lines

I have been using the DataFrames package which is excellent by the way! I don't know if its a problem but if i try to read a .csv file with only one line (which is supposed to be the header of the dataframe), it returns the error: LoadError: ArgumentError: Columns and column index must be the s

[julia-users] CUDArt main example does not run smoothly

I could not run the main example of CUDArt from : https://github.com/JuliaGPU/CUDArt.jl The code and error follow bellow. running: #MyCudaModule previously defined as in: https://github.com/JuliaGPU/CUDArt.jl using CUDArt, MyCudaModule A = rand(10,5) result = devices(dev->capability(dev)[1]>=

[julia-users] CUDArt example

Hi, I am trying to use CUDArt library but i could not run the example in the git repository. I ran the pkg.test("CUDArt") and it seems to be well installed. Does anyone has a simple example. Something like summing two vectors with vadd.cu (of CUDA.jl git repository) would be awesome! Thanks!