Closed aniketsen closed 1 year ago
@aniketsen Can you change the source of this pull-request to the quda_work_eigensolver branch of etmc/tmLQCD (which I just created based on aniketsen:quda_work_eigensolver)? That way we can all push into it for minor tweaks and changes. (I wanted to add some documentation, for example)
I've now also merged with the latest quda_work
there as aniketset:quda_work_eigensolver
was not kept up to date with etmc:quda_work
and was missing a number of commits (several of which related to the continuous integration stuff in .github)
eigsolveQuda
sets the parameter and calls the eigensolver in QUDA.initQudaforEig
initializes QUDA. This currently only works for theTwoFalvourSolver
, but should be easy to generalize. Ineigenvalues_bi.cpp
, the interface is called ifTM_USE_QUDA
is defined, probably not ideal and can be changed to a different flag.Currently the interface passes
NULL
for the eigenvectors, since they are not necessary. On QUDA side, support is added to deal with this.