added lapeq.h, computes A+=alpha*B