Function transformTriangle

plugin/seq/plotPDF.cpp:1874–1899

    if( fabs(lambda2) > EPS )
	F -= E*E/(4*lambda2);

    return;
}

void transformTriangle( double *const Vx, double *const Vy,
			const double *const vx, const double *const vy, const double *const PHI )
{
    const double &lambda1 = PHI[0]; const double &lambda2 = PHI[1];
    const double &D = PHI[6]; const double &E = PHI[7];

    const double &ev1x = PHI[2]; const double &ev1y = PHI[3];
    const double &ev2x = PHI[4]; const double &ev2y = PHI[5];
    const double  P[2][2] = { { ev1x, ev2x }, { ev1y, ev2y } };
    const double PT[2][2] = { { P[0][0], P[1][0] }, { P[0][1], P[1][1] } };

    const int nVertices = 3;
    
    const double EPS = 1e-10;

    for(int i = 0; i < nVertices; i++){
	Vx[i] = PT[0][0] * vx[i] + PT[0][1] * vy[i];
	if( fabs(lambda1) > EPS )
	    Vx[i] += D/(2*lambda1);
	Vy[i] = PT[1][0] * vx[i] + PT[1][1] * vy[i];
	if( fabs(lambda2) > EPS )
	    Vy[i] += E/(2*lambda2);
    }