friend pvector O(const pvector& in1) friend pvector Linv(const pvector& in1)Return value:
- pvector O: Application of overlap matrix to in1
- pvector Linv: Application of inverse of L to in1
Input:
- in1: pvector
Uncomment the function L() in pvector.C and use it as a guide to
construct functions of the above prototypes, recalling that 
, where V is the volume of the cell and 
is the identity matrix, and that Linv divides all components by
, except for 
, which it multiplies by zero.
Note that the code provided for L() already takes proper account of
aliasing, as discussed in lecture. To debug, uncomment the
tests for O and L, and run the program.