Function plot_P2_isoline

plugin/seq/plotPDF.cpp:2700–2751

    } // for triangle

    st << "Q\n";

    return;
}

void plot_P2_isoline( std::stringstream &Content, const Fem2D::Mesh &Th, const KN<double> &f_P2_,
		      const KNM<double> &palette,
		      const int sizex, const int sizey, const double scale, const double ar,
		      const double x0, const double y0, const double y1,
		      const int marginl, const int marginb,
		      const double textfontsize, const bool monochrome,
		      const bool legend, const int prec, const bool logscale,
		      const double withmesh,
		      const int NISOLINES, const KN<double>*const viso,
		      const double linewidth )
{
    const double EPS = 1e-10;

    //------------------------------
    // values to plot
    //------------------------------
    const double fmax = (viso)? viso->max(): f_P2_.max();
    const double fmin = (viso)? viso->min(): f_P2_.min();
    std::vector<double> isoline_val;

    find_isoline_values( isoline_val, fmax, fmin, NISOLINES, viso, logscale );

    // color/gray-scale range to plot
#if 1
    // If user specifies an irrelevant color range (in viso array),
    // we do not take care of it.
    const double cmax = *max_element( isoline_val.begin(), isoline_val.end() );
    const double cmin = *min_element( isoline_val.begin(), isoline_val.end() );
#else
    // map color/grayscale range to function values
    const double cmax = fmax;
    const double cmin = fmin;
#endif

    std::stringstream &st = Content;
    st.str("");

    if( 0 < withmesh )
	overlayMesh( st, Th, withmesh, scale, ar, x0, y0, marginl, marginb );
    
    plot_P2_isoline_body( st, Th, f_P2_, isoline_val, cmin, cmax, palette, scale, ar, x0, y0,
			  marginl, marginb, monochrome, logscale, NISOLINES, linewidth );

    if( legend )
	drawLegend_contour( st, isoline_val, prec, palette, cmin, cmax,
			    monochrome, logscale, sizex, textfontsize, scale,
			    y1, y0, marginl, marginb);

    drawBoundary( st, Th, scale, ar, x0, y0, marginl, marginb );