336         real*8  tmp1,   tmp2,   tmp3,   tmp4  local
350         tmp1 = 0
352             tmp1 = tmp1 + x(i,1) * y(i)
354         c(1) = tmp1
358         tmp1 = 0
361             tmp1 = tmp1 + x(i,1) * y(i)
364         c(1) = tmp1
369         tmp1 = 0
373             tmp1 = tmp1 + x(i,1) * y(i)
377         c(1) = tmp1
383         tmp1 = 0
388             tmp1 = tmp1 + x(i,1) * y(i)
393         c(1) = tmp1
400         tmp1 = 0
406             tmp1 = tmp1 + x(i,1) * y(i)
412         c(1) = tmp1
420         tmp1 = 0
427             tmp1 = tmp1 + x(i,1) * y(i)
434         c(1) = tmp1
443         tmp1 = 0
451             tmp1 = tmp1 + x(i,1) * y(i)
459         c(1) = tmp1
473             tmp1 = 0
482                 tmp1 = tmp1 + x(i,j+0) * y(i)
491             c(j+0) = tmp1
518         real*8  tmp1,   tmp2,   tmp3,   tmp4  local
531         tmp1 = c(1)
534      1           + tmp1 * x(i,1)
539         tmp1 = c(1)
543      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
548         tmp1 = c(1)
554      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
560         tmp1 = c(1)
566      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
572         tmp1 = c(1)
579      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
586         tmp1 = c(1)
594      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
601         tmp1 = c(1)
610      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
622             tmp1 = c(j+0)
632      1               + tmp1 * x(i,j+0) + tmp2 * x(i,j+1)
656         real*8  tmp1,   tmp2,   tmp3,   tmp4  local
669         tmp1 = c(1)
672      1           + tmp1 * x(i,1)
677         tmp1 = c(1)
681      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
686         tmp1 = c(1)
691      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
697         tmp1 = c(1)
703      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
709         tmp1 = c(1)
716      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
723         tmp1 = c(1)
731      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
738         tmp1 = c(1)
747      1           + tmp1 * x(i,1) + tmp2 * x(i,2)
765             tmp1 = c(j+0)
775      1               + tmp1 * x(i,j+0) + tmp2 * x(i,j+1)
799         real*8  tmp1,   tmp2,   tmp3,   tmp4  local
812         tmp1 = c(1)
815      1           - tmp1 * x(i,1)
820         tmp1 = c(1)
824      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
829         tmp1 = c(1)
834      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
840         tmp1 = c(1)
846      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
852         tmp1 = c(1)
859      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
866         tmp1 = c(1)
874      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
881         tmp1 = c(1)
890      1           - tmp1 * x(i,1) - tmp2 * x(i,2)
902             tmp1 = c(j+0)
912      1               - tmp1 * x(i,j+0) - tmp2 * x(i,j+1)