As part of my code for julia-1.6 I try to calculate an NxN symmetric ‘distance’ matrix, where ‘distance’ is a function of 4 vectors (extracted as columns from a DataFrame),and two parameters (‘R2’ and ‘range’). While the code below worked with julia-1.0, it now gets stuck without an error message. …

What is N? What is d2? What is z_pred? I’ll stop there, my point is basically make an MWE so that people can reproduce your issues, and ideally read: [image] Please read: make it easier to help you Meta Discussion Welcome to the Julia Discourse! We are enthusiastic about …

Sorry for not putting my question in the right way. I’ll try to do it better now by providing a simplified version of the problematic part of my code. See below. # A simplified version of my calculation of a distance matrix 'd2': N = 5 d2 = zeros(N,N) R2 = 0.36 range = 582 z_pred = (1, 8, 3, 2, 4) …

[image] jjdegruijter: To my surprise this small example worked and produced a good 5x5 matrix. So I’ve still no idea why the real code doesn’t work. This is the reason MWE’s are so important, because they teach the person afflicted by the problem that their problem is not where they think it…

I wouldn’t be that harsh. But MWEs are great! They make your questions better, and often solve the problem outright. I’ll add one piece of advice: wrap your code in a function. It makes the code more self-contained and performant, and removes problems caused by various hidden state you may have han…

Do you mean all the code in one function? What I have now is 4 files, say A, B, C and D. The codes in B, C and D are each wrapped in a function. Any run follows A > B > C or A > B > D, depending on a parameter setting in A. This is the reason why I did not see how to make a MWE

I meant put the main part of your MWE in a function.

[image] jjdegruijter: So I’ve still no idea why the real code doesn’t work. You could still show us the error message at least?

As others have told, without further information is difficult to know with certainty, but let’s try a guess: If the error really happens in the code you have posted, it may be either that: The vectors used inside the loop cannot be indexed with the i, j values that are defined (e.g. they have les…

Part of the problem is that there isn’t any error message. Julia simply halts

The way to debug something like this can be: Be sure that you don’t have any global variable anywhere in your functions, particularly that one that halts. This is a good advice not only because of the peculiarity of the Julia performance that depends on this, but as a general advice in any progra…

How to upgrade code for a symmetric distance matrix from julia-1.0 to julia-1.6

General Usage

nilshg March 2, 2022, 2:07pm 2

What is N? What is d2? What is z_pred? I’ll stop there, my point is basically make an MWE so that people can reproduce your issues, and ideally read:

Topic		Replies	Views
Distance matrix General Usage	5	332	July 22, 2022
Looping in Julia New to Julia question	6	574	September 2, 2022
Symmetric matrix is not symmetric General Usage question , linearalgebra , distances	2	792	February 17, 2021
Pairwise distances from a single column or vector New to Julia	7	4182	October 2, 2019
Improving performance of a nested for loop New to Julia performance , loops	14	6180	October 10, 2019

How to upgrade code for a symmetric distance matrix from julia-1.0 to julia-1.6

Related topics