Function plot_P1_isoline_body

plugin/seq/plotPDF.cpp:1360–1419

    }

    return;
}

void plot_P1_isoline_body( std::stringstream &Content, const Fem2D::Mesh &Th, const KN<double> &f_P1,
			   const std::vector<double> &isoline_val, const double cmin, const double cmax,
			   const KNM<double> &palette, const double scale, const double ar,
			   const double x0, const double y0,
			   const int marginl, const int marginb,
			   const bool monochrome, const bool logscale,
			   const int NISOLINES, const double linewidth )
{
    std::stringstream &st = Content;

    st << "q\n";
    st << linewidth << " w\n"; // setlinewidth
    st << "1 0 0 1 " << PLOTPDFVAR::PADDING+marginl << " " << PLOTPDFVAR::PADDING+marginb << " cm\n";

    const int &nTriangles = Th.nt;
    for(int k = 0; k < nTriangles; k++){

	const int &v0 = Th(k,0);
	const int &v1 = Th(k,1);
	const int &v2 = Th(k,2);

	const double vx[] = { Th(v0).x,     Th(v1).x,     Th(v2).x };
	const double vy[] = { Th(v0).y,     Th(v1).y,     Th(v2).y };
	const double vf[] = { f_P1[3*k+0],  f_P1[3*k+1],  f_P1[3*k+2] };

	for(size_t m = 0; m < isoline_val.size(); m++){

	    const double &value = isoline_val[m];

	    std::vector<double> px, py;
	    trackP1isoline( px, py, vx, vy, value, vf );

	    assert( px.size() == py.size() );

	    if( px.size() == 0 ) continue; // goto next isoline_val

	    setrgbcolor(st, value, palette, cmin, cmax, monochrome, logscale);

	    if( px.size() > 3 ){ 

		// f(x,y)==value on all vertices, i.e. f(x,y) \equiv value on the triangle
		st << "rg\n";
		st << scale*ar*(vx[0] - x0) << ' ' << scale*(vy[0] - y0) << " m "
		   << scale*ar*(vx[1] - x0) << ' ' << scale*(vy[1] - y0) << " l "
		   << scale*ar*(vx[2] - x0) << ' ' << scale*(vy[2] - y0) << " l f\n";
		
	    } else {

		// assert( (px.size() == 2) || (px.size() == 3) );
		st << "RG\n";
		st << scale*ar*(px[0] - x0) << ' ' << scale*(py[0] - y0) << " m "
		   << scale*ar*(px[1] - x0) << ' ' << scale*(py[1] - y0) << " l\n";
		st << "S\n";
	    }
	}