typedef double l_double;  //precision for lanczos