@marc and @xavier and anyone else courageous...
projdb is a horrible routine. I hope it is optimized and well benchmarked, in which case we can leave its internals as is, but we should:
1) break it into 2 routines, to calculated scprod and apply the orthogonalization, respectively, avoiding having intent(inout) variables and huge useless if() statements
2) reuse a _lot_ of code - this is absurd
3) document what the ortalg flag does.
cheers!
Matthieu
projdb cleaning.
projdb cleaning.
Matthieu Verstraete
University of Liege, Belgium
University of Liege, Belgium