124 FR_DOUBLE qsin,qcos,sincos,sigr,sigt,taurt,sigv;  in scsh88()  local
344             qcos= cosphi*cosphi;  in scsh88()
347             sigr= sig[1]*qcos + sig[2]*qsin + 2.*sig[3]*sincos;  in scsh88()
348             sigt= sig[1]*qsin + sig[2]*qcos - 2.*sig[3]*sincos;  in scsh88()
349             taurt=(sig[2]-sig[1])*sincos + sig[3]*(qcos-qsin);  in scsh88()
387             qcos= cosphi*cosphi;  in scsh88()
390             sigr= sig[1]*qcos + sig[2]*qsin + 2.*sig[3]*sincos;  in scsh88()
391             sigt= sig[1]*qsin + sig[2]*qcos - 2.*sig[3]*sincos;  in scsh88()
392             taurt=(sig[2]-sig[1])*sincos + sig[3]*(qcos-qsin);  in scsh88()
596           qcos= cosphi*cosphi;  in scsh88()
599           sigr= sig[1]*qcos + sig[2]*qsin + 2.*sig[3]*sincos;  in scsh88()
600           sigt= sig[1]*qsin + sig[2]*qcos - 2.*sig[3]*sincos;  in scsh88()
601           taurt=(sig[2]-sig[1])*sincos + sig[3]*(qcos-qsin);  in scsh88()