The implementation of a variational quantum eigensolver, VQE, is simple: Provided with an optimization routine define the loss function as the expectation value of your Hamiltonian with the ansatz and plug it into the optimizer. Our current implementation of the VQE doesn't necessarily reflect this.
For research applications it can be easier to re-implement it from scratch rather than use the