[julia-users] Multi-threaded BLAS limiting at 16 cores

2015-04-09 Thread Chris Rackauckas
Hello, I am new to Julia and am trying to get it to use multi-threaded BLAS via openBLAS on an HPC node. I compiled it on the HPC and in Julia I check julia versioninfo() Julia Version 0.3.8-pre+2 Commit af187c5 (2015-03-25 14:33 UTC) Platform Info: System: Linux (x86_64-unknown-linux-gnu)

[julia-users] Private Forks of Julia Repositories

2016-06-02 Thread Chris Rackauckas
Hey, I was wondering if anyone has experience creating a private fork of a Julia repository. I would like to implement some research codes via modifying (my own) public repository, keep it private until publication, and then merge it into the main repository. I found a StackOverflow solution

[julia-users] Re: Status of FEM packages

2016-06-08 Thread Chris Rackauckas
You have those packages and DifferentialEquations.jl which implements solvers for Poisson and Heat problems. I am looking to do more (and multithread the vectorized calls it currently has), but right now the landscape is kind of bare. DifferentialEquations.jl will most likely do what's in the

Re: [julia-users] Differential Equations Package

2016-06-08 Thread Chris Rackauckas
Thanks for chiming in! There's a lot to say about this topic, and I think it's great! I felt like Julia's differential equation solvers were behind where they should be, which was crazy to me since there is so many state-of-the-art mathematical packages in other areas (optimization, stats,

[julia-users] juliabloggers.com: Using Julia’s Type System For Hidden Performance Gains

2016-06-08 Thread Chris Rackauckas
It should all be fixed now.

[julia-users] Re: Differential Equations Package

2016-06-06 Thread Chris Rackauckas
Hey, This is the author here. Let me address a few things. Giving the solution is optional for any problem type. So you can just declare a PoissonProblem(f). If you give it a solution, then it calculates errors and a bunch of extra things. I'll clear up that part of the tutorial. You can

Re: [julia-users] Differential Equations Package

2016-06-06 Thread Chris Rackauckas
It also solves ODEs and SDEs, and I will be adding DDEs (delay equations) and differential algebraic equations ASAP. It implements more solvers for ODEs (with more tests) than ODE.jl, and implements more modern solvers for SDEs than SDE.jl. It also implements solvers for common (S)PDEs, and I

[julia-users] Re: Memory corruption when using BigFloats

2016-06-07 Thread Chris Rackauckas
> string( bigpi )# if possible >> >> >> >> >> On Tuesday, June 7, 2016 at 2:05:34 PM UTC-4, Pablo Zubieta wrote: >>> >>> It would help to have all the fields of the BigFloat, not only d (but >>> also prec, sign and exp). I think

[julia-users] Re: Memory corruption when using BigFloats

2016-06-07 Thread Chris Rackauckas
Forgot to mention that if you change the deepcopy in the loop to copy, it works.

Re: [julia-users] Differential Equations Package

2016-06-07 Thread Chris Rackauckas
There is a chance I will be adding some bifurcation analysis soon. One of the projects I'm working on may need to do a bifurcation analysis soon, and when that happens I am going to make some arclength and homotopy continuation bifurcation diagram(mers?). It will take the same ODEProblem type

Re: [julia-users] Differential Equations Package

2016-06-09 Thread Chris Rackauckas
ithub.com/ChrisRackauckas/DifferentialEquations.jl/blob/master/src/ode/ODECoefficientTypes.jl. I just checked and am surprised you guys don't have a DP78 method, but it shouldn't take more than minutes to plop that tableau over there. > A few more comments in-line below. Cheers! Mauro &g

Re: [julia-users] Re: Private Forks of Julia Repositories

2016-06-03 Thread Chris Rackauckas
I can keep a private local branch, but then it's only on one computer and I can't develop/test on any other computer (/hpc). On Thursday, June 2, 2016 at 11:18:09 PM UTC-7, Mauro wrote: > > On Fri, 2016-06-03 at 07:58, Chris Rackauckas <rack...@gmail.com > > wrote: > >

[julia-users] Re: Private Forks of Julia Repositories

2016-06-02 Thread Chris Rackauckas
I think I will need both versions available, since the majority of the work is public, while the private work will tend to sit around longer (i.e. waiting to hear back from reviewers). So I'd want to be able to easily work with the public repository, but basically switch over to a private

Re: [julia-users] Re: Private Forks of Julia Repositories

2016-06-03 Thread Chris Rackauckas
ge is not really compatible with orignal intent, > so I don't have multiple remotes. > (Well I do for other reasons) > > > On Friday, 3 June 2016 14:18:09 UTC+8, Mauro wrote: >> >> On Fri, 2016-06-03 at 07:58, Chris Rackauckas <rack...@gmail.com> wrote: >>

Re: [julia-users] Re: Private Forks of Julia Repositories

2016-06-03 Thread Chris Rackauckas
Make PR on github. > > I hope this helps. > > > On Thursday, June 2, 2016 at 11:18:09 PM UTC-7, Mauro wrote: > >> > >> On Fri, 2016-06-03 at 07:58, Chris Rackauckas <rack...@gmail.com > >> > wrote: > >> > I think I will need bot

[julia-users] Re: I feel that on the syntax level, Julia sacrificed too much elegancy trying to be compatible with textbook math notations

2016-05-24 Thread Chris Rackauckas
While I don't think it's true that numerical computing is only 5% matrix math, most user facing code isn't matrix math. At its core, just about every numerical algorithm is matrix math. Every nonlinearity or connection between terms becomes a matrix, and so every equation is either solving Ax

Re: [julia-users] Lack of an explicit return in Julia, heartache or happiness?

2016-05-24 Thread Chris Rackauckas
I like Stefan's idea. In fact, this is how I've already been using Julia, precisely for scalability. I think that not having an explicit return value in a long function is just unreadable, so for long-form functions I always add return statements at the bottom. At the same time, there's no

[julia-users] Re: Fun Fact: Julia is space sensitive (in some sense)

2016-05-24 Thread Chris Rackauckas
Yes, spaces denote a new column in an array definition, so a[2 +1] is size 2 while a[2 + 1] has the parser not care about the spaces because it tries to complete the expression (like it does with newlines). My favorite case is this piece of code I wrote a few months ago: Du(x) =

[julia-users] Re: What is Julia and what it is not?

2016-05-26 Thread Chris Rackauckas
As mentioned, the goal is for it to be a general purpose language with a scientific focus: Scientific language - It already has that down. It's easy to do mathematics with it, it's fast to prototype algorithms, and the code has very high performance. This is what most scientists are looking

[julia-users] Julia Memory Management

2016-05-26 Thread Chris Rackauckas
I see mentions like this one every once in awhile: "D language is a special case, as it has GC, but it's also optional (as with Julia)" Is GC optional? I thought the only way to discard something from memory was to set it to zero and call garbage control (which then runs the whole garbage

Re: [julia-users] Differential Equations Package

2016-06-10 Thread Chris Rackauckas
> > I agree, keyword-args are a maintenance headache (and also potentially > bad for performance). That was indeed one of the reason to make > Parameters.jl to generate the keyword constructors for the types > automatically. Then use the types instead of keyword functions. > > Concerning

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

2016-06-14 Thread Chris Rackauckas
I agree that more should be backed out to packages. Even as a math person who would probably never Julia without linear algebra, I think that some of the linear algebra should move out to a JuliaMath package, and then part of the standard startup script should be to "using LinAlg". That way if

[julia-users] Re: Memory corruption when using BigFloats

2016-06-07 Thread Chris Rackauckas
That's the default precision. get_bigfloat_precision() 256 On Tuesday, June 7, 2016 at 11:01:05 AM UTC-7, Pablo Zubieta wrote: > > What is the precision of the BigFloat you wrote above? > > On Tuesday, June 7, 2016 at 7:26:02 PM UTC+2, Chris Rackauckas wrote: >> >> Th

Re: [julia-users] Differential Equations Package

2016-06-07 Thread Chris Rackauckas
d, I did the same for myself in a non repo > package. > > Any one last time, your are a hero for getting this going. I am dying to > start using it in my own work as I love diff.eqs :) > > On Tuesday, June 7, 2016 at 9:23:55 AM UTC-7, Chris Rackauckas wrote: >> >>

Re: [julia-users] Julia Memory Management

2016-05-26 Thread Chris Rackauckas
Thanks, that's good to know. Can you explain a little bit about why it would hurt performance? On Thursday, May 26, 2016 at 8:23:34 AM UTC-7, Yichao Yu wrote: > > On Thu, May 26, 2016 at 11:06 AM, Chris Rackauckas <rack...@gmail.com > > wrote: > > I see mentions lik

Re: [julia-users] Julia Memory Management

2016-05-26 Thread Chris Rackauckas
Ahh, that makes sense. Thanks for clearing that up. On Thursday, May 26, 2016 at 8:59:14 AM UTC-7, Yichao Yu wrote: > > On Thu, May 26, 2016 at 11:44 AM, Chris Rackauckas <rack...@gmail.com > > wrote: > > Thanks, that's good to know. Can you explain a little bit about why

[julia-users] Require a package not in METADATA

2016-06-14 Thread Chris Rackauckas
Is there a way to add a package to REQUIRES that's not in METADATA?

[julia-users] Re: Applied NLA course in Julia

2016-06-21 Thread Chris Rackauckas
@Crivelli You can a packaged form of these types of algorithms (like Jacobi iteration) at IterativeSolvers.jl Thanks for these teaching materials in Julia. They look great! On Tuesday, June 21, 2016 at 1:46:58 PM UTC+1, Dawid Crivelli wrote: > > Thanks for the notes, they look interesting. > >

Re: [julia-users] Re: Deepcopy of BigInt and BigFloat (#16999)

2016-06-21 Thread Chris Rackauckas
; Julia's stdlib, please do file issues. We want to get to the point where > *everything* just works. > > On Mon, Jun 20, 2016 at 4:31 PM, Chris Rackauckas <rack...@gmail.com > > wrote: > >> I see. I surely can't wait for native bigs, but there are definitely more >&g

[julia-users] Re: Fast vector element-wise multiplication

2016-06-20 Thread Chris Rackauckas
I think that for medium size (but not large) arrays in v0.5 you may want to use @threads from the threadding branch, and then for really large arrays you may want to use @parallel. But you'd have to test some timings. On Monday, June 20, 2016 at 11:38:15 AM UTC+1, chobb...@gmail.com wrote: > >

[julia-users] Re: Unexpected error on 1st order ODE using ODE package

2016-06-19 Thread Chris Rackauckas
I wouldn't call this a bug, it's a standard Julia thing for a reason. You get an InexactError() because you start with an Int and you do an operation which turns the Int into a Float so the program gets mad at the type instability. You can just change everything to floats, but then you're

[julia-users] Re: pyplot plots in atom Plots pane

2016-06-19 Thread Chris Rackauckas
You can do it via Plots.jl using the PyPlot backend. You just have to switch an environment variable. However, currently not everything syncs up all the time, so it's more of in a trial phase right now. On Friday, June 17, 2016 at 9:14:22 PM UTC+1, Lewis Lehe wrote: > > Hi I am using the

[julia-users] Re: Unexpected error on 1st order ODE using ODE package

2016-06-19 Thread Chris Rackauckas
simple, but I worry being to harsh about > types when it feels unessary. > > > On Sunday, June 19, 2016 at 5:28:39 PM UTC-7, Chris Rackauckas wrote: > >> I wouldn't call this a bug, it's a standard Julia thing for a reason. You >> get an InexactError() because you start

[julia-users] Re: Unexpected error on 1st order ODE using ODE package

2016-06-19 Thread Chris Rackauckas
a dynamic language and not purely static. Suggesting > users, especially new users, to be worry about these situations doesn't > jive with my understanding of Julia. It feels like premature optimization > to me. > > > On Sunday, June 19, 2016 at 6:38:37 PM UTC-7, Chris Racka

[julia-users] Re: Fast vector element-wise multiplication

2016-06-20 Thread Chris Rackauckas
ighted > above still the best choice for doing the job? > > > On Monday, June 20, 2016 at 1:57:25 PM UTC+1, Chris Rackauckas wrote: >> >> I think that for medium size (but not large) arrays in v0.5 you may want >> to use @threads from the threadding branch, and then

Re: [julia-users] Re: Deepcopy of BigInt and BigFloat (#16999)

2016-06-20 Thread Chris Rackauckas
I see. I surely can't wait for native bigs, but there are definitely more pressing issues to work on first. Tim's proof of principle results look really nice. I tried DoubleDouble and ran into some errors. I don't think it worked with all standard library math. I'll probably take another look

Re: [julia-users] Conditional Definition of Functions?

2016-06-16 Thread Chris Rackauckas
n (if you ask for the first one), or you get nothing! This is on Commit 59d1539 (4 days old master). On Thursday, June 16, 2016 at 1:34:34 AM UTC-7, Mauro wrote: > > On Thu, 2016-06-16 at 07:32, Chris Rackauckas <rack...@gmail.com > > wrote: > > There seems to be an in

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

2016-06-17 Thread Chris Rackauckas
Uwe, take a look here . You can use mkdocs build then mkdocs serve to do so (I think serving will actually build the docs too so that might be redundant). Glad to see this package

Re: [julia-users] Intel Xeon Phi support?

2016-01-18 Thread Chris Rackauckas
Any updates on Julia for the Phi? I know that MKL automatic offload works, but am looking for the whole thing. I do research in stochastic dynamical systems (easy to parallelize) and have a 5110 working, so I am ready to test it once Julia has it!

Re: [julia-users] ODE.jl - backwards (negative in time) integration

2016-06-26 Thread Chris Rackauckas
Won't some of the Runge-Kutta solvers be really unstable with a small negative Delta t? What happens when you use a fixed time step method: does it run and diverge? On Friday, June 24, 2016 at 8:50:42 PM UTC+1, Mauro wrote: > > This is a bug, could you file an issue? Thanks! Note that some

Re: [julia-users] ODE.jl - backwards (negative in time) integration

2016-06-26 Thread Chris Rackauckas
? On Monday, June 27, 2016 at 6:15:45 AM UTC+1, Chris Rackauckas wrote: > > Won't some of the Runge-Kutta solvers be really unstable with a small > negative Delta t? What happens when you use a fixed time step method: does > it run and diverge? >

[julia-users] Conditional Definition of Functions?

2016-06-15 Thread Chris Rackauckas
There seems to be an interesting change between v0.4.5 and v0.5. Take a look at the following code: if largeFluctuations function σ(y,t) dσ = zeros(19) dσ[1] = noiseLevel*1.5y[1] dσ[18]= noiseLevel*6y[18] return(dσ) end else function σ(y,t) dσ =

[julia-users] SIMD long inner expression

2016-02-04 Thread Chris Rackauckas
Hi, I am trying to optimize a code which just adds a bunch of things. My first instinct was to unravel the loops and run it as SIMD, like so: dx = 1/400 addprocs(3) imin = -6 jmin = -2 #Some SIMD @time res = @sync @parallel (+) for i = imin:dx:0 tmp = 0 for j=jmin:dx:0 ans = 0 @simd

[julia-users] Calling your own C Function

2016-02-02 Thread Chris Rackauckas
Hey, sorry if this is a really noob question, but how does one use ccall on your own C function? Since the binding to the Monte Carlo VEGAS integration libraries in GSL.jl currently need work (if anyone wants to fix that up / guide me on how to fix it up, I'll give you 10 cookies!), I was

Re: [julia-users] Calling your own C Function

2016-02-03 Thread Chris Rackauckas
Thanks for the link. That explained what's going on. I thought I could just double check by calling from the command line, but that was wrong. When I called main from Julia, it worked and I got back 0 (since that's what main returns. I guess it doesn't get the stdout as well? I can find that

Re: [julia-users] Calling your own C Function

2016-02-03 Thread Chris Rackauckas
Thanks for the replies. I mixed the compile commands from here with the information on the libraries to link

[julia-users] Some blog posts on Julia

2016-02-22 Thread Chris Rackauckas
Hey, I wanted to let everyone know I put a few blog posts out there on Julia. You can find them here: http://www.stochasticlifestyle.com/ I am by no means an expert on Julia, but I put some tutorials for some ways I have been using Julia, like how to use Julia on the XSEDE Comet HPC

[julia-users] Make the computer make a beep?

2016-02-27 Thread Chris Rackauckas
Hey, I was wondering about a quick way to make the computer beep when code is finished. Usually there's a beep() function. The one in Julia at Base.Terminals.beep takes as input a terminal. Where do I find this terminal object? If it matters, I am using Juno with Atom.

Re: [julia-users] Make the computer make a beep?

2016-02-27 Thread Chris Rackauckas
cks, see > http://docs.julialang.org/en/release-0.4/manual/running-external-programs/ > > Am 28.02.2016 um 00:39 schrieb Chris Rackauckas: > > Hey, > > I was wondering about a quick way to make the computer beep when code > > is finished. Usually there's a beep() fu

Re: [julia-users] Some blog posts on Julia

2016-02-23 Thread Chris Rackauckas
better, sometimes.) > > See http://docs.julialang.org/en/stable/manual/performance-tips/ > > --Tim > > On Monday, February 22, 2016 11:34:20 PM Chris Rackauckas wrote: > > Hey, > > I wanted to let everyone know I put a few blog posts out there on > Julia. > &

[julia-users] Re: Google SoC '16: Native Julia solvers for ordinary differential equations, Julia|Quantum>, and Native Julia implementation of Holonomic Functions

2016-03-18 Thread Chris Rackauckas
Taking a quick look at the paper you linked makes big use of symbolic computing algorithms to compute things like a Groebner basis. To get the infrastructure for what's proposed you would need a whole symbolic computing library, which is why they used Mathematica. I think implementing a fast

[julia-users] Re: Your experiences with Julia's parallel computing

2016-04-07 Thread Chris Rackauckas
> > >- What is the main reason that you started using Julia? > > To more quickly develop code for HPCs. I do stochastic numerical methods in developmental biology. When I am making methods, I want to show that it works, get a relatively fast prototype, and move on (I'll let "library

[julia-users] Re: [ANN] Cuba.jl: library for multidimensional numerical integration

2016-04-12 Thread Chris Rackauckas
Nice! I was looking for a Vegas function awhile ago and the GSL.jl one isn't bound correctly yet. This will be a nice addition to the Julia package list. Good job! On Sunday, April 10, 2016 at 1:34:53 PM UTC-7, Mosè Giordano wrote: > > Dear all, > > I am proud to announce Cuba.jl

[julia-users] Re: Tools for creating presentation slides with Julia code?

2016-03-25 Thread Chris Rackauckas
I've used Beamer with LaTeX and a custom program listing for syntax highlighting with success.

[julia-users] Re: Where to publish numerical work using julia

2016-03-21 Thread Chris Rackauckas
Why not just do Journal of Computational Physics or SIAM Journal on Scientific Computing? I recently submitted numerical work done with Julia to the latter and no one cares that it's Julia. In fact, with numerical work you could use whitespace and no one would care: it's all about the

[julia-users] Re: When will juno-atom bundle become available?

2016-03-21 Thread Chris Rackauckas
It's available. https://github.com/JunoLab/atom-julia-client/tree/master/manual Or do you mean a Julia/Atom bundle with the client pre-installed? On Monday, March 21, 2016 at 1:08:45 AM UTC-7, Yao Lu wrote: > > I 'm expecting this bundle. >

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

2016-03-03 Thread Chris Rackauckas
3, 2016 at 10:54:48 AM UTC-8, Yichao Yu wrote: > > On Thu, Mar 3, 2016 at 1:36 PM, Chris Rackauckas <rack...@gmail.com > > wrote: > > Oh, meant the line with the mic environment variable: > > > > printf("%s\n",getenv("MIC_OMP_NUM_THREADS&quo

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

2016-03-03 Thread Chris Rackauckas
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

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

2016-03-03 Thread Chris Rackauckas
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 soluti

Re: [julia-users] Re: Julia is a great idea, but it disqualifies itself for a big portion of its potential consumers

2016-04-03 Thread Chris Rackauckas
I came to this thinking it would be about how Julia pedagogy is a little more difficult than other scripting languages because there are a lot of details in the advanced features which are normally the points of conversation for more experienced users. That's a legitimate concern. However,

[julia-users] Re: Improvements to GPU integration with Julia

2016-04-26 Thread Chris Rackauckas
I think the GPU integration libraries for Julia are already really good if you're using CUDA. CUDArt.jl and Arrayfire.jl work quite well. I don't know too much about the OpenCL side, but I don't tend to have a use for it. On Monday, April 25, 2016 at 6:06:10 PM UTC-7, Michael Jin wrote: > >

[julia-users] Re: GPU capabilities

2016-04-29 Thread Chris Rackauckas
2016 at 4:31:32 PM UTC-4, Chris Rackauckas wrote: >> >> Works great for me. Here's a tutorial where I describe something I did >> on XSEDE's Comet >> <http://www.stochasticlifestyle.com/julia-on-the-hpc-with-gpus/> which >> has Tesla K80s. It works great. I

Re: [julia-users] Suppressing plot windows by default - is this purposeful?

2016-05-23 Thread Chris Rackauckas
But it does plot it, and use *display()* to display it. And I didn't know the REPL doesn't let you edit plots. It works just fine (i.e. like MATLAB) in Juno. On Monday, May 23, 2016 at 7:37:34 AM UTC-7, NotSoRecentConvert wrote: > > I agree with Dave. Having come from working with Matlab myself

Re: [julia-users] Suppressing plot windows by default - is this purposeful?

2016-05-23 Thread Chris Rackauckas
But it does plot it, and use *display()* to display it. And I didn't know the REPL did that itself didn't let you edit plots. It works just fine (i.e. like MATLAB) in Juno. On Monday, May 23, 2016 at 7:37:34 AM UTC-7, NotSoRecentConvert wrote: > > I agree with Dave. Having come from working

Re: [julia-users] Suppressing plot windows by default - is this purposeful?

2016-05-23 Thread Chris Rackauckas
I still think the Julia way matches better to reality. You can go up to a kid and ask "can you draw me a cow?", and they draw it, and later ask "can you show me the cow?" and they'll show it to you. If you're standing right there (i.e. using the iterative REPL at its highest scope) you'll see

Re: [julia-users] Suppressing plot windows by default - is this purposeful?

2016-05-23 Thread Chris Rackauckas
But it is consistent. If you type 2+2 in the REPL you get 4, and if you put the same code in a loop it won't display 4 anymore. This is highlighted in the documentation. Julia just treats plots like any other type. I think it's just a documentation issue. Somehow it should be noted in the

[julia-users] Re: Let's Bridge the IRC and Gitter

2016-05-23 Thread Chris Rackauckas
The bot worked really well today. Answered questions from people on IRC without a hitch. I say we keep it. On Sunday, May 22, 2016 at 10:08:26 PM UTC-7, Lyndon White wrote: > > Alright for purpose of testing, > I have created an instance of gitter-irc-bot >

[julia-users] Re: Parallel computing

2016-05-22 Thread Chris Rackauckas
You need to use @everywhere to have the function defined on the workers. Can we see some code? On Saturday, May 21, 2016 at 6:33:17 PM UTC-7, SHORE SHEN wrote: > > Hi > > Im running a loop, and system shows that only 20% of CPU were used. I then > think that due to the fact that my CPU is 4

Re: [julia-users] Suppressing plot windows by default - is this purposeful?

2016-05-22 Thread Chris Rackauckas
As I noted when I explained this solution to him in the Gitter, I think this should be explained here: http://docs.julialang.org/en/release-0.4/manual/noteworthy-differences/. This is a feature for a good reason. It's the same feature that makes it so that you don't have to add semicolins

[julia-users] Re: Let's Bridge the IRC and Gitter

2016-05-22 Thread Chris Rackauckas
I like this idea. I am fine if it's a different channel: it's easy to sit on different channels on Gitter and get notified to switch between them. I'd IRC but I wouldn't even know a good client to use: I kind of think of it as antiquated and from a previous generation (with no ill will attached

[julia-users] Package Documentation Tools

2016-05-11 Thread Chris Rackauckas
Hi, I was wondering if there's any documentation/tutorials for generating documentation for Julia packages. I would like to make one of those Read the Docs things but I don't know where to start (or if that's still the preferred method) and a quick Google / Julia-users search didn't hit a

[julia-users] Re: Package Documentation Tools

2016-05-11 Thread Chris Rackauckas
jl/latest/ is good imo. > > You can look at packages that are using Documenter here: > http://michaelhatherly.github.io/Documenter.jl/latest/man/examples/ > > On Wednesday, May 11, 2016 at 4:04:58 PM UTC+2, Chris Rackauckas wrote: >> >> Hi, >> I was wonder

[julia-users] Re: Package Documentation Tools

2016-05-11 Thread Chris Rackauckas
tion (see > http://docs.julialang.org/en/release-0.4/manual/documentation/). As far > as I know they have to be written in Markdown. I do not know how well > Sphinx and Markdown work together (Sphinx usually uses files in > reStructuredText (rst) format). > > On Wednesday, May

[julia-users] Re: Package Documentation Tools

2016-05-11 Thread Chris Rackauckas
; > Documenter.jl itself has nothing to do with LaTeX or not. No one prevents > you from using LaTeX in your dostrings or Markdown files. > > > On Wednesday, May 11, 2016 at 5:20:36 PM UTC+2, Chris Rackauckas wrote: >> >> Forgot to mentioned that LaTeX is a requirement. I d

[julia-users] Re: How to vectorize a function which has a tuple as return value to a function which returns a tuple of matrices instead of a matrix with tuple elements?

2016-05-16 Thread Chris Rackauckas
I make this function to solve this problem in one of m codes: function shapeResult(res) #Input is a vector of tuples #Output the "columns" as vectors out = cell(length(res[1])) for i = 1:length(res[1]) out[i] = Vector{typeof(res[1][i])}(length(res)) end for i = 1:length(res), j =

[julia-users] Re: Julia and Slurm 'Hello World'

2016-05-16 Thread Chris Rackauckas
XSEDE's Comet uses slurm. Here's a tutorial that will get you started. While cluster managers will work, the machinefile system is dead simple and I found it integrated the easiest into many different machines. On

Re: [julia-users] LQ decomposition.

2016-05-16 Thread Chris Rackauckas
Lol. On Monday, May 16, 2016 at 12:21:49 AM UTC-7, Sheehan Olver wrote: > > Especially when it's a discussion board that the lecturer posts > regularly... > > On Sunday, May 15, 2016 at 10:56:14 PM UTC+10, Stefan Karpinski wrote: >> >> This appears to be a homework question. Just as on

[julia-users] Re: Julia Utopia: Share your tips and tricks to efficient coding in Julia

2016-05-15 Thread Chris Rackauckas
I work on numerical/scientific code, so my experience may be different than more traditional programming uses. Atom is great. I tried the original JunoLT, and it put me off from Julia for awhile, but I feel at home in Atom. Nice tip: there's a package for hidpi if you have a 4K screen. It

[julia-users] Re: PLots kills the kernel ?

2016-05-18 Thread Chris Rackauckas
With what version of matplotlib? I used the normal install of matplotlib on CentOS 7 but it was matplotlib 1.2, it needs matplotlib 1.5 or it will have this error. On Wednesday, May 18, 2016 at 6:13:04 AM UTC-7, Henri Girard wrote: > > using Plots > pyplot(reuse=true) > x=y=linspace(-5,5,30) >

Re: [julia-users] Too many processes for non-parallel Julia computation

2016-05-06 Thread Chris Rackauckas
No, -p 0 makes there be no worker processes which is the default, but that won't change the linear algebra threads. Instead, if you want to have BLAS use 1 thread, use blas_set_num_threads(1). On Friday, May 6, 2016 at 6:02:01 AM UTC-7, Jorge Fernández de Cossío Díaz wrote: > > I am doing

[julia-users] Re: SIMD multicore

2016-04-16 Thread Chris Rackauckas
BLAS functions are painstakingly developed to be beautiful bastions of parallelism (because of how ubiquitous their use is). The closest I think you can get is ParallelAccelerator.jl's @acc which does a lot of optimizations all together. However, it still won't match BLAS in terms of its

[julia-users] New Package Looking for Contributors: VectorizedRoutines.jl

2016-07-14 Thread Chris Rackauckas
Hey, After some discussion (and letting the idea hang around for a long time) I decided to create VectorizedRoutines.jl. The idea behind this package is to include the useful and familiar vectorized routines that one knows and loves from MATLAB/R/Python. An example of such functions include

[julia-users] Re: What packages, features, other strengths would you recommend when showing Julia to molecular neurobiologists?

2016-07-25 Thread Chris Rackauckas
It seems like most of what they do is biostatistics/bioinformatics. I would show them PyCall and RCall. Knowing that you easily have all of those libraries (and your previous libraries) is great. Also show them the JuliaStats stuff. In fact, ask them what they'd want to add to Julia if they

[julia-users] Re: GR subplots

2016-07-23 Thread Chris Rackauckas
One way to easily do this would be to use Plots.jl with GR as the backend. This way you get the API of Plots.jl and the speed of GR. This may not be the answer you were looking for, but I'd give it a try. Plots.jl is really well documented and has plenty of examples. On Saturday, July 23, 2016

[julia-users] Re: Julia at the SIAM Annual Conference

2016-07-21 Thread Chris Rackauckas
Wait until next year. I can't present on Julia codes there until they're published (adviser's rules). I think in general that would apply a to a lot of people: started getting involved during v0.3, did a project in v0.4 which is now submitted, but won't be out there to present until next year.

[julia-users] ANN: DifferentialEquations.jl

2016-08-01 Thread Chris Rackauckas
I am pleased to announce the first release of DifferentialEquations.jl . An accompanying blog post explains the motivation and philosophy of the package in more detail.

[julia-users] Issues with defining get: how exactly is arr[end] defined?

2016-07-30 Thread Chris Rackauckas
I thought arr[end] just called endof(arr), but that seems to not be the case. For example: type TestType data::AbstractArray end Base.length(t::TestType) = length(t.data) Base.size(t::TestType) = (length(t.data)) Base.endof(t::TestType) = t.data[2,:,:] Base.getindex(t::TestType,i::Int) =

[julia-users] Re: [ANN] Nemo 0.5 released

2016-07-31 Thread Chris Rackauckas
Hey, I am really interested in trying Nemo since I want to use the ArbFloats and ArbReals in DifferentialEquations.jl for "faster BigFloats". However, I am getting the following error when trying to build Nemo: julia> Pkg.build("Nemo") INFO: Building Nemo WARNING: `@windows` is deprecated,

[julia-users] Re: Problems with A\b and BigFloat

2016-08-10 Thread Chris Rackauckas
Though I don't know if they have sparse algorithms. But they have a good base something there to help you get started making one... On Wednesday, August 10, 2016 at 2:20:54 PM UTC-7, Chris Rackauckas wrote: > > GenericSVD.jl <https://github.com/simonbyrne/GenericSVD.jl> has line

[julia-users] Re: Making a New Repl Mode

2016-08-14 Thread Chris Rackauckas
Another one to look at may be SJulia. On Saturday, August 13, 2016 at 11:24:18 PM UTC-7, Nathan Smith wrote: > > Hey all, > > I'm looking at making a new repl mode. I can't find much documentation on > this process. Does anyone have any resources or suggestions to get me > started before I dig

Re: [julia-users] Help please: how to expand dictionary into (local ?) variables with the same name as the keys and the corresponding values

2016-08-14 Thread Chris Rackauckas
Any issues with @materialize? It seems like it does exactly what you're asking for... On Sunday, August 14, 2016 at 3:13:37 AM UTC-7, Adrian Salceanu wrote: > > Variables contained in a module and then parsed Julia code included within > a function using include_string(). > > Any obvious

[julia-users] Re: Help please: how to expand dictionary into (local ?) variables with the same name as the keys and the corresponding values

2016-08-13 Thread Chris Rackauckas
There's a PR for this on Parameters.jl. See this: https://github.com/mauro3/Parameters.jl/pull/13 On Saturday, August 13, 2016 at 1:04:19 AM UTC-7, Adrian Salceanu wrote: > > Hi, > > This seems to be a recurring question per my googlings, but still I could > not find a satisfactory answer.

[julia-users] Re: Problems with A\b and BigFloat

2016-08-10 Thread Chris Rackauckas
Yes textbook answer is, why do you want to use `\`? Iterative techniques are likely better suited for the problem. There's no need to roll you own, the package IterativeSolvers.jl has a good number of techniques implemented which are well-suited for the problem since A is a large sparse matrix.

[julia-users] Re: Problems with A\b and BigFloat

2016-08-10 Thread Chris Rackauckas
e usually need a preconditioner to get any convergence with the > non-symmetric ones for interesting ranks. (I've been struggling for months > to find a good preconditioner for an application of GMRES in my work, so > this is a sore point.) > > On Wednesday, August 10, 2016 at

[julia-users] Re: Problems with A\b and BigFloat

2016-08-10 Thread Chris Rackauckas
GenericSVD.jl has linear solver routines which work for generic number types (like BigFloat). You can use an SVD to solve the linear system. It's not as fast as other methods, but you may find this useful. On Wednesday, August 10, 2016 at 12:47:10

[julia-users] Package Structure for Conditional Dependencies

2016-07-12 Thread Chris Rackauckas
Hey, I wanted to look more into using conditional dependencies, and am looking for insight from both developers who have done this successfully and feedback from users on what is more intuitive. Just so we're all on the same page, a conditional dependency is a package which your package

[julia-users] Re: Coveralls and coverage issues

2016-07-21 Thread Chris Rackauckas
updating? > > Best > Simon > > On Thursday, July 21, 2016 at 9:06:51 AM UTC+1, Chris Rackauckas wrote: >> >> Look at the files it's trying to cover... it's DataFrames.jl :) >> >> I sent you a pull request to fix your travis.yml to be for your package. >>

[julia-users] Re: Help Julia win a performance comparison!

2016-07-21 Thread Chris Rackauckas
if this allocation can be dealt with better (the R code has both of these same issues, which is why we're virtually tied). On Thursday, July 21, 2016 at 5:49:23 AM UTC-7, Steven G. Johnson wrote: > > > > On Thursday, July 21, 2016 at 5:37:12 AM UTC-4, Chris Rackauckas wrote: >> >>

[julia-users] Help Julia win a performance comparison!

2016-07-20 Thread Chris Rackauckas
Hey, there's a performance comparison on Code Review . I gave it a go and with help from @TotalVerb via Gitter, I got Julia really close to being

[julia-users] Re: Help Julia win a performance comparison!

2016-07-20 Thread Chris Rackauckas
Most of the arrays are changing size each time though, since they represent a population which changes each timestep. On Wednesday, July 20, 2016 at 6:47:39 PM UTC-7, Steven G. Johnson wrote: > > It looks like you are allocating lots of arrays in your doStep inner-loop > function, so I'm sure

  1   2   3   >