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
0 answers
28 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 ...
user avatar
4 votes
1 answer
171 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 ...
user avatar
0 votes
0 answers
29 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 ...
user avatar
  • 215
2 votes
0 answers
39 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 ...
user avatar
3 votes
3 answers
130 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. ...
user avatar
2 votes
2 answers
138 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. ...
user avatar
  • 333
1 vote
0 answers
77 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: ...
user avatar
  • 455
2 votes
1 answer
99 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 ...
user avatar
  • 143
5 votes
2 answers
178 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 ...
user avatar
  • 247
3 votes
1 answer
128 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 ...
user avatar
  • 31
3 votes
1 answer
132 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: ...
user avatar
2 votes
1 answer
70 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 ...
user avatar
  • 195
6 votes
1 answer
239 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 ...
user avatar
  • 225
3 votes
2 answers
213 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....
user avatar
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 ...
user avatar
  • 221
6 votes
0 answers
131 views

Julia set in Julia (and other fractals)

I've made the following to draw a series of zooming in images of a fractal. ...
user avatar
  • 3,512
2 votes
0 answers
71 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 '...
user avatar
  • 121
6 votes
0 answers
125 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 ...
user avatar
  • 161
6 votes
1 answer
663 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). ...
user avatar
  • 161
2 votes
1 answer
47 views

Add optional distance range to function

I have a simple brute-force search function: ...
user avatar
  • 1,423
2 votes
1 answer
171 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 ...
user avatar
  • 1,423
1 vote
1 answer
71 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 ...
user avatar
  • 111
0 votes
1 answer
80 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 ...
user avatar
  • 109
1 vote
0 answers
204 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 ...
user avatar
  • 11
3 votes
1 answer
66 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 ...
user avatar
  • 832
2 votes
2 answers
149 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 ...
user avatar
3 votes
0 answers
252 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 ...
user avatar
  • 131
3 votes
0 answers
662 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 ...
user avatar
1 vote
0 answers
123 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 ...
user avatar
  • 333
3 votes
1 answer
81 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 ...
user avatar
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 ...
user avatar
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 ...
user avatar
  • 362
2 votes
1 answer
188 views

Recursively walk directories with a callback function

This code works, but I want to refine it: ...
user avatar
1 vote
0 answers
83 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 ...
user avatar
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: ...
user avatar
  • 295
3 votes
2 answers
967 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: <...
user avatar
  • 295
4 votes
2 answers
220 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 ...
user avatar
5 votes
1 answer
158 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: ...
user avatar
1 vote
0 answers
489 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 ...
user avatar
4 votes
2 answers
178 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. ...
user avatar
  • 53
5 votes
0 answers
157 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 ...
user avatar
  • 551
2 votes
0 answers
33 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 ...
user avatar
  • 153
5 votes
1 answer
525 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 ...
user avatar
2 votes
1 answer
731 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: ...
user avatar
8 votes
1 answer
224 views

Individuals reproduce and mutate

How can I improve the performance (in terms of computation time) of this code? ...
user avatar
  • 607
6 votes
2 answers
702 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: ...
user avatar
4 votes
1 answer
376 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 ...
user avatar
2 votes
1 answer
245 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: ...
user avatar
  • 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. ...
user avatar
  • 362