I have the following code: if sol[5,i] - sol_1[5,i] >= 0.0 s[i] = sol[5,i] - sol_1[5,i] else s[i] = sol_1[5,i] - sol[5,i] end end Now, here, sol and sol_1 correspond to 2 solutions I have from DifferentialEquations.jl (basically, I would like to plot the absolute value of the difference at each timestep; my solution vectors are same sized and have 6 variables - here, I am only considering the 5th variable). Now, problem with this is that the code doesn’t work. I get both positive and negative values for s which certainly not what I intended. In fact, I ran the following line, s = sol[5,i] - sol_1[5,i] This and the above code give me the same s. Any help would be appreciated!

Can you explain why you are using if/else instead of abs? It seems like a very particular, but odd, choice.

[image] Sh_166: the absolute value of the difference s = abs(sol[5,i] - sol_1[5,i]) You haven’t provided enough information for anyone to tell why your original code doesn’t work (suggested in the link in the other response). Maybe there’s something else in your loop overwriting s[i]. Sometimes unexpected branching can be due to comparing with NaN, but if that were the case you’d likely be complaining about s getting filled with NaN.

Ah, sorry. My code is quite large, so I skipped it (it takes time to run, and contains a lot of irrelevant portions). However, I did manage to figure it out (i.e. I managed to compute what I wanted), so there’s that.

Actually, I did at first. Must’ve made a typo there because it kept showing me error, so I went a more naive way. Anyhow, using abs.() works fine now, so I have just skipped the loop. Thanks a lot.

If-else loop failure for a large array

General Usage

stevengj April 6, 2023, 2:00pm 2

See Please read: make it easier to help you — it’s hard to help you without a minimal code that we can actually run to reproduce your problem.

Topic		Replies	Views
Callback function in DifferentialEquations General Usage question , diffeq	1	327	February 15, 2021
This if statements is not work properly New to Julia question	2	484	January 6, 2020
Algorithm with if...else New to Julia question	8	1279	June 16, 2021
Automatic differentiation and `abs()` (or comparison operators in general) issues with JuMP/Ipopt General Usage	4	736	October 1, 2019
Submit a bug report while trying to use DifferentialEquations.jl General Usage windows , bug	1	400	May 31, 2024

If-else loop failure for a large array

Related topics