🤖 Format .jl files

Author: amontoison

src/ReadFiles.jl

@@ -18,41 +18,41 @@ Expected file structure:
   - 3D points: 3 lines per point
 """
 function readfile(filename::String; T::Type = Float64)
-    open(filename) do io
-        f = Bzip2DecompressorStream(io)
-
-        # First line: number of cameras, points, and observations
-        ncams, npnts, nobs = map(x -> parse(Int, x), split(readline(f)))
-        @debug "$filename: reading $ncams cameras, $npnts points, $nobs observations"
-
-        # Preallocate vectors
-        cam_indices = Vector{Int}(undef, nobs)
-        pnt_indices = Vector{Int}(undef, nobs)
-        pt2d = Vector{T}(undef, 2 * nobs)
-        x0 = Vector{T}(undef, 3 * npnts + 9 * ncams)
-
-        # Read observations
-        for i = 1:nobs
-            cam, pnt, x, y = split(readline(f))
-            cam_indices[i] = parse(Int, cam) + 1  # Convert to 1-based indexing
-            pnt_indices[i] = parse(Int, pnt) + 1
-            pt2d[2 * i - 1] = parse(T, x)
-            pt2d[2 * i] = parse(T, y)
-        end
-
-        # Read camera parameters (9 values per camera)
-        for i = 1:ncams
-            offset = 3 * npnts + 9 * (i - 1)
-            for j = 1:9
-                x0[offset + j] = parse(T, readline(f))
-            end
-        end
-
-        # Read 3D points (3 coordinates per point)
-        for k = 1:(3 * npnts)
-            x0[k] = parse(T, readline(f))
-        end
-
-        return cam_indices, pnt_indices, pt2d, x0, ncams, npnts, nobs
+  open(filename) do io
+    f = Bzip2DecompressorStream(io)
+
+    # First line: number of cameras, points, and observations
+    ncams, npnts, nobs = map(x -> parse(Int, x), split(readline(f)))
+    @debug "$filename: reading $ncams cameras, $npnts points, $nobs observations"
+
+    # Preallocate vectors
+    cam_indices = Vector{Int}(undef, nobs)
+    pnt_indices = Vector{Int}(undef, nobs)
+    pt2d = Vector{T}(undef, 2 * nobs)
+    x0 = Vector{T}(undef, 3 * npnts + 9 * ncams)
+
+    # Read observations
+    for i = 1:nobs
+      cam, pnt, x, y = split(readline(f))
+      cam_indices[i] = parse(Int, cam) + 1  # Convert to 1-based indexing
+      pnt_indices[i] = parse(Int, pnt) + 1
+      pt2d[2 * i - 1] = parse(T, x)
+      pt2d[2 * i] = parse(T, y)
     end
+
+    # Read camera parameters (9 values per camera)
+    for i = 1:ncams
+      offset = 3 * npnts + 9 * (i - 1)
+      for j = 1:9
+        x0[offset + j] = parse(T, readline(f))
+      end
+    end
+
+    # Read 3D points (3 coordinates per point)
+    for k = 1:(3 * npnts)
+      x0[k] = parse(T, readline(f))
+    end
+
+    return cam_indices, pnt_indices, pt2d, x0, ncams, npnts, nobs
+  end
 end

test/testBundleAdjustmentAllocations.jl

@@ -6,28 +6,27 @@ meta_nls = nls_meta(model)
 S = typeof(model.meta.x0)
 
 @testset "residual allocations" begin
-    F = S(undef, meta_nls.nequ)
+  F = S(undef, meta_nls.nequ)
 
-    residual!(model, model.meta.x0, F)
-    residual_alloc(model) = @allocated residual!(model, model.meta.x0, F)
-    @test residual_alloc(model) == 0
+  residual!(model, model.meta.x0, F)
+  residual_alloc(model) = @allocated residual!(model, model.meta.x0, F)
+  @test residual_alloc(model) == 0
 end
 
 @testset "jac_structure_residual allocations" begin
-    rows = Vector{Int}(undef, meta_nls.nnzj)
-    cols = Vector{Int}(undef, meta_nls.nnzj)
+  rows = Vector{Int}(undef, meta_nls.nnzj)
+  cols = Vector{Int}(undef, meta_nls.nnzj)
 
-    jac_structure_residual!(model, rows, cols)
-    jac_structure_residual_alloc(model, rows, cols) =
-        @allocated jac_structure_residual!(model, rows, cols)
-    @test jac_structure_residual_alloc(model, rows, cols) == 0
+  jac_structure_residual!(model, rows, cols)
+  jac_structure_residual_alloc(model, rows, cols) =
+    @allocated jac_structure_residual!(model, rows, cols)
+  @test jac_structure_residual_alloc(model, rows, cols) == 0
 end
 
 @testset "jac_coord_residual allocations" begin
-    vals = S(undef, meta_nls.nnzj)
+  vals = S(undef, meta_nls.nnzj)
 
-    jac_coord_residual!(model, model.meta.x0, vals)
-    jac_coord_residual_alloc(model, vals) =
-        @allocated jac_coord_residual!(model, model.meta.x0, vals)
-    @test jac_coord_residual_alloc(model, vals) == 0
+  jac_coord_residual!(model, model.meta.x0, vals)
+  jac_coord_residual_alloc(model, vals) = @allocated jac_coord_residual!(model, model.meta.x0, vals)
+  @test jac_coord_residual_alloc(model, vals) == 0
 end