`OPRD`,*result*,*matrix*,*orb1*,*orb2*,*factor*

Takes the column vectors *v1* and *v2* from *matrix*
and adds their outer product to *result*. *v1* and *v2*
must be given in the form *icol.isym*, e.g., 3.2 means the third
vector in symmetry 2. The result is

If *result* has not been used before, it is zeroed before performing
the operation.

molpro@molpro.net 2018-09-19