Function findFillValue

plugin/seq/plotPDF.cpp:2912–2939

	    break;
	}
    }

    return;
}

double findFillValue( const std::vector<double> &cx, const std::vector<double> &cy, const double *const phi )
{
    const double &a = phi[0]; const double &b = phi[1]; const double &c = phi[2];
    const double &d = phi[3]; const double &e = phi[4]; const double &f = phi[5];

    int npoints = 0;

    // find average of phi at nodal points
    double sum = 0;
    for(size_t i = 0; i < cx.size(); i += 3){

	const double &x = cx[i];
	const double &y = cy[i];

	sum += (a*x*x + b*x*y + c*y*y + d*x + e*y + f);
	npoints++;

	if( (i+3 < cx.size()) && isSegment(cx,cy,i) ){

	    const double mx = ( cx[i]+cx[i+3] ) / 2;
	    const double my = ( cy[i]+cy[i+3] ) / 2;
	    sum += (a*mx*mx + b*mx*my + c*my*my + d*mx + e*my + f);
	    npoints++;
	}
    }