submitted by /u/heyheyhey27 [visit reddit] [comments ..read more

I'm trying to run the following code for this practice question Solve this problem using value function iteration, for x ∈ [10−4, 10]. Use an evenly spaced grid for x with 50 points. Start with an initial guess of v0(x) = 0, and use linear interpolation to approximate your guess so that you can evaluate it off of the grid. Keep iterating until the value functions stop changing within a tolerance of 10−6.3 Use three different values for β: 0.9, 0.95, 0.99. Plot the log errors against the number of iterations for each solution. What are the slopes? Interpret these results. using Interpolation ..read more

Hi , I am new to julia and I was wondering how do I find documentation for packages to know how to apply functions. I checked julia documentation but I am unable to find documentation for functions such as Roots.jl , Optim.jl , Parameters.jl and such would any one be able to tell me how to go about this. Thanks! submitted by /u/SquareZebra610 [visit reddit] [comments ..read more

submitted by /u/QuasarInfection [visit reddit] [comments ..read more

TL;DR Great community. Excellent at expressing "math". Very fast language. Almost great REPL. Immature ecosystem. Inconvenient debugging. Bad code organization. Background Recently I finished a decently sized Julia project (~3100 LoC), and I'd like to share my experience of using the language. It's mostly about developer experience, so I hope these will provide insights for language developers and users. Before this project, I've already used Julia plenty of times, mostly to analyze some experiment data generated by other programs. They are generally fairly clean, because I want the data to be ..read more

Grassmann.jl v0.8.14 https://github.com/chakravala/Grassmann.jl Note that Quaternions in Grassmann.jl are oriented with k = -v13 ```Julia julia> using Grassmann; basis"3" (⟨×××⟩, v, v₁, v₂, v₃, v₁₂, v₁₃, v₂₃, v₁₂₃) julia> s,i,j,k = v,v12,-v13,v23 (v, v₁₂, -1v₁₃, v₂₃) julia> i*j == k true julia> quatvalues(s+2i+3j+4k) 4-element StaticVectors.Values{4, Int64} with indices SOneTo(4): 1 2 3 4 julia> v1 ⊘ exp(π*i/2) -1.0v₁ + 1.2246467991473532e-16v₂ + 0.0v₃ julia> v1 ⊘ exp(π*i/4) 2.220446049250313e-16v₁ + 1.0v₂ + 0.0v₃ julia> v1 ⊘ exp(π*i/8) 0.7071067811865475v₁ + 0.7071067811 ..read more

It has been mentioned in a talk at the juliacon that for obtaining the first image of a black hole they had to use 1 week of c++ compute on a cluster and using julia and differentiable programming someone got this down to 1h on 1 thread on a laptop. See https://youtube.com/clip/Ugkxd34TdNe4uzAdOEVPkg0-Hc-kF1t-RfIn?si=zVi1yZ0HQxLdUXqi This is pretty crazy to me. Does someone know any details? submitted by /u/AlrikBunseheimer [visit reddit] [comments ..read more

Hello, So I have just installed julia 1.10.2 on a mac and I try running A = [1 2 3;4 5 6;7 8 9] inv(A) and I expect a singular error. However, I actually get a matrix as an output. What is weirder is that if I do the same on a Windows or a Linux machine running arch, I do indeed get an error. Is this something to do with Mac or is there a problem with how I am running my code? For reference, I installed julia 1.6.4 on mac and the same code does in fact throw an error. Is there any intuition behind this? submitted by /u/Swagdalfthegrey [visit reddit] [comments ..read more

In short they released TidierDB.jl which brings TidierData.jl/Tidier.jl (tidyverse) syntax to databases. The default backend is Duckdb so u can read in arrow/parquet/json/csvs directly, but it also supports 6 others (clickhouse, postgres, mysql etc). This is similar to the Ibis python package and dbplyr R package, for those familiar. I played around it with, and its pretty slick. Perhaps some of you saw this thru twitter/discourse, but I thought I'd share here since the Tidier.jl team isn't on the reddit it seems. ​ submitted by /u/Suspicious-Oil6672 [visit reddit] [comments ..read more

Howdy! I'm trying to run a mixed effects model in Julia, and the one question I have is about convergence issues/warnings. In R, the lme4 and lmerTest packages will give convergence warnings for singularity and degenerate Hessians. I've had instances where the same model throws a degenerate Hessian warning in R but nothing in Julia. Now, I'm pretty new to models, and my advisor is very used to R (I'm mostly just poking around out of interest), so I guess I'm not sure if the model does converge in Julia but for some reason R's implementation of the BOBYQA or NM optimizers doesn't? I'm hoping s ..read more