Re: [julia-users] Julia's current state of web-programming

2016-09-26 Thread Alexey Cherkaev
andle requests, which is pretty > cool. > > If you check here: http://genieframework.com/packages (built with Genie > btw) you’ll find a few more options - like Merly and Bukdu. > > Other than that you have the low(er) level components: Mux, HttpServer, > WebSockets. &g

[julia-users] Julia's current state of web-programming

2016-09-22 Thread Alexey Cherkaev
Hi all, What is the current state of web programming with Julia? http://juliawebstack.org/ seems quite out of date (still suggesting to use Morsel, which is marked as deprecated). Escher looks quite nice (https://github.com/shashi/Escher.jl), but still fails to build on both 0.4 and 0.5

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

2016-10-03 Thread Alexey Cherkaev
Use x and u[7], u[8], u[9], and u[10] to write directly to du > nothing > end > > should be a good way to write the function for Sundials.jl, > DifferentialEquations.jl, ODE.jl (after the iterator PR). > > On Sunday, October 2, 2016 at 5:43:01 AM UTC-7, Alexey Cherkaev wrote:

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

2016-10-02 Thread Alexey Cherkaev
I have the model where it is convenient to represent one of the variables as a matrix. This variable, however, is obtained as a solution of ODEs (among other variables). I'm using Sundials to solve ODEs and `Sundials.cvode` requires the ODE RHS function to take a vector. So, it seems logical

[julia-users] Cost of @view and reshape

2016-10-30 Thread Alexey Cherkaev
I'm writing RadauIIA (for now, fixed order 5 with 3 points) method for ODE BVP (basically, collocation method). In the process, I construct an overall non-linear equation that needs to be solved. It takes "mega-vector" x[j] as an argument. However, internally it is more convenient to reshape it

[julia-users] Julia and Fortran: passing floating point numbers

2016-11-09 Thread Alexey Cherkaev
I have encountered a rather strange error as I'm trying to create a wrapper for Fortran's TWPBVP (ODE BVP) solver. Full code is available at https://github.com/mobius-eng/TWPBVP.jl After downloading (to julia's package directory), run Pkg.build("TWPBVP") and Pkg.test("TWPBVP") The issue

Re: [julia-users] Julia and Fortran: passing floating point numbers

2016-11-09 Thread Alexey Cherkaev
at 3:49:55 PM UTC+2, Ángel de Vicente wrote: > > Hi, > > Alexey Cherkaev <alexey@gmail.com > writes: > > To cut long story short, I have a problem of passing floating-point > numbers from > > Julia to Fortran subroutine. Documentation suggests, that each parame

Re: [julia-users] Julia and Fortran: passing floating point numbers

2016-11-09 Thread Alexey Cherkaev
Found the problem: "double precision" on my machine (presumably, any 64-bit machine) with "-fdefault-real-8" flag will be "real*16" (something like non-existent "Float128" in Julia). On Wednesday, November 9, 2016 at 4:05:59 PM UTC+2, Alexey Cherkaev w

[julia-users] Does new "dot" syntax fuse?

2016-10-20 Thread Alexey Cherkaev
Hi all! Consider example: # collect to produce vectors - not strictly necessary though y = collect(linspace(0.0, 1.0, 1000)) z = collect(linspace(0.5,3.0, 1000)) x = zeros(Float64, 1000) If I do @time x .= sin.(y) The timing output is 0.57 seconds (7 allocations: 208 bytes) So, OK,

[julia-users] Re: Does new "dot" syntax fuse?

2016-10-21 Thread Alexey Cherkaev
Ay, how many times I've done it! Still forget about global variables! OK, I see now about `.+` and friends. Thanks! On Friday, October 21, 2016 at 3:59:59 AM UTC+2, Steven G. Johnson wrote: > > Putting things into a function to avoid benchmarking in global scope: > > > foo!(x,y) = x .=

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

2016-11-01 Thread Alexey Cherkaev
eshape 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