Make sure init is always inside the domain

src/Cubes/Rechunker.jl

@@ -57,16 +57,19 @@ function get_copy_buffer_size(incube, outcube;writefac=4.0, maxbuf = YAXDefaults
     if nd == 1
         return (min(maxbuf,length(incube)),)
     end
-    init = fill(maxbuf^(1/nd),nd-1)
-    incs = DiskArrays.approx_chunksize(eachchunk(incube))
-    outcs = DiskArrays.approx_chunksize(eachchunk(outcube))
     insize = size(incube)
     outsize = size(outcube)
+    totsize = prod(insize)
+    incs = DiskArrays.approx_chunksize(eachchunk(incube))
+    outcs = DiskArrays.approx_chunksize(eachchunk(outcube))
     #Catch case where buffer is larger than cube
     if maxbuf > prod(insize)
         return insize
     end
+    init = [(Float64(maxbuf)*insize[i]^(nd)/totsize)^(1/nd) for i in 1:(nd-1)]
+    @debug "Init: ", init
     r = optimize(sz->optifunc(sz,maxbuf,incs,outcs,insize, outsize,writefac),init, iterations=100)
+    @debug "Optimization result: $(r.minimizer)"
     bufnow = (r.minimizer...,maxbuf/prod(r.minimizer))
 
     bufcorrected = if align_output
@@ -91,6 +94,9 @@ end
 function copy_diskarray(incube,outcube;writefac=4.0, maxbuf = YAXDefaults.max_cache[], align_output=true)
     size(incube) == size(outcube) || throw(ArgumentError("Input and output cubes must have the same size"))
     bufcorrected = get_copy_buffer_size(incube, outcube;writefac,maxbuf,align_output)
+    @debug "Copying with buffer size $bufcorrected"
+    @debug "Input chunk size: $(eachchunk(incube).chunks)"
+    @debug "Output chunk size: $(eachchunk(outcube).chunks)"
     copybuf = DiskArrays.GridChunks(size(outcube),bufcorrected)
     copydata(outcube,incube,copybuf)
 end