Re: [julia-users] Is FMA/Muladd Working Here?

25910 Should I open an issue? Note that this is on v0.6 Windows. On Linux the sysimg isn't rebuilding for some reason, so I may need to just build from source. On Wednesday, September 21, 2016 at 6:22:06 AM UTC-7, Erik Schnetter wrote: > > On Wed, Sep 21, 2016 at 1:56 AM, Chris Rac

Re: [julia-users] Is FMA/Muladd Working Here?

m0, %xmm0 movabsq $568231040, %rax# imm = 0x21DE8480 vaddsd (%rax), %xmm0, %xmm0 popq%rbp retq nopw%cs:(%rax,%rax) On Wednesday, September 21, 2016 at 6:29:44 PM UTC-7, Erik Schnetter wrote: > > On Wed, Sep 21, 2016 at 9:22 PM, C

Re: [julia-users] Is FMA/Muladd Working Here?

t; function attribute unsafe-fp-math=true. > > > > > -erik > > > > On Wed, Sep 21, 2016 at 9:33 PM, Yichao Yu > wrote: > >> > >> On Wed, Sep 21, 2016 at 9:29 PM, Erik Schnetter > > >> wrote: > >> > On Wed, Sep 21, 201

[julia-users] Re: Is there a way to download a copy of Plots' documentation?

You have to do what Sundara describes. This is a limitation of Documenter.jl with the mkdocs render. That would be a feature request for Documenter.jl . I think it's planned for the new native renderer, but I couldn't find the issue (and Plots would h

[julia-users] Re: Benchmarking Julia

Is your code still type-stable and the type are inferred correctly? There are some type-inference issues that came up when I transitioned a lot of things to v0.5. After I got rid of those, the performance on v0.5 ended up being a bit better (most likely due to -O3 optimization). I transitioned

[julia-users] Re: PkgDev.tag issues

It happened to me as well when I tagged yesterday. I just changed computers over to Linux... v0.5 Windows 10. On Monday, September 26, 2016 at 8:26:00 PM UTC-7, Tony Kelman wrote: > > The "no changes to commit" issue sounds like > https://github.com/JuliaLang/PkgDev.jl/issues/28 >

[julia-users] Re: Why does Julia 0.5 keep complaining about method re-definitions?

Or just a standard way to suppress warnings of a given type (say, surpress("MethodDefinition")). For now, Suppressor.jl does well. On Tuesday, September 27, 2016 at 12:13:00 PM UTC-7, Andrew wrote: > > It seems like a lot of people are complaining abo

Re: [julia-users] Re: Why does Julia 0.5 keep complaining about method re-definitions?

You could've just used Suppressor.jl ... On Tuesday, September 27, 2016 at 9:55:53 PM UTC-7, K leo wrote: > > > On Wednesday, September 28, 2016 at 12:53:12 PM UTC+8, K leo wrote: >> >> This a very heavy install. It's fetching tons of things that I hav

Re: [julia-users] Re: [ANN] UnicodeFun

what about the other direction? It would be a sick to write symbolic code for SymEngine or SymPy in that unicode form, and have it convert to the appropriate code. On Wednesday, September 28, 2016 at 6:43:48 AM UTC-7, Simon Danisch wrote: > > That's the short form which works with sub/superscrip

[julia-users] Southern California Julia Users Kickoff Meeting - October 7th

Hello, We are excited to announce a kickoff meeting for the Southern California Julia Users on October 7th at 6PM at UCLA (Engineering 4). There will be presentations by local Julia users showcasing their work. If you are interested in presenting, please let us know through email, Gitter, or t

[julia-users] Re: ANN: CUDAdrv.jl, and CUDA.jl deprecation

Thanks for the update. On Thursday, September 29, 2016 at 6:31:29 PM UTC-7, Tim Besard wrote: > > Hi all, > > CUDAdrv.jl is Julia wrapper for > the CUDA driver API -- not to be confused with its counterpart CUDArt.jl >

[julia-users] ANN: ParallelDataTransfer.jl

ParallelDataTransfer.jl is a library for sending and receiving data among processes defined using Julia's parallel computing framework. This library can be used to: - Send variables to worker processes - Directly define variabl

[julia-users] Changing Repositories Without Deleting METADATA

Does anyone have a good way to change repositories? A common example for me is, Sundials is in JuliaDiffEq, so I fork it to my Github account for an extended PR, but to work on it I need to remove my current Sundials install and clone from my own repository. However, METADATA does not like this

[julia-users] Re: Changing Repositories Without Deleting METADATA

is cool-feature > > I'll then create a pull request from Github's ui. > > On Sunday, 2 October 2016 18:01:28 UTC+1, Chris Rackauckas wrote: >> >> Does anyone have a good way to change repositories? A common example for >> me is, Sundials is in JuliaDiffEq, so I

[julia-users] Re: View (a portion of a) vector as a matrix

Fengyang's reshape((@view x[1:6]), (3, 2)) will work well and will be essentially cost-free since reshape creates a view, and a view of a view is still just a view (no copy). Another way to write it is reshape(view(x,1:6), (3, 2)). For example: function f(t,u,du) x = reshape(view(x,1:6), (3,

[julia-users] Re: Changing Repositories Without Deleting METADATA

Thanks for the suggestion. For some reason I never checked to see if there were other GUIs, but once you mentioned it I did quite a bit of Googling. I am now using GitKraken and while it will take a bit to get used to, it's already improving my productivity. Thanks for the suggestion. On Sunday

[julia-users] Re: multiple processes question

See this blog post . If your code is perfectly efficient, yes then processes equal to the number of cores (so for something like BLAS where it's written as the most efficient threaded algorithms you could image). But for your simple homework assignment

[julia-users] Re: Julia vs Seed7, why languages succeed or fail

Who's the audience for Seed7? I googled Seed7 BLAS, Seed7 Linpack, Seed7 FFT and nothing came up. So a large portion of Julia users are not the Seed7 audience. To me, there is almost no similarity between Julia and Seed7, even if the syntax or features were similar. But for the reasons you say,

[julia-users] New Julia-Focused Blog: Pkg Update

Hey, I am happy to announce Pkg Update, a new blog focused on summarizing the large changes throughout the Julia package ecosystem. As I feel it's always easier to get people involved once you already have things up and running, I went ahead and created the first blog post

[julia-users] Re: Julia and the Tower of Babel

Create a repo where we can all bikeshed different names, agree upon some, and then standardize. I honestly don't care which conventions are chosen and will just find/replace with whatever people want, but there has to be a "whatever people want" to do that. On Saturday, October 8, 2016 at 1:47:

[julia-users] Re: Julia and the Tower of Babel

Conventions would have to be arrived at before this is possible. On Saturday, October 8, 2016 at 3:39:55 AM UTC-7, Traktor Toni wrote: > > In my opinion the solutions to this are very clear, or would be: > > 1. make a mandatory linter for all julia code > 2. julia IDEs should offer good intellisen

[julia-users] Re: New SPEC - open to Julia[applications]?

>From your second link: >- Submissions for the first step in the search program will be >accepted by SPEC beginning 11 November 2008 and ending 30 June 2010 (11:59 >pm, Pacific Standard Time). > > On Saturday, October 8, 2016 at 12:18:53 PM UTC-7, Páll Haraldsson wrote: > > > https:

[julia-users] Re: macro: with

What about Parameters.jl ? On Wednesday, September 7, 2016 at 7:23:47 AM UTC-7, Tom Breloff wrote: > > Hey all... I just threw together a quick macro to save some typing when > working with the fields of an object. Disclaimer: this should not be used > i

[julia-users] Re: Setting/getting field-like members with side effects

Why not use Symbols instead of strings here? On Sunday, October 9, 2016 at 8:26:57 AM UTC-7, Bart Janssens wrote: > > Hi all, > > I'm thinking about how to translate a Python interface that makes > extensive use of __getattr__ and __setattr__ overloading to allow chaining > a series of option va

Re: [julia-users] Re: macro: with

anything like this in Parameters. > > On Sun, Oct 9, 2016 at 1:27 PM, Chris Rackauckas > wrote: > >> What about Parameters.jl <https://github.com/mauro3/Parameters.jl>? >> >> >> On Wednesday, September 7, 2016 at 7:23:47 AM UTC-7, Tom Breloff wrote: >>>

[julia-users] Re: Setting/getting field-like members with side effects

Missed that one, though it should be mysolver[:linear_system][:parameters][:solver_type] = :GMRES There's no reason for an algorithm choice to be a String. But this entire thing seems wrong-headed. The choice of the solver method should likely be done by dispatching on the solve method somehow.

Re: [julia-users] Re: Setting/getting field-like members with side effects

a lot > less headaches for students I hope), but rewriting the entire framework is > just too much work right now, so I think we are stuck with this overall > structure for the time being. > > On Sun, Oct 9, 2016 at 9:13 PM Chris Rackauckas > wrote: > >> Mis

[julia-users] Re: Root finding package

I don't know if NLsolve handles complex roots but I've always found it to be very good. Maybe you can just act like the problem is on a vector of two points (the real and imaginary parts) and solve for where the norm of f(x) is zero. On Saturday, October 15, 2016 at 4:56:23 PM UTC-7, digxx wrot

[julia-users] Re: matrix multiplication in Matlab is much faster

Take a look at the performance tips . The first time you run it, the function will compile. Then the compiled function is cached. On my computer I did: a = rand(1000,1000) y=similar(a) @time a*a @time a*a @time A_mul_B!(y,a,a)

[julia-users] Re: Status of Plots.jl?

using Plots #Pkg.add("GR") gr() # Change the backend plot(rand(4,4)) There's a bug with the plot pane where you might need to hit it twice. If that's not working, then it's not setup correctly. On Sunday, October 16, 2016 at 9:45:07 AM UTC-7, missp...@gmail.com wrote: > > Hi folks, > > I don't

[julia-users] Re: ANN: ParallelAccelerator v0.2 for Julia 0.5 released.

Thank you for all of your amazing work. I will be giving v0.2 a try soon. But I have two questions: 1) How do you see ParallelAccelerator integrating with packages? I asked this in the chatroom, but I think having it here might be helpful for others to chime in. If I want to use ParallelAcceler

[julia-users] Re: ANN: ParallelAccelerator v0.2 for Julia 0.5 released.

e former > has the potential to be a lot easier to use particularly for scientific > programmers who are more scientist than programmer. In general, I feel > there is room for all approaches to be supported across a range of > programming ability. > > On Thursday,

[julia-users] Re: so many plotting packages

Just use Plots.jl. JuliaPlots and Plots.jl is essentially a metapackage/org which puts this all together into one convenient package. It works very well and should be recommended as the standard plotting package to almost everyone. On Friday, October 28, 2016 at 12:17:56 PM UTC-7, Ben Arthur wr

[julia-users] Re: Defining a type that inherits methods

Functions shouldn't be written for concrete types and instead for abstract types. If you write your function for AbstractMatrix and then make your type <: AbstractMatrix, this will work naturally. Making the type declarations on a function too strict doesn't help performance anyways. On Saturda

[julia-users] Re: 3D plot (Plots)

Plotly/PlotlyJS can also do some 3d. So can GR. Just look at the chart of supported attributes: https://juliaplots.github.io/supported/ On Saturday, October 29, 2016 at 5:29:06 PM UTC-7, Sheehan Olver wrote: > > surface(x,y,z) > > (or maybe surface(x,y,z')). > > Note that if x and y are matrices

[julia-users] Re: Cost of @view and reshape

reshape makes a view, and views are cheap. Don't worry about this. BTW, I would love to add a collocation method to JuliaDiffEq. Would you consider making this a package? On Sunday, October 30, 2016 at 3:52:37 AM UTC-7, Alexey Cherkaev wrote: > > I'm writing RadauIIA (for now, fixed order 5 with

[julia-users] Re: What's julia's answer to tapply or accumarray?

For reference I've been gathering these kinds of "vectorized" functions in, well, VectorizedRoutines.jl . I am just trying to get an implementation of all of those vectorized routines you know and love since, in some cases, they lead to

[julia-users] Re: Changing label fontsize in Plots.jl

Wonderful to see the learning process in action haha. For future reference, to see which commands are available in which packages, you can check the supported attributes page of the documentation. On Monday, October 31, 2016 at 3:05:30 PM UTC-7, Nitin A

Re: [julia-users] inconsistent 'unique' in Atom

Just click on the number and it will expand it. On Sunday, October 30, 2016 at 7:28:47 PM UTC-7, missp...@gmail.com wrote: > > Hi Yichao, > > thanks a lot, > it does display it correctly if I use dump, but it's annoying that Atom is > inconsistent while displaying the results > > thanks a lot, >

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

This is pretty much obsolete by the . fusing changes: A .= A.*B should be an in-place update of A scaled by B (Tomas' solution). On Tuesday, November 1, 2016 at 4:39:15 PM UTC-7, Sheehan Olver wrote: > > Should this be added to a package? I imagine if the arrays are on the GPU > (AFArrays) the

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

ve used it... > > On Tue, Nov 1, 2016 at 10:06 PM, Sheehan Olver > wrote: > >> Ah, good point. Though I guess that won't work til 0.6 since .* won't >> auto-fuse yet? >> >> Sent from my iPhone >> >> On 2 Nov. 2016, at 12:55, Chris

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

do this, but I don't think they can fuse. On Tuesday, November 1, 2016 at 8:06:12 PM UTC-7, Sheehan Olver wrote: > > Ah thanks! > > Though I guess if I want the same code to work also on a GPU array then > this won't help? > > Sent from my iPhone > > On 2 Nov.

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

Nice work! On Wednesday, November 2, 2016 at 1:14:06 PM UTC-7, Michael Hatherly wrote: > > I’m pleased to announce the initial 0.1 release of Highlights.jl > — a Julia package for > highlighting source code similar to the well-known Python package cal

[julia-users] Re: PSA: Documenter.jl deprecations

Once again, thank you very much for this change! It's highly appreciated by those of us who are unfamiliar with gem, Ruby, travis, and the whole install feature. As you know, I gave this new access token setup a go and had it working in minutes with no extra software to install. Just generate s

[julia-users] Re: Sharing experience on packages

For now, stars are the best bad measurement we have. On Tuesday, November 15, 2016 at 8:05:16 AM UTC-8, Jérôme Collet wrote: > > Hi all, > > > I am new to Julia, I used to use R. And using R packages, the main > difficulty for me is the choice of a package for a given task. Most of the > time, t

[julia-users] Re: Use of special characters in labels of plots.jl

Use LaTeXStrings.jl's L_str macro: L"$10^{10}" I know this works with the PyPlot backend, and I think GR. On Wednesday, November 16, 2016 at 4:37:16 AM UTC-8, Ferran Mazzanti wrote: > > ...and please notice that I'm not looking for explicit LaTeX suppport. > Just wanted to know if there is an

Re: [julia-users] How to tell if Val{T} is using a fast route

The top level scope of a module is still a global (interactive scope). On Wednesday, November 16, 2016 at 7:59:04 AM UTC-8, FANG Colin wrote: > > Typo, should be > > module ... > > ff(x::Type{Val{1}}) = 1 > > x = 1 > a = ff(Val{x}) > > const y = 1 > a = ff(Val{y}) > > end > >> >>

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

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] Re: Julia Utopia: Share your tips and tricks to efficient coding in Julia

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 will

[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?

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 = 1

Re: [julia-users] LQ decomposition.

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 StackOverf

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

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 M

[julia-users] Re: Package Documentation Tools

g for Fedora version, and the download > link seems broke. Others seeing this as well (at > https://copr.fedoraproject.org/coprs/nalimilan/julia/) > > On Wednesday, May 11, 2016 at 10:04:58 AM UTC-4, Chris Rackauckas wrote: >> >> Hi, >> I was wondering if there&#x

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

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) > t

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

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

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

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 every

[julia-users] Re: Parallel computing

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 core

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

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 with

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

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?

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 it

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

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 spots

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

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: I feel that on the syntax level, Julia sacrificed too much elegancy trying to be compatible with textbook math notations

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 o

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

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 reaso

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

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) = [cos(2*π.*x(:,

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

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 for

[julia-users] Julia Memory Management

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 cont

Re: [julia-users] Julia Memory Management

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 > wrote: > > I see mentions like this one every once in

Re: [julia-users] Julia Memory Management

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 > wrote: > > Thanks, that's good to know. Can you explain a little bit about why it > woul

[julia-users] Private Forks of Julia Repositories

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: Private Forks of Julia Repositories

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 branch

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

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 > wrote: > > I think I will nee

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

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 wrote: >> > I think I will need

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

-to-default> to set the default remotes properly, and can merge/rebase between them on my own local repo. After setting this up, it even seems to work with the GUI which is nice. Thanks for the idea! On Friday, June 3, 2016 at 12:55:38 AM UTC-7, Mauro wrote: > > On Fri, 2016-06-03 at

[julia-users] Re: Differential Equations Package

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 also

Re: [julia-users] Differential Equations Package

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 am

Re: [julia-users] Differential Equations Package

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 t

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

That's the main issue here. I do have a way of creating them... in a private branch that I cannot share right now (give it a month?). It's not easy to create them either: essentially you can solve a really stiff SDE and get them, so think about repeatedly taking random numbers in a loop, multip

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

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: >> >&g

Re: [julia-users] Differential Equations Package

gt; isn't that hard to switch around, 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,

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

>> 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 thi

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

Forgot to mention that if you change the deepcopy in the loop to copy, it works.

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

2016 at 6:24 PM, Chris Rackauckas > wrote: > > Sorry, was making benchmarks and accidentally "fixed it". It turns out > to be > > an issue with deepcopying BigFloats. Here's a small test case: > > https://github.com/JuliaLang/julia/issues/16667 > &

Re: [julia-users] Differential Equations Package

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, etc.

[julia-users] Re: Status of FEM packages

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 r

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

Yes, there was an error. Scott Jones helped me fix it. I am working with Randy to make sure the updated version is the one that's re-blogged. Sorry for the error! On Wednesday, June 8, 2016 at 1:14:44 PM UTC-7, Andre Bieler wrote: > > "Switching Float to Int in `test2()` makes the vector of arra

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

It should all be fixed now.

Re: [julia-users] Differential Equations Package

s://github.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. Cheer

Re: [julia-users] Differential Equations Package

> > 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 the

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

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] Require a package not in METADATA

Is there a way to add a package to REQUIRES that's not in METADATA?

[julia-users] Conditional Definition of Functions?

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σ = zeros

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

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 > wrote: > > There seems to be an interesting change between v0.4.5 a

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

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 announcem

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

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 getti

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

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 PyPlot

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

I am not going to comment on the implementation issues because it's outside of my expertise and I will have nothing of substance to contribute. But as a user, I do find it troubling that the only reliable high-precision number type is clearly not aiming for performance. It seems that what would

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

unction Base method. > > Maybe I am missing something 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 thi

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

nguages, and more a > feature of it being 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 Sund

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

at I an use Rational inputs for an ode > solver > > ;), and that all my code would be expected to be filled with eltypes and > > inexact errors. We shall see. > > > > On Sunday, June 19, 2016 at 7:32:55 PM UTC-7, Chris Rackauckas wrote: > >> > >> But

  1   2   3   >