Re: [julia-users] Re: quadgk with 2 arguments

2016-10-18 Thread Mosè Giordano
Hi Steven, 2016-10-19 0:36 GMT+02:00 Steven G. Johnson : > For example: > > quadgk(z -> 1/z, 1, 1im, -1, -1im) > > integrates 1/z over a closed counter-clockwise diamond-shaped contour around > the origin in the complex plane, returning 2πi by the residue theorem. Did

Re: [julia-users] Re: julia-i18n: Translators and reviewer needed!

2016-10-07 Thread Mosè Giordano
Hi Ismael, I can confirm the problems in http://julialang.org/downloads/platform.html are now fixed and I added the missing string. Thank you! However I just noticed that the description of Gadfly in http://julialang.org/downloads/plotting.html is back in English (this happens for both Italian

[julia-users] Re: julia-i18n: Translators and reviewer needed!

2016-10-02 Thread Mosè Giordano
Hi, last week I completely translated http://julialang.org/downloads/platform.html to Italian (or at least I thought I did it), but now most of the translations disappeared, without recent changes to https://github.com/JuliaLang/julialang.github.com/blob/master/downloads/platform.md. What

Re: [julia-users] Re: Idea: Julia Standard Libraries and Distributions

2016-09-14 Thread Mosè Giordano
2016-09-15 0:10 GMT+02:00 Chris Rackauckas: > If you use Reexport.jl's @reexport in the top-scope of a module for such a > metapackage it will actually export all of the functions from inside: > > module MyMetapackage # Setup t > using Reexport > @reexport using DifferentialEquations >

Re: [julia-users] Re: Idea: Julia Standard Libraries and Distributions

2016-09-14 Thread Mosè Giordano
2016-09-14 23:46 GMT+02:00 Chris Rackauckas: > I too am weary about how different distributions would play together. I'd > assume you'd just install with one and that would be it. There are many people working in many fields that may want to install different distributions. I don't think that

[julia-users] Re: Idea: Julia Standard Libraries and Distributions

2016-09-14 Thread Mosè Giordano
Hi Chris, what would be the difference between a distribution and a meta-package like those used, for example, in Debian and derivatives: a package without code, that only requires other packages? In this sense you can create right now a meta-package: just create a repository for a Julia

[julia-users] [ANN] Deconvolution.jl: library for deconvolution of digital signals

2016-09-08 Thread Mosè Giordano
Hi all, I'm pleased to announce a new package: Deconvolution.jl . It is part of the JuliaDSP project and aims at providing a set of tools to deconvolve digital signals. The package is released under the terms of MIT

[julia-users] Re: How to shuffle the order of matrix by rows

2016-09-07 Thread Mosè Giordano
Hi Ahmed, Hello, > l have a matrix as follow : > k=rand(1:100,3,3) > 3x3 Array{Int64,2}: > 97 26 77 > 58 64 86 > 87 85 99 > > > l want to shuffle randomly the order by lines for instance line 1 becomes > 3. l have a matrix of 1 by 700 . > The following function is based on the

Re: [julia-users] AutoGrad port for Julia

2016-08-29 Thread Mosè Giordano
Hi Tom, Mose: what version of julia are you on? Anonymous functions and closures > are much faster on 0.5... In fact there should be no performance penalty vs > regular functions, which allows you to rethink your paradigm. > It was Julia 0.4.6, but I get similar results also with Julia

[julia-users] Re: AutoGrad port for Julia

2016-08-28 Thread Mosè Giordano
Hi Deniz, Announcing AutoGrad.jl : an > automatic differentiation package for Julia. It is a Julia port of the > popular Python autograd package. It > can differentiate regular Julia code that includes loops,

[julia-users] Re: [ANN] LombScargle.jl: compute periodogram for unevenly sampled data

2016-08-19 Thread Mosè Giordano
Hi all, an update about the situation of LombScargle.jl <https://github.com/giordano/LombScargle.jl> 2016-07-18 10:15 GMT+02:00 Mosè Giordano: > In the future I may implement another much-faster Lomb-Scargle algorithm > by Press & Rybicki (1989, ApJ, 338, 277), which however r

[julia-users] Re: overwrite function value

2016-08-11 Thread Mosè Giordano
> > so when I define > > function f!(x) > x=x+1 > end > > and then > x=1 > f!(x) > I get 2 returned but I also want x to be replaced by the newly calculated > value :/ > There is an entry of the FAQ about this:

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

2016-08-11 Thread Mosè Giordano
Another relevant update for this package: Cuba.jl can now be used on Windows (GNU/Linux and Mac OS were already supported). I managed to cross compile the shared library for Windows on GNU/Linux, so during building of the package a prebuilt dll is downloaded for Windows users. If you want to

[julia-users] [ANN] LombScargle.jl: compute periodogram for unevenly sampled data

2016-07-18 Thread Mosè Giordano
Dear all, I'm pleased to announce the first release of LombScargle.jl , a package to compute the Lomb-Scargle periodogram . Differently from standard FFT, this can be used to find

[julia-users] Re: A Very Simple Benchmark for Brutal-force Loops in Several Languages: revised, Julia is fast!

2016-07-11 Thread Mosè Giordano
Hi Zhong, you may want to check out Julia 0.5, on my box your benchmark is ~13% faster with Julia 0.5: 3.367658 seconds with Julia 0.4.6 and 2.898068 seconds with Julia 0.5. Bye, Mosè Hi, > > Sorry I have to post a revision so quickly. Just after I posted my > previous benchmark, I found

[julia-users] Re: [ANN] Measurements.jl v0.1.1: Error propagation library

2016-06-27 Thread Mosè Giordano
Hi, 2016-06-27 4:32 GMT+02:00 g: > > Hello, > This looks very nice. I will use it next time I have a need. > Thank you :-) > I took a brief look at the code, and I have one question. Why did you > choose to "tag" your Measurement types with a random Float64, instead of an >

[julia-users] [ANN] Measurements.jl v0.1.1: Error propagation library

2016-06-26 Thread Mosè Giordano
Hi all, I'm really happy to announce availability of Measurements.jl , v0.1.1, a full-fledged uncertainty propagation library. This is the list of features: - Support for most

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

2016-06-23 Thread Mosè Giordano
Hi Mauro, 2016-06-22 21:13 GMT+02:00 Mauro: > > A REPL session of mine trying to figure out how to make that macro would > look something like this: > > julia> ex = :(f(x, y, z, zz)) # this is what is passed into the macro > :(f(x,y,z,zz)) > > julia> xdump(ex) # xdump is nice to get an

[julia-users] Re: Help with a macro

2016-06-22 Thread Mosè Giordano
Ok, here is a simpler example, not specific to Measuments.jl. From @foo f(x, y, z, ...) get f([x.re, y.re, z.re, ...].^2 + [x.im, y.im, z.im, ...]) The expression provided to the macro should can take any number of arguments. For example: x = complex(3.2, -5.2) y = complex(8.9, 0.4)

[julia-users] Re: Help with a macro

2016-06-21 Thread Mosè Giordano
2016-06-21 23:02 GMT+02:00 Kristoffer Carlsson: > > > Your macro should not evaluate anything per say. It should generate > expressions that when the code is run will call functions. > > Say a user types @uncertain log(x,y), what do you want the resulting > expression after macro expansion has

[julia-users] Re: Help with a macro

2016-06-21 Thread Mosè Giordano
Hi Kristoffer, 2016-06-21 22:19 GMT+02:00 Kristoffer Carlsson: > > Using @eval in a macro is generally the wrong thing to do. A macro should > simply transform an expression to another expression. Ok, that's probably what I'm missing, thank you! > If you give an example function body and

[julia-users] Help with a macro

2016-06-21 Thread Mosè Giordano
Hi all, in Measurements.jl , an uncertainty propagation library, I added a macro to propagate uncertainty for an arbitrary function. Here you can see some examples of use of the macro:

Re: [julia-users] Re: [ANN & RFC] Measurements.jl: Uncertainty propagation library

2016-06-17 Thread Mosè Giordano
2016-06-16 10:12 GMT+02:00 Mosè Giordano: > Hi all, > > an update: I continued the development of this package, you can have a look > at the available features here: > https://github.com/giordano/Measurements.jl#features-list > > In addition, I implemented support for correlat

[julia-users] Re: [ANN & RFC] Measurements.jl: Uncertainty propagation library

2016-06-16 Thread Mosè Giordano
Hi all, an update: I continued the development of this package, you can have a look at the available features here: https://github.com/giordano/Measurements.jl#features-list In addition, I implemented support for correlation between variables

[julia-users] Regression of promote with 6 or more arguments in Julia 0.5

2016-06-01 Thread Mosè Giordano
Hi all, I'm noticing a large (factor of >100) regression of performance of promote function when fed with more than 5 arguments ("julia" is Julia 0.4.5; "./julia" is the official current nightly build, version 0.5.0-dev+4438, commit aa1ce87): % julia -e 'using Benchmarks;print(@benchmark

Re: [julia-users] Re: Is Julia slow with large arrays?

2016-06-01 Thread Mosè Giordano
Oh my bad, this was really easy to fix! I usually do inspect the code with @code_warntype but missed to do the same this time. Lesson learned. Now the same loop is about ~30 times faster in Julia than in Fortran, really impressive. Thank you all for the prompt comments! Bye, Mosè 2016-06-01

[julia-users] Is Julia slow with large arrays?

2016-06-01 Thread Mosè Giordano
Hi all, I'm working on a Fortran77 code, but I'd much prefer to translate it to Julia. However, I've the impression that when working with large arrays (of the order of tens of thousands elements) Julia is way slower than the equivalent Fortran code. See the following examples: Julia:

[julia-users] Re: Quadratic equation solver ?

2016-05-31 Thread Mosè Giordano
> > which one is a good solver(open source) for quadratic equation using jump > ? Thank you > What is jump? PolynomialRoots.jl is a package to find all roots (real and complex, they will have always Complex type) of polynomials of any degree,

Re: [julia-users] Re: [ANN] PolynomialRoots.jl: Fast Complex Polynomial Root Finder

2016-05-27 Thread Mosè Giordano
gt; On Thu, May 5, 2016 at 7:39 PM, Mosè Giordano <mose.g...@gmail.com > > wrote: > >> Hi Steven and Stefan, >> >> 2016-05-05 22:29 GMT+02:00 Stefan Karpinski <ste...@karpinski.org >> >: >> > The phrasing of this licensing is unclear

Re: [julia-users] Passing a function outputted by a higher order function to C - function not yet c-callable

2016-05-22 Thread Mosè Giordano
Hi Michael, > Thank you for the quick reply. Could you please elaborate on how to pass > the data explicitly as you suggest? I don't quite understand. > This is a very useful post about closures in general: http://julialang.org/blog/2013/05/callback Bye, Mosè

Re: [julia-users] Re: [ANN & RFC] Measurements.jl: Uncertainty propagation library

2016-05-19 Thread Mosè Giordano
Hi Andre, 2016-05-19 9:05 GMT+02:00 Andre Bieler: > As an experimental physicist I thank you very much for this package! :D You're welcome :-) If you have wishes or suggestions, please share them! Bye, Mosè

[julia-users] Re: [ANN & RFC] Measurements.jl: Uncertainty propagation library

2016-05-17 Thread Mosè Giordano
Hi Lucas, Looks great! > > How about overloading the ± operator? > using Measurements > > ±(μ, σ) = Measurement(μ, σ) > # => ± (generic function with 1 method) > > a = 4.5 ± 0.1 > # 4.5 ± 0.1 > > > I don't know if it's defined in other packages, but I've never seen it. > This is a great idea,

[julia-users] [ANN & RFC] Measurements.jl: Uncertainty propagation library

2016-05-17 Thread Mosè Giordano
Dear all, I'd like to share with you my new simple package Measurements.jl , a library for propagation of uncertainty. The package provides a new type, Measurement, that is composed by the value and its associated error, assumed to be a standard

Re: [julia-users] Re: [ANN] PolynomialRoots.jl: Fast Complex Polynomial Root Finder

2016-05-05 Thread Mosè Giordano
Hi Steven and Stefan, 2016-05-05 22:29 GMT+02:00 Stefan Karpinski : > The phrasing of this licensing is unclear but essentially the same as the > original Fortran library: > >> The authors release the source codes associated with the Paper under terms >> of the GNU Lesser

[julia-users] [ANN] PolynomialRoots.jl: Fast Complex Polynomial Root Finder

2016-05-05 Thread Mosè Giordano
Hi, a couple of days ago I released version 0.0.2 of PolynomialRoots.jl (formerly known as CmplxRoots.jl, some people suggested to change the name to something more descriptive), a fast root finder of real and complex polynomials. This is a

Re: [julia-users] Re: [ANN] CmplxRoots.jl: Fast Complex Polynomial Root Finder

2016-04-30 Thread Mosè Giordano
2016-04-30 14:22 GMT+02:00 Sheehan Olver: > > Hi Mose, > > I think AMVW is specifically designed for large degree polynomials: it's > an O(n^2) algorithm that works via orthogonal transformations on the > companion matrix > > How do the timings compare on, say, a degree 200 polynomial? > Just

[julia-users] Re: [ANN] CmplxRoots.jl: Fast Complex Polynomial Root Finder

2016-04-30 Thread Mosè Giordano
Hi Sheehan, 2016-04-30 8:12 GMT+02:00 Sheehan Olver: > > How does this compare to AMVW.jl? > > https://github.com/andreasnoack/AMVW.jl > I compared both programs with this script (requires development version of PolynomialRoots): using Benchmarks, PolynomialRoots, AMVW for polynomial in

Re: [julia-users] Re: [ANN] CmplxRoots.jl: Fast Complex Polynomial Root Finder

2016-04-29 Thread Mosè Giordano
Hi David, 2016-04-30 0:17 GMT+02:00 David P. Sanders <dpsand...@gmail.com>: > > > El jueves, 28 de abril de 2016, 16:09:54 (UTC-4), Mosè Giordano escribió: >> >> Hi David, >> >> 2016-04-28 14:10 GMT+02:00 David P. Sanders <dpsa...@gmail.com>: >>

Re: [julia-users] Re: [ANN] CmplxRoots.jl: Fast Complex Polynomial Root Finder

2016-04-29 Thread Mosè Giordano
2016-04-29 9:58 GMT+02:00 2n : > So translating header files for bindings is considered a derived work also? > And also the GNU guys and the porting of unix sofware to linux, that can't > be considered derived work. I guess it comes down to the definition of the > word

Re: [julia-users] Re: [ANN] CmplxRoots.jl: Fast Complex Polynomial Root Finder

2016-04-28 Thread Mosè Giordano
2016-04-28 22:34 GMT+02:00 Steven G. Johnson <stevenj@gmail.com>: > > > On Thursday, April 28, 2016 at 4:09:54 PM UTC-4, Mosè Giordano wrote: >> >> > +1 for ComplexRoots.jl >> >> Ok, then I'll see how to change the name of the package in METADA

Re: [julia-users] Re: [ANN] CmplxRoots.jl: Fast Complex Polynomial Root Finder

2016-04-28 Thread Mosè Giordano
Hi Steven, 2016-04-28 14:58 GMT+02:00 Steven G. Johnson : > On Thursday, April 28, 2016 at 8:46:06 AM UTC-4, Steven G. Johnson wrote: >> >> Out of curiosity, how does the performance and accuracy compare to the >> "roots" function in the Polynomials.jl package, which uses

Re: [julia-users] Re: [ANN] CmplxRoots.jl: Fast Complex Polynomial Root Finder

2016-04-28 Thread Mosè Giordano
Hi David, 2016-04-28 14:10 GMT+02:00 David P. Sanders : > As far as I could see, the original library has an Apache license, so you > should be able to use MIT. > > I believe that you need to include a copy of the original license in your > package? Original code is

Re: [julia-users] Re: [ANN] CmplxRoots.jl: Fast Complex Polynomial Root Finder

2016-04-28 Thread Mosè Giordano
Hi, 2016-04-28 12:57 GMT+02:00 DNF : > Just a question: Is there some particular reason you chose to to skip the > letters 'o' and 'e' in the name? > > It seems like the kind of thing that would make the package hard to find, > and cause confusion and frustration. In fact, while

Re: [julia-users] Re: [ANN] CmplxRoots.jl: Fast Complex Polynomial Root Finder

2016-04-28 Thread Mosè Giordano
2016-04-27 15:10 GMT+02:00 Mosè Giordano <mose.giord...@gmail.com>: > Hi Kristoffer, > > 2016-04-27 14:32 GMT+02:00 Kristoffer Carlsson <kcarlsso...@gmail.com>: >> This looks interesting. > > Thank you! > >> It is in my opinion a bit unfortunat

Re: [julia-users] Re: [ANN] CmplxRoots.jl: Fast Complex Polynomial Root Finder

2016-04-27 Thread Mosè Giordano
Hi Kristoffer, 2016-04-27 14:32 GMT+02:00 Kristoffer Carlsson : > This looks interesting. Thank you! > It is in my opinion a bit unfortunate that you chose > to use the GPL license for this since it means that it will not be usable by > most packages in Julia which are

[julia-users] [ANN] CmplxRoots.jl: Fast Complex Polynomial Root Finder

2016-04-27 Thread Mosè Giordano
Dear all, I am happy to announce the first version of CmplxRoots.jl , a fast root finder of real and complex polynomials. This is a Julia implementation of the algorithm described in - J. Skowron & A. Gould, 2012, "General Complex Polynomial Root

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

2016-04-13 Thread Mosè Giordano
Hi Kaj, thanks for your question, it served me as inspiration to add another example to the manual: https://cubajl.readthedocs.org/en/latest/#passing-data-to-the-integrand-function Bye, Mosè 2016-04-13 5:14 GMT+02:00 Steven G. Johnson : > > > On Tuesday, April 12, 2016

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

2016-04-13 Thread Mosè Giordano
Hi Oliver, 2016-04-12 21:38 GMT+02:00 Oliver Schulz : > Very nice, Mose! I thought I might have to write a Julia wrapper for Cuba at > some point, but you beat me to it (which I'm grateful for!). :-) Thanks! > I noticed that you're maintaining a modified version of

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

2016-04-13 Thread Mosè Giordano
com>: > 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, >

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

2016-04-12 Thread Mosè Giordano
A quick update: with new version 0.0.5 the syntax of the integrand function has been changed and now the integrand function should be passed as in Cubature,jl, for the case of vector-valued integrands. So, for example, to integrate log(x)/sqrt(x) for x in [0, 1], one can use one of the

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

2016-04-10 Thread Mosè Giordano
Dear all, I am proud to announce Cuba.jl a library for multidimensional numerical integration with four independent algorithms: Vegas, Suave, Divonne, and Cuhre (this algorithm is the same used in Cubature.jl). This package is just a wrapper around Cuba