Lines Matching refs:a1
251 double a0,a1,b0,b1;// coeffs in a[k] and b[k]
256 // a1*x(1)+b1*y(1)=r1 & 2*a1*a(1)+2*b1*b(1)=rr1
259 a1 = r1/dot(v1,V.at1())*v1[0]-rr1/2*v1[1];
262 a[k] = Linear(a0,a1);
557 * Below are basic functions dedicated to solving this assuming a0 and a1 !=0.
564 double a0=aa0,a1=aa1,c0=cc0,c1=cc1;
566 bool flip = a1<0;
567 if (a1<0){a1=-a1; c1=-c1;}
569 double a = (a0<a1 ? a0 : a1);
590 solve_lambda0(double a0,double a1,double c0,double c1,
594 p[0] = Linear( a1*c0*c0+c1, a1*a0*(a0+ 2*c0) +a1*c0*c0 +c1 -1 );
595 p[1] = Linear( -a1*a0*(a0+2*c0), -a1*a0*(3*a0+2*c0) );
596 p[2] = Linear( a1*a0*a0 );
598 OptInterval domain = find_bounds_for_lambda0(a0,a1,c0,c1,insist_on_speeds_signs);
655 // lambda0 = a1 lambda1^2 + c1
656 double a0,c0,a1,c1;
659 a1 = -d2M1xdM1/2/dM1xdM0;
664 lambda0.push_back( a1*c0*c0 + c1 );
665 }else if (fabs(a1)<epsilon){
670 vector<double> solns=solve_lambda0(a0,a1,c0,c1,insist_on_speed_signs);