Questions tagged [julia]

Julia is a high-level dynamic programming language designed to address the requirements of high-performance numerical and scientific computing while also being effective for general purpose programming.

Filter by
Sorted by
Tagged with
0 votes
1 answer
32 views

streaming / low memory groupby in julia

I am working on a groupby-aggregation function that will work without RAM overflow issues. Essentially, I want it to run as fast as possible, while not necessarily loading the entire data structure ...
  • 123
0 votes
0 answers
38 views

Ewald summation in Julia

I am writing a small MD (molecular dynamics) engine in Julia and I noticed that 67% of the time spent running the following Ewald summation function (to calculate coulomb energy of an infinite crystal ...
4 votes
1 answer
178 views

Minimize memory allocations while calculating the matrix representation of a Hamiltonian

I am writing a script using Julia (calling my functions from the REPL) to compute the matrix representation of certain Hamiltonian. After I run my script I see that the memory allocations count is ...
0 votes
0 answers
33 views

Numerical solution to the Lane-Emden equation

I am new to Julia programming, but for fun and practice, I wrote a program that numerically solves the Lane-Emden equation. I even wrote an interactive version that changes the values of ...
  • 215
2 votes
0 answers
46 views

Non-allocating k-nearest-neighbours between two matrices in Julia

I played around with some optimizations when answering this SO question. Sought is a function which given matrices x of size (m, l) and ...
3 votes
3 answers
148 views

Project Euler #12 in Julia slower than Python?

My code in Julia, almost identical as the Python code (see below), runs in 4.6 s while the Python version runs in 2.4 s. Obviously there is a lot or room for improvement. ...
2 votes
2 answers
152 views

Sort three numbers

So, I wrote this function to sort three numbers in Julia. Any feedback on readability, performance, and Julian-ness would be appreciated. ...
  • 333
1 vote
0 answers
92 views

How to make my 1D EM algorithm written in Julia faster?

I wrote this code that's supposed to implement a 1-dimensional (!) EM algorithm with k-means initialization: ...
  • 455
2 votes
1 answer
109 views

Reducing allocations in arbitrary precision calculations with Julia

I am new to Julia and have only written in python before which is most likely also reflected by my coding style. Unfortunately, my applications require very high numerical precision which is why I am ...
  • 143
5 votes
2 answers
206 views

Fast hamming distance function in Julia that returns a distance matrix

I would like to use Julia to compute the hamming distance on a very large dataset. I need to get back a distance matrix between rows in order to run further ...
  • 247
3 votes
1 answer
161 views

Discrete event simulation framework for discrete times in Julia

I would like to receive constructive feedback on my simulation framework (this is for my thesis work). I have no formal background in programming, and I am wondering if there is a more effective and ...
  • 31
3 votes
1 answer
148 views

Parsing values from string into struct using match in Julia

My goal I am parsing from a string which contains token:value pairs into a type. Example: ...
2 votes
1 answer
80 views

Return n largest nonzero entries of an array

For something I'm working on, I need a function that takes in a one-dimensional array vec of integers and returns a boolean array of the same shape indicating where ...
  • 207
6 votes
1 answer
271 views

Speed of Python vs Julia for Game of Life. Why is Julia slow?

To compare the speed of Julia to Python+numba, I implemented Game of Life in both languages. For anyone who is not familiar with ...
  • 225
3 votes
2 answers
218 views

Solving 15Puzzle with Julia

Asking here instead of SO as suggested. I'm trying to use Julia to solve the common tile game 15 Puzzle using Julia using A* algorithm. I am quite new to the language and my style may seem very C like....
2 votes
1 answer
52 views

Restructuring loop logic and its effect on allocations

I came across in a code base # TODO this seems very clumsy to duplicate the loop code like this? so I figured I'd give it a shot as a learning exercise. I'd like ...
  • 221
6 votes
0 answers
135 views

Julia set in Julia (and other fractals)

I've made the following to draw a series of zooming in images of a fractal. ...
  • 3,512
2 votes
0 answers
75 views

Sampling and Reconstructing a signal in Julia

I am a newcomer to the Julia world, so I started trying to implement a very simple program: Sample a signal Reconstruct it (various ways, just rectangular for now) Plot it Everything is fine and '...
  • 121
6 votes
0 answers
128 views

How to speed up simple linear algebra optimization probelm in Julia?

I implemented the LSDD changepoint detection method decribed in [1] in Julia, to see if I could make it faster than the existing python implementation [2], which is based on a grid search that looks ...
  • 161
6 votes
1 answer
710 views

Optimizing my 2D Ising model code in Julia

I'm just starting to learn Julia, I work primarily in physics and am used to writing most of my code in Fortran90 and occasionally Python for Tensorflow (also Mathematica but that's less relevant). ...
  • 161
2 votes
1 answer
48 views

Add optional distance range to function

I have a simple brute-force search function: ...
  • 1,453
2 votes
1 answer
173 views

Time-series search with early stopping

I want to search for a pattern in a time-series while either ignoring the mean/shift/bias or the scale/standard deviation. Consequently, I've written two functions. The first function searches ...
  • 1,453
1 vote
1 answer
75 views

Experience replay memory

I have just begun learning Julia. Here is an implementation of experience replay memory for use in a reinforcement learning algorithm. It is pretty simple, essentially a ring buffer with the following ...
  • 111
0 votes
1 answer
82 views

Building an equation [closed]

I've developed a function in the Julia language, using the SymEngine library, which builds an equation. This function works quite quickly, returning a lambdified ...
  • 109
1 vote
0 answers
205 views

WebSocket based API library in Julia

I made a very bare-bones API "library" in Julia which should abstract the interaction with the API. The API is fairly complex - it's based around 2 core layers an HTTPS/REST API for general ...
  • 11
3 votes
1 answer
67 views

Bernstein polynomial type for quadratic clipping

I would like to implement a Bernstein polynomial class in Julia for learning purposes. The goal in the end is to implement the quadratic clipping algorithm from this paper. I am completely new to ...
  • 842
2 votes
2 answers
152 views

For names sake - Project Euler 22

Project Euler - 22 Names scores Using names.txt (right click and 'Save Link/Target As...'), a 46K text file containing over five-thousand first names, begin ...
3 votes
0 answers
256 views

Weighted logistic regression in Julia

I'm trying to estimate a weighted logistic regression as part of a bigger project. I have an implementation in Matlab2015b, but I wanted to give Julia a try to see if I could speed up the estimation ...
  • 131
3 votes
0 answers
678 views

Parsing simple log files in Julia vs Python

I wrote a log parser in Python and Julia. The Python version takes about 8 seconds to process 350 files, while the Julia version takes about 8 seconds to compile, after which it takes over 20 seconds ...
1 vote
0 answers
125 views

Idiomatic string from/to integer conversions in Julia (Project Euler #33)

I am starting to learn Julia. For getting fluent in the syntax, I am solving some of the easier Project Euler challenges. Problem 33 is about so-called curious fractions, which keep their value after ...
  • 333
3 votes
1 answer
86 views

Saving the values of all variabled declared in a block

Often when running exprimental code to empirically investigate an algorithm, you tweak a parameter, then rerun the code, and generate a bunch of saved models for various parameters. To make sure I ...
23 votes
5 answers
2k views

Benchmarks of various scientific programming languages for theoretical modelling

For a theoretical modelling course for biology students, I am trying to decide which would be the best technical programming language for doing evolutionary simulations in terms of elegance and ...
4 votes
1 answer
1k views

Speeding up OpenCL matrix-vector multiplication

So I'd like to get a performance boost above and beyond standard Julia matrix-vector multiply, using my Intel HD Graphics 4000 1536 MB GPU, but I can't do better than an order of magnitude worse ...
  • 362
2 votes
1 answer
190 views

Recursively walk directories with a callback function

This code works, but I want to refine it: ...
1 vote
0 answers
85 views

Using RemoteChannels for an parallel iterable map

RemoteChannels hit with julia v0.5 There are no examples online, but I managed to get them to work. They are still very new, and all this whole area of multithreading code continues to have running ...
4 votes
1 answer
1k views

Sort Algorithms in Julia

Similar to this question here, I am trying to implement sort algorithms in Julia, as part of a package. The code for the implementation is as follows: ...
  • 295
3 votes
2 answers
994 views

Search algorithms in julia

I've been trying to implement the basic search algorithms: Sequential search for ordered and unordered arrays and the binary search algorithms, as part of a package. So, this is the implementation: <...
  • 295
4 votes
2 answers
223 views

Computing the double Integral using MonteCarlo techniques using Julia

I decided to try and learn Julia for doing scientific computing, and I decided to tackle the problem of finding $$ \int_{D_{\frac{1}{4}}} x^4 + y^2 dA $$ where \$ D_{\frac{1}{4}} \$ is the part of ...
5 votes
1 answer
162 views

Implementing cat function in Julia

I wanted to implement the cat function in Julia and I ran a few test cases but I am not sure if I am missing anything: ...
1 vote
0 answers
493 views

Flattening lists of Vectors and Matrices into arrays

For many optimisation methods the parameters must be given as a single vector, such as Optim.jl. But for many interesting functions the parameters/coeffients to be optimised are matrices, or a mixture ...
4 votes
2 answers
180 views

Latent Dirichlet Allocation Posterior Inference

I have implemented Latent Dirichlet Allocation in Julia. Since the code is rather long, I copied the essential part here but the complete code can be found on GitHub. ...
  • 53
5 votes
0 answers
159 views

Kernel density regression in Julia

Here is a naive implementation of a kernel density regression algorithm in Julia. The code works, but there is quite a bit of room for improvement though. First, there is not an optimal bandwidth ...
  • 551
2 votes
0 answers
34 views

Iterative procedure to get 3D coordinates from distance constraints

Imagine you have a series of n points randomly generated in a box in 3D space. You also have a list of distance bounds, e.g. points 5 and 3 should be between 1.0 and 2.0 Angstroms apart. There are ...
  • 153
5 votes
1 answer
535 views

Terry Feagin's 10th order explicit Runge-Kutta method

The following Julia code implements Terry Feagin's 10th order explicit Runge-Kutta method (a more accurate cousin of RK4). Though the structure of the code is quite simple (i.e. no cyclomatic ...
2 votes
1 answer
735 views

Functions for outputting the rows and columns of a matrix

I often want to iterate either row-wise or column-wise though a matrix in Julia, so I've created a pair of functions to help: ...
8 votes
1 answer
226 views

Individuals reproduce and mutate

How can I improve the performance (in terms of computation time) of this code? ...
  • 617
6 votes
2 answers
706 views

Julia RPN calculator algorithm ported from Python version

I've just finished the Reverse Polish Notation task in Julia from rosettacode.com by porting the existing Python code. There are only a few lines I don't get how to port: ...
4 votes
1 answer
378 views

Neural Network in Julia

I am currently trying to implement a Neural Net in Julia with the goal of eventually implementing a stacked autoencoder. My code seems to work but I would appreciate any constructive criticism. If ...
2 votes
1 answer
246 views

Speed optimization for girth-counting algorithm

I have the following girth-counting algorithm on a large 0-1 matrix and I need to optimize it for speed: ...
  • 129
12 votes
2 answers
1k views

Neural Network appears to be slow

I have some neural network Julia code which I'm hoping to speed up. It's possibly simply poorly-designed, but I'm not convinced. ...
  • 362