Select Git revision
GModelIO_CGNS.cpp
Forked from
gmsh / gmsh
Source project has a limited visibility.
Scale.cpp 11.24 KiB
// $Id: Scale.cpp,v 1.64 2006-11-27 22:22:16 geuzaine Exp $
//
// Copyright (C) 1997-2007 C. Geuzaine, J.-F. Remacle
//
// This program is free software; you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation; either version 2 of the License, or
// (at your option) any later version.
//
// This program is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
// GNU General Public License for more details.
//
// You should have received a copy of the GNU General Public License
// along with this program; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307
// USA.
//
// Please report all bugs and problems to <gmsh@geuz.org>.
#include "Gmsh.h"
#include "GmshUI.h"
#include "Numeric.h"
#include "Draw.h"
#include "Context.h"
#include "Views.h"
#include "gl2ps.h"
extern Context_T CTX;
// Even if all computations in these routines are made in window
// coordinates, double precision is used to work at subpixel accuracy
void draw_scale(Post_View * v,
double xmin, double ymin, double width, double height,
double tic, int horizontal)
{
gl_font(CTX.gl_font_enum, CTX.gl_fontsize);
double font_h = gl_height(); // total font height
double font_a = gl_height() - gl_descent(); // height above ref pt
double box;
if(horizontal){
box = width / (v->NbIso ? v->NbIso : 1); // colorscale box width
}
else{
box = height / (v->NbIso ? v->NbIso : 1); // colorscale box height
}
glPointSize(v->PointSize);
gl2psPointSize(v->PointSize * CTX.print.eps_point_size_factor);
glLineWidth(v->LineWidth);
gl2psLineWidth(v->LineWidth * CTX.print.eps_line_width_factor);
double ValMin, ValMax;
if(v->ExternalViewIndex >= 0){
ValMin = v->ExternalMin;
ValMax = v->ExternalMax;
}
else{
switch(v->RangeType){
case DRAW_POST_RANGE_CUSTOM:
ValMin = v->CustomMin;
ValMax = v->CustomMax;
break;
case DRAW_POST_RANGE_PER_STEP:
if(v->TimeStepMin && v->TimeStepMax){
ValMin = v->TimeStepMin[v->TimeStep]; ValMax = v->TimeStepMax[v->TimeStep];
}
else{ // possible of the view is empty
ValMin = v->Min;
ValMax = v->Max;
}
break;
case DRAW_POST_RANGE_DEFAULT:
default:
ValMin = v->Min;
ValMax = v->Max;
break;
}
}
switch (v->ScaleType) {
case DRAW_POST_LINEAR:
v->GIFV = GiveIndexFromValue_Lin;
v->GVFI = GiveValueFromIndex_Lin;
break;
case DRAW_POST_LOGARITHMIC:
v->GIFV = GiveIndexFromValue_Log;
v->GVFI = GiveValueFromIndex_Log;
break;
case DRAW_POST_DOUBLELOGARITHMIC:
v->GIFV = GiveIndexFromValue_DoubleLog;
v->GVFI = GiveValueFromIndex_DoubleLog;
break;
}
// colorscale
for(int i = 0; i < v->NbIso; i++) {
if(v->IntervalsType == DRAW_POST_DISCRETE) {
PaletteDiscrete(v, v->NbIso, i);
glBegin(GL_QUADS);
if(horizontal){
glVertex2d(xmin + i * box, ymin);
glVertex2d(xmin + (i + 1) * box, ymin);
glVertex2d(xmin + (i + 1) * box, ymin + height);
glVertex2d(xmin + i * box, ymin + height);
}
else{
glVertex2d(xmin, ymin + i * box);
glVertex2d(xmin + width, ymin + i * box);
glVertex2d(xmin + width, ymin + (i + 1) * box);
glVertex2d(xmin, ymin + (i + 1) * box);
}
glEnd();
}
else if(v->IntervalsType == DRAW_POST_CONTINUOUS) {
glBegin(GL_QUADS);
PaletteContinuousLinear(v, ValMin, ValMax,
ValMin + i * (ValMax - ValMin) /
(v->NbIso ? v->NbIso : 1));
if(horizontal){
glVertex2d(xmin + i * box, ymin + height);
glVertex2d(xmin + i * box, ymin);
}
else{
glVertex2d(xmin, ymin + i * box);
glVertex2d(xmin + width, ymin + i * box);
}
PaletteContinuousLinear(v, ValMin, ValMax,
ValMin + (i + 1) * (ValMax - ValMin) /
(v->NbIso ? v->NbIso : 1));
if(horizontal){
glVertex2d(xmin + (i + 1) * box, ymin);
glVertex2d(xmin + (i + 1) * box, ymin + height);
} else{
glVertex2d(xmin + width, ymin + (i + 1) * box);
glVertex2d(xmin, ymin + (i + 1) * box);
}
glEnd();
}
else {
PaletteDiscrete(v, v->NbIso, i);
glBegin(GL_LINES);
if(horizontal){
glVertex2d(xmin + i * box + 0.5 * box, ymin);
glVertex2d(xmin + i * box + 0.5 * box, ymin + height);
}
else{
glVertex2d(xmin, ymin + i * box + 0.5 * box);
glVertex2d(xmin + width, ymin + i * box + 0.5 * box);
}
glEnd();
}
}
// valuescale
char label[1024];
int nbv;
double cv_box;
if(horizontal){
sprintf(label, v->Format, -M_PI/1.e4);
double estim = gl_width(label);
nbv = (width/estim > v->NbIso) ? v->NbIso : ((width/estim > 2.) ? -2 : -1);
cv_box = width / nbv;
}
else{
nbv = (height/font_h > v->NbIso) ? v->NbIso : ((height/font_h > 3.) ? -2 : -1);
cv_box = height / nbv;
}
glColor4ubv((GLubyte *) & CTX.color.text);
if(nbv < 0) { // only min and max (+ possibly mid value) if not enough room
if(v->IntervalsType == DRAW_POST_DISCRETE ||
v->IntervalsType == DRAW_POST_CONTINUOUS) {
sprintf(label, v->Format, ValMin);
if(horizontal){
glRasterPos2d(xmin, ymin + height + tic);
Draw_String_Center(label);
}
else{
glRasterPos2d(xmin + width + tic, ymin - font_a / 3.);
Draw_String(label);
}
sprintf(label, v->Format, ValMax);
if(horizontal){
glRasterPos2d(xmin + width, ymin + height + tic);
Draw_String_Center(label);
}
else{
glRasterPos2d(xmin + width + tic, ymin + height - font_a / 3.);
Draw_String(label);
}
}
else {
sprintf(label, v->Format, ValMin);
if(horizontal){
glRasterPos2d(xmin + (box / 2), ymin + height + tic);
Draw_String_Center(label);
}
else{
glRasterPos2d(xmin + width + tic, ymin + (box / 2) - font_a / 3.);
Draw_String(label); }
sprintf(label, v->Format, ValMax);
if(horizontal){
glRasterPos2d(xmin + width - (box / 2), ymin + height + tic);
Draw_String_Center(label);
}
else{
glRasterPos2d(xmin + width + tic, ymin + height - (box / 2) - font_a / 3.);
Draw_String(label);
}
}
if(nbv == -2){
double Val = v->GVFI(ValMin, ValMax, v->NbIso, v->NbIso/2);
sprintf(label, v->Format, Val);
if(horizontal){
glRasterPos2d(xmin + width/2., ymin + height + tic);
Draw_String_Center(label);
}
else{
glRasterPos2d(xmin + width + tic, ymin + height/2 - font_a / 3.);
Draw_String(label);
}
}
}
else {
if(v->IntervalsType == DRAW_POST_DISCRETE ||
v->IntervalsType == DRAW_POST_CONTINUOUS) {
for(int i = 0; i < nbv + 1; i++) {
double Val = v->GVFI(ValMin, ValMax, nbv + 1, i);
sprintf(label, v->Format, Val);
if(horizontal){
glRasterPos2d(xmin + i * cv_box, ymin + height + tic);
Draw_String_Center(label);
}
else{
glRasterPos2d(xmin + width + tic, ymin + i * cv_box - font_a / 3.);
Draw_String(label);
}
}
}
else {
for(int i = 0; i < nbv; i++) {
double Val = v->GVFI(ValMin, ValMax, nbv, i);
sprintf(label, v->Format, Val);
if(horizontal){
glRasterPos2d(xmin + (2 * i + 1) * (cv_box / 2), ymin + height + tic);
Draw_String_Center(label);
}
else{
glRasterPos2d(xmin + width + tic, ymin + (2 * i + 1) * (cv_box / 2) - font_a / 3.);
Draw_String(label);
}
}
}
}
// the label
if((v->ShowTime == 1 && List_Nbr(v->Time) > 1) ||
(v->ShowTime == 2 && List_Nbr(v->Time) > 0)){
char tmp[256];
sprintf(tmp, v->Format, *(double *)List_Pointer(v->Time, v->TimeStep));
sprintf(label, "%s (%s)", v->Name, tmp);
}
else if((v->ShowTime == 3 && List_Nbr(v->Time) > 1) ||
(v->ShowTime == 4 && List_Nbr(v->Time) > 0)){
sprintf(label, "%s (%d)", v->Name, v->TimeStep);
}
else
sprintf(label, "%s", v->Name);
if(horizontal){
glRasterPos2d(xmin + width / 2., ymin + height + tic + 1.4*font_h);
Draw_String_Center(label);
}
else{
glRasterPos2d(xmin, ymin - 2 * font_h);
Draw_String(label);
}
}
void Draw_Scales(void)
{
static List_T *todraw = NULL;
if(!CTX.post.list)
return;
// scales to draw?
if(!todraw)
todraw = List_Create(5, 5, sizeof(Post_View *));
else
List_Reset(todraw);
for(int i = 0; i < List_Nbr(CTX.post.list); i++) {
Post_View *v = *(Post_View **) List_Pointer(CTX.post.list, i);
if(v->Visible && !v->Dirty && v->ShowScale &&
v->Type == DRAW_POST_3D && !v->TextOnly)
List_Add(todraw, &v);
}
if(!List_Nbr(todraw)) {
return;
}
const double tic = 10., bar_size = 16.;
double width = 0., width_prev = 0., width_total = 0.;
gl_font(CTX.gl_font_enum, CTX.gl_fontsize);
char label[1024];
double largest_number = 0.;
for(int i = 0; i < List_Nbr(todraw); i++) {
Post_View *v = *(Post_View **) List_Pointer(todraw, i);
sprintf(label, v->Format, -M_PI/1.e4);
if(largest_number < gl_width(label))
largest_number = gl_width(label);
}
for(int i = 0; i < List_Nbr(todraw); i++) {
Post_View *v = *(Post_View **) List_Pointer(todraw, i);
if(!v->AutoPosition) {
double x = v->Position[0], y = v->Position[1];
int center = Fix2DCoordinates(&x, &y);
draw_scale(v,
x - ((center & 1) ? v->Size[0]/2. : 0.),
y - v->Size[1] + ((center & 2) ? v->Size[1]/2. : 0.),
v->Size[0], v->Size[1], tic, CTX.post.horizontal_scales);
}
else{
if(CTX.post.horizontal_scales){
double ysep = 20.;
if(List_Nbr(todraw) == 1){
double ww = (CTX.viewport[2] - CTX.viewport[0]) / 2.;
draw_scale(v,
(CTX.viewport[2] - CTX.viewport[0]) / 2. - ww / 2,
CTX.viewport[1] + ysep,
ww, bar_size,
tic, 1); }
else{
double xsep = largest_number/4. + (CTX.viewport[2] - CTX.viewport[0])/10.;
double ww = (CTX.viewport[2] - CTX.viewport[0] - 4 * xsep) / 2.;
if(ww < 20) ww = 20;
draw_scale(v,
(CTX.viewport[2] - CTX.viewport[0])/2. - (i%2 ? -xsep/1.5 : ww+xsep/1.5),
CTX.viewport[1] + ysep + (i/2) * (bar_size+tic+2*gl_height()+ysep),
ww, bar_size,
tic, 1);
}
}
else{
double xsep = 20.;
double dy = 2. * gl_height();
if(List_Nbr(todraw) == 1){
double ysep = (CTX.viewport[3] - CTX.viewport[1]) / 6.;
double hh = CTX.viewport[3] - CTX.viewport[1] - 2 * ysep - dy;
draw_scale(v,
CTX.viewport[0] + xsep,
CTX.viewport[1] + ysep + dy,
bar_size, hh,
tic, 0);
}
else{
double ysep = (CTX.viewport[3] - CTX.viewport[1]) / 15.;
double hh = (CTX.viewport[3] - CTX.viewport[1] - 3 * ysep - 2.5 * dy) / 2.;
draw_scale(v,
CTX.viewport[0] + xsep + width_total + (i / 2) * xsep,
CTX.viewport[1] + ysep + dy + (1 - i % 2) * (hh + 1.5 * dy + ysep),
bar_size, hh,
tic, 0);
}
// compute width
width_prev = width;
sprintf(label, v->Format, -M_PI/1.e4);
width = bar_size + tic + gl_width(label);
if(List_Nbr(v->Time) > 1 && v->ShowTime){
char tmp[256];
sprintf(tmp, v->Format, *(double *)List_Pointer(v->Time, v->TimeStep));
sprintf(label, "%s (%s)", v->Name, tmp);
}
else
sprintf(label, "%s", v->Name);
if(gl_width(label) > width)
width = gl_width(label);
if(i % 2)
width_total += DMAX(bar_size + width, bar_size + width_prev);
}
}
}
}