apply_Dense_Adiabatic_base Subroutine

    subroutine apply_Dense_Adiabatic_base(A, b, isUniform, Edge, Dval, perm)
        implicit none
        real(real64), intent(inout), optional :: A(:, :)
        real(real64), intent(inout) :: b(:)
        logical(logical32), intent(in) :: isUniform
        integer(int32), intent(in) :: Edge(2)
        real(real64), intent(in) :: Dval
        integer(int32), intent(in) :: perm(:)

        integer(int32) :: i, ind, ps, pe
        integer(int32) :: p1, p2

        ! if (isUniform) then
        !     p1 = perm(Edge(1))
        !     p2 = perm(Edge(2))

        !     if (present(A)) then
        !         A(p1, :) = 0.0d0
        !         A(p1, p1) = 1.0d0

        !         A(p2, :) = 0.0d0
        !         A(p2, p2) = 1.0d0
        !     end if

        !     b(p1) = Dval
        !     b(p2) = Dval
        ! end if

    end subroutine apply_Dense_Adiabatic_base