Fix an array bound mismatch.

Description

Fix an array bound mismatch.

In the case of k-point symmetries, together with a k-point path, the copy routine was leading to an array bound mismatch, as the arrays of the original grids are allocated to the number of points in the full Brillouin zone, and not to the irreducible one.

The range of the DOS was also computed taking into account k-points in the path, which is now solved.

News snippet

Fix an array bound mismatch. The range of the DOS is now independent of the k-points in the path.

Checklist

• I have checked that my code follows the Octopus coding standards
• I have added tests for all the new features added in this request.

src/system/dos.F90

@@ -208,7 +208,6 @@ contains
           ! sum up Lorentzians
           do ik = 1, st%d%nik, ns
             do is = 0, ns-1
-             print *, ik,  is , st%d%kweights(ik+is)
              dos(ie, ist, is) = dos(ie, ist, is) + st%d%kweights(ik+is) * M_ONE/M_Pi * &
                 this%gamma / ( (energy - st%eigenval(ist, ik+is))**2 + this%gamma**2 )
             end do