multiply_hamil_and_vector_full_complex Subroutine

    subroutine multiply_hamil_and_vector_full_complex(this, input_vector, output_vector)
        type(HamiltonianCalcType), intent(inout) :: this
        complex(dp), intent(in) :: input_vector(:)
        complex(dp), intent(out) :: output_vector(:)

        associate(space_size => this%space_size)

            call zgemv('N', &
                       space_size, &
                       space_size, &
                       1.0_dp, &
                       hamiltonian, &
                       space_size, &
                       input_vector, &
                       1, &
                       0.0_dp, &
                       output_vector, &
                       1)

        end associate
    end subroutine multiply_hamil_and_vector_full_complex