REPOS_ERICCA/mesh/imprime.cpp

//---------------------------------------------------------------------------
#include <vcl\vcl.h>
#include <vcl\printers.hpp>
#include <math.h>
#include <string.h>
#include <stdlib.h>
#pragma hdrstop

#include "fenetre.h"
#include "s_valeur.h"
//---------------------------------------------------------------------------

 extern void trirea(float *a,int *la,int *ir);


void __fastcall TFenp::imprime(void)
{
int i,xc,yc;
float zoomx,zoomy;
int xcmax,ycmin,ycmax,xcmin;
int Height,Height2,Width;
char mess[500];
float zoomtmp,xcentretmp,ycentretmp;

Width=(int)(Printer()->PageWidth*rapport);
Height=(int)(Printer()->PageHeight*rapport*0.95);
Height2=(int)(Printer()->PageHeight*rapport);
/* normalisation sur la page */
for (i=0;i<nb_point;i++)
        {
    projection(point[i].x,point[i].y,point[i].z,&xc,&yc,0);
    if (i==0)
        {
        xcmax=xc;
        xcmin=xc;
        ycmax=yc;
        ycmin=yc;
        }
    if (xc>xcmax) xcmax=xc;
    if (yc>ycmax) ycmax=yc;
    if (xc<xcmin) xcmin=xc;
    if (yc<ycmin) ycmin=yc;
    }
xcmax++;xcmin--;
ycmax++;ycmin--;
if (xcmax!=xcmin) zoomx=Width*1./(xcmax-xcmin);  else zoomx=1E30;
if (ycmax!=ycmin) zoomy=Height*1./(ycmax-ycmin); else zoomy=1E30;
if (zoomx>zoomy) zoom_p=zoomy; else zoom_p=zoomx;
zoom_p=zoom_p/1.15;
xcentre_p=(Width/2.-zoom_p*(xcmin+xcmax)/2.);
ycentre_p=(Height/2.-zoom_p*(ycmin+ycmax)/2.);
zoomtmp=zoom;
xcentretmp=xcentre;
ycentretmp=ycentre;
norme();
zoom_p=zoom_p*zoomtmp/zoom;
zoom=zoomtmp;
xcentre_p=xcentre_p+(xcentretmp-xcentre)/Fenp->Width*Width;
ycentre_p=ycentre_p+(ycentretmp-ycentre)/Fenp->Height*Height;
xcentre=xcentretmp;
ycentre=ycentretmp;


if (Geometrie1->Checked==true) strcpy(mess,"Geometrie de ");
if (Maillage1D1->Checked==true) strcpy(mess,"Maillage 1D de ");
if (Maillage2D1->Checked==true) strcpy(mess,"Maillage 2D de ");
if (Maillage3D1->Checked==true) strcpy(mess,"Maillage 3D de ");
if (Remaillage2D1->Checked==true) strcpy(mess,"Remaillage 2D de ");
if (Remaillage3D1->Checked==true) strcpy(mess,"Remaillage 3D de ");
strcat(mess,Fenp->Fichier);
strcpy(info,"Entrez la l�gende du dessin ?");
Valeur->Edit1->Text=mess;
Valeur->ShowModal();
strcpy(mess,Valeur->Edit1->Text.c_str());


/* debut d 'impression */
Printer()->Title=Fenp->Fichier;
Printer()->BeginDoc();
Printer()->Canvas->Pen->Color=clBlack;
Printer()->Canvas->MoveTo(1,1);
Printer()->Canvas->LineTo(1,Height2-2);
Printer()->Canvas->LineTo(Width-2,Height2-2);
Printer()->Canvas->LineTo(Width-2,1);
Printer()->Canvas->LineTo(1,1);
Printer()->Canvas->MoveTo(1,(int)(Height));
Printer()->Canvas->LineTo(Width-2,(int)(Height));

if (Geometrie1->Checked==true) print_geo();
if (Maillage1D1->Checked==true) print_1D();
if (Maillage2D1->Checked==true) print_2D();
if (Maillage3D1->Checked==true) print_3D();
if (Remaillage2D1->Checked==true) print_2D();
if (Remaillage3D1->Checked==true) print_3D();
if (Front3D1->Checked==true) print_2D();
print_axe();
Printer()->Canvas->Font->Color=clBlack;
Printer()->Canvas->Font->Size=(int)(16*rapport);
Printer()->Canvas->Brush->Color=clWhite;
Printer()->Canvas->Pen->Color=clBlack;
Printer()->Canvas->Rectangle(1,(int)(Height),Width-2,Height2-2);
Printer()->Canvas->TextOut((int)(0.05*Width),(int)(Height+0.3*(Height2-Height)),mess);
Printer()->EndDoc();
}


void __fastcall TFenp::print_geo(void)
{
int i,n1,n2;
float x1,x2,y1,y2,z1,z2;
int xc1,xc2,yc1,yc2;

Printer()->Canvas->Pen->Color=clBlack;
for (i=0;i<nb_segment;i++)
        {
        n1=segment[i].n1;
        n2=segment[i].n2;
        x1=point[n1].x;
        y1=point[n1].y;
        z1=point[n1].z;
        x2=point[n2].x;
        y2=point[n2].y;
        z2=point[n2].z;
        projection(x1,y1,z1,&xc1,&yc1,2);
        projection(x2,y2,z2,&xc2,&yc2,2);
    Printer()->Canvas->MoveTo(xc1,yc1);
        Printer()->Canvas->LineTo(xc2,yc2);
        }
}


void __fastcall TFenp::print_1D(void)
{
int i,n1,n2;
float x1,x2,y1,y2,z1,z2;
int xc1,xc2,yc1,yc2;
char mess[20];

Printer()->Canvas->Pen->Color=clRed;
Printer()->Canvas->Font->Color=clFuchsia;
for (i=0;i<nb_point;i++)
        {
    x1=point[i].x;
        y1=point[i].y;
        z1=point[i].z;
        projection(x1,y1,z1,&xc1,&yc1,2);
    Printer()->Canvas->MoveTo(xc1-2,yc1);
        Printer()->Canvas->LineTo(xc1+2,yc1);
        Printer()->Canvas->MoveTo(xc1,yc1-2);
        Printer()->Canvas->LineTo(xc1,yc1+2);
    if (Noeud1->Checked==true)
        {
        sprintf(mess,"%d",i);
        Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
        }
        }
Printer()->Canvas->Font->Color=clYellow;
for (i=0;i<nb_segment;i++)
        {
        if ((segment[i].select==0)&& (Fenp->Debug1->Checked)) continue;
        n1=segment[i].n1;
        n2=segment[i].n2;
        x1=point[n1].x;
        y1=point[n1].y;
        z1=point[n1].z;
        x2=point[n2].x;
        y2=point[n2].y;
        z2=point[n2].z;
        projection(x1,y1,z1,&xc1,&yc1,2);
        projection(x2,y2,z2,&xc2,&yc2,2);
    Printer()->Canvas->MoveTo(xc1,yc1);
        Printer()->Canvas->LineTo(xc2,yc2);
        if (Maille2->Checked==true)
        {
        sprintf(mess,"%d",i);
        Printer()->Canvas->TextOut((int)((xc1+xc2)*0.5),(int)((yc1+yc2)*0.5),mess);
        }
        }
if (numero_noeud>=0)
        {
    Printer()->Canvas->Font->Color=clBlack;
        x1=point[numero_noeud].x;
        y1=point[numero_noeud].y;
        z1=point[numero_noeud].z;
        projection(x1,y1,z1,&xc1,&yc1,2);
    sprintf(mess,"%d",numero_noeud);
    Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
        }
if (numero_maille>=0)
        {
    Printer()->Canvas->Font->Color=clBlack;
        n1=segment[numero_maille].n1;
        n2=segment[numero_maille].n2;
        x1=point[n1].x;
        y1=point[n1].y;
        z1=point[n1].z;
        x2=point[n2].x;
        y2=point[n2].y;
        z2=point[n2].z;
        projection(x1,y1,z1,&xc1,&yc1,2);
        projection(x2,y2,z2,&xc2,&yc2,2);
    sprintf(mess,"%d",numero_maille);
    Printer()->Canvas->TextOut((int)((xc1+xc2)*0.5),(int)((yc1+yc2)*0.5),mess);
    }
}
//---------------------------------------------------------------------------
void __fastcall TFenp::print_2D(void)
{
int i,n1,n2,n3;
float x1,x2,y1,y2,z1,z2,x3,y3,z3;
float ab[3],ac[3],g[3],crit;
int xc1,xc2,yc1,yc2,xc3,yc3;
TRect r;
char mess[20];

if (fc==1)
        {
    print_reel();
    return;
    }
Printer()->Canvas->Pen->Color=clRed;
if (Eclatee1->Checked==true) Printer()->Canvas->Pen->Color=clOlive;
for (i=0;i<nb_triangle;i++)
        {
        if ((triangle[i].num_ent==numface) || (numface==(-1)))
        {
        if ((triangle[i].select==0) && (Fenp->Debug1->Checked)) continue;
        n1=triangle[i].n1;
                n2=triangle[i].n2;
        n3=triangle[i].n3;
        x1=point[n1].x;
                y1=point[n1].y;
                z1=point[n1].z;
                x2=point[n2].x;
                y2=point[n2].y;
                z2=point[n2].z;
                x3=point[n3].x;
                y3=point[n3].y;
                z3=point[n3].z;
        if (Normal1->Checked==false)
                {
            if (Qualite1->Checked==true)
                {
                crit=crit_cal2d(n1,n2,n3);
                    if (crit>0.5) Printer()->Canvas->Pen->Color=clLime;
                    else if (crit>0.2) Printer()->Canvas->Pen->Color=clBlue;
                    else if (crit>0.1) Printer()->Canvas->Pen->Color=clRed;
                    else Printer()->Canvas->Pen->Color=clMaroon;
                    }
                        ab[0]=x2-x1;ab[1]=y2-y1;ab[2]=z2-z1;
                        ac[0]=x3-x1;ac[1]=y3-y1;ac[2]=z3-z1;
            g[0]=x1+0.33333*(ab[0]+ac[0]);
                        g[1]=y1+0.33333*(ab[1]+ac[1]);
                        g[2]=z1+0.33333*(ab[2]+ac[2]);
            x1=g[0]+0.8*(x1-g[0]);y1=g[1]+0.8*(y1-g[1]);z1=g[2]+0.8*(z1-g[2]);
                        x2=g[0]+0.8*(x2-g[0]);y2=g[1]+0.8*(y2-g[1]);z2=g[2]+0.8*(z2-g[2]);
                        x3=g[0]+0.8*(x3-g[0]);y3=g[1]+0.8*(y3-g[1]);z3=g[2]+0.8*(z3-g[2]);
                        }
                projection(x1,y1,z1,&xc1,&yc1,2);
                projection(x2,y2,z2,&xc2,&yc2,2);
            projection(x3,y3,z3,&xc3,&yc3,2);
            Printer()->Canvas->MoveTo(xc1,yc1);
                Printer()->Canvas->LineTo(xc2,yc2);
                Printer()->Canvas->LineTo(xc3,yc3);
                Printer()->Canvas->LineTo(xc1,yc1);
        Printer()->Canvas->Font->Color=clFuchsia;
            if (Noeud1->Checked==true)
                {
                sprintf(mess,"%d",n1);
                projection(point[n1].x,point[n1].y,point[n1].z,&xc1,&yc1,2);
                        Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
                sprintf(mess,"%d",n2);
                projection(point[n2].x,point[n2].y,point[n2].z,&xc2,&yc2,2);
                        Printer()->Canvas->TextOut(xc2+2,yc2+2,mess);
                sprintf(mess,"%d",n3);
                projection(point[n3].x,point[n3].y,point[n3].z,&xc3,&yc3,2);
                        Printer()->Canvas->TextOut(xc3+2,yc3+2,mess);
                }
        Printer()->Canvas->Font->Color=clYellow;
                if (Maille2->Checked==true)
                        {
                sprintf(mess,"%d",i);
                Printer()->Canvas->TextOut((int)((xc1+xc2+xc3)*0.333333),(int)((yc1+yc2+yc3)*0.33333333),mess);
                }
                }
    }
if (numero_noeud>=0)
        {
    Printer()->Canvas->Font->Color=clBlack;
        x1=point[numero_noeud].x;
        y1=point[numero_noeud].y;
        z1=point[numero_noeud].z;
        projection(x1,y1,z1,&xc1,&yc1,2);
    sprintf(mess,"%d",numero_noeud);
    Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
        }
if (numero_maille>=0)
        {
    Printer()->Canvas->Font->Color=clBlack;
        n1=triangle[numero_maille].n1;
    n2=triangle[numero_maille].n2;
        n3=triangle[numero_maille].n3;
        x1=point[n1].x;
        y1=point[n1].y;
        z1=point[n1].z;
        x2=point[n2].x;
        y2=point[n2].y;
        z2=point[n2].z;
        x3=point[n3].x;
        y3=point[n3].y;
        z3=point[n3].z;
    projection(x1,y1,z1,&xc1,&yc1,2);
        projection(x2,y2,z2,&xc2,&yc2,2);
    projection(x3,y3,z3,&xc3,&yc3,2);
    sprintf(mess,"%d",numero_maille);
    Printer()->Canvas->TextOut((int)((xc1+xc2+xc3)*0.333333),(int)((yc1+yc2+yc3)*0.33333333),mess);
    }

}
//---------------------------------------------------------------------------
void __fastcall TFenp::print_3D(void)
{
int i,n1,n2,n3,n4;
float x1,x2,y1,y2,z1,z2,x3,y3,z3,x4,y4,z4;
float ab[3],ac[3],ad[3],g[3],crit;
int xc1,xc2,yc1,yc2,xc3,yc3,xc4,yc4;
char mess[20];
TRect r;

if ((Activerlerendu3D1->Checked==true)&&(fc==1))
        {
    print_reel();
    return;
    }
Printer()->Canvas->Pen->Color=clRed;
if (Eclatee1->Checked==true) Printer()->Canvas->Pen->Color=clOlive;
for (i=0;i<nb_tetra;i++)
        {
        if ((tetra[i].select==0)&& (Fenp->Debug1->Checked)) continue;
        n1=tetra[i].n1;
        n2=tetra[i].n2;
    n3=tetra[i].n3;
    n4=tetra[i].n4;
    x1=point[n1].x;
        y1=point[n1].y;
        z1=point[n1].z;
        x2=point[n2].x;
        y2=point[n2].y;
        z2=point[n2].z;
        x3=point[n3].x;
        y3=point[n3].y;
        z3=point[n3].z;
    x4=point[n4].x;
    y4=point[n4].y;
    z4=point[n4].z;
    if (Normal1->Checked==false)
        {
        if (Qualite1->Checked==true)
                {
                crit=crit_cal3d(n1,n2,n3,n4);
            if (crit>0.5) Printer()->Canvas->Pen->Color=clLime;
            else if (crit>0.2) Printer()->Canvas->Pen->Color=clBlue;
            else if (crit>0.1) Printer()->Canvas->Pen->Color=clRed;
            else Printer()->Canvas->Pen->Color=clMaroon;
            }
                ab[0]=x2-x1;ab[1]=y2-y1;ab[2]=z2-z1;
                ac[0]=x3-x1;ac[1]=y3-y1;ac[2]=z3-z1;
        ad[0]=x4-x1;ad[1]=y4-y1;ad[2]=z4-z1;
        g[0]=x1+0.25*(ab[0]+ac[0]+ad[0]);
                g[1]=y1+0.25*(ab[1]+ac[1]+ad[1]);
                g[2]=z1+0.25*(ab[2]+ac[2]+ad[2]);
        x1=g[0]+0.8*(x1-g[0]);y1=g[1]+0.8*(y1-g[1]);z1=g[2]+0.8*(z1-g[2]);
                x2=g[0]+0.8*(x2-g[0]);y2=g[1]+0.8*(y2-g[1]);z2=g[2]+0.8*(z2-g[2]);
                x3=g[0]+0.8*(x3-g[0]);y3=g[1]+0.8*(y3-g[1]);z3=g[2]+0.8*(z3-g[2]);
                x4=g[0]+0.8*(x4-g[0]);y4=g[1]+0.8*(y4-g[1]);z4=g[2]+0.8*(z4-g[2]);
                }
        projection(x1,y1,z1,&xc1,&yc1,2);
        projection(x2,y2,z2,&xc2,&yc2,2);
    projection(x3,y3,z3,&xc3,&yc3,2);
    projection(x4,y4,z4,&xc4,&yc4,2);
    Printer()->Canvas->MoveTo(xc1,yc1);
        Printer()->Canvas->LineTo(xc2,yc2);
        Printer()->Canvas->LineTo(xc3,yc3);
        Printer()->Canvas->LineTo(xc1,yc1);
        Printer()->Canvas->MoveTo(xc4,yc4);
        Printer()->Canvas->LineTo(xc1,yc1);
        Printer()->Canvas->MoveTo(xc4,yc4);
        Printer()->Canvas->LineTo(xc2,yc2);
        Printer()->Canvas->MoveTo(xc4,yc4);
        Printer()->Canvas->LineTo(xc3,yc3);
        Printer()->Canvas->Font->Color=clFuchsia;
        if (Noeud1->Checked==true)
        {
        sprintf(mess,"%d",n1);
        Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
        sprintf(mess,"%d",n2);
        Printer()->Canvas->TextOut(xc2+2,yc2+2,mess);
        sprintf(mess,"%d",n3);
        Printer()->Canvas->TextOut(xc3+2,yc3+2,mess);
        sprintf(mess,"%d",n4);
        Printer()->Canvas->TextOut(xc4+2,yc4+2,mess);
        }
    Printer()->Canvas->Font->Color=clYellow;
        if (Maille2->Checked==true)
        {
        sprintf(mess,"%d",i);
        Printer()->Canvas->TextOut((int)((xc1+xc2+xc3+xc4)*0.25),(int)((yc1+yc2+yc3+yc4)*0.25),mess);
        }
        }
if (numero_noeud>=0)
        {
    Printer()->Canvas->Font->Color=clBlack;
        x1=point[numero_noeud].x;
        y1=point[numero_noeud].y;
        z1=point[numero_noeud].z;
        projection(x1,y1,z1,&xc1,&yc1,2);
    sprintf(mess,"%d",numero_noeud);
    Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
        }
if (numero_maille>=0)
        {
    Printer()->Canvas->Font->Color=clBlack;
        n1=tetra[numero_maille].n1;
        n2=tetra[numero_maille].n2;
    n3=tetra[numero_maille].n3;
    n4=tetra[numero_maille].n4;
    x1=point[n1].x;
        y1=point[n1].y;
        z1=point[n1].z;
        x2=point[n2].x;
        y2=point[n2].y;
        z2=point[n2].z;
        x3=point[n3].x;
        y3=point[n3].y;
        z3=point[n3].z;
    x4=point[n4].x;
    y4=point[n4].y;
    z4=point[n4].z;
    projection(x1,y1,z1,&xc1,&yc1,2);
        projection(x2,y2,z2,&xc2,&yc2,2);
    projection(x3,y3,z3,&xc3,&yc3,2);
    projection(x4,y4,z4,&xc4,&yc4,2);
    sprintf(mess,"%d",numero_maille);
    Printer()->Canvas->TextOut((int)((xc1+xc2+xc3+xc4)*0.25),(int)((yc1+yc2+yc3+yc4)*0.25),mess);
    }


}
//---------------------------------------------------------------------------
void __fastcall TFenp::print_reel(void)
{
int *iordre;
float *distance;
int xpoint[8];
TPoint ypoint[8];
float xe,ye,ze,xee,yee,zee;
float xa,ya,za;
float xb,yb,zb;
float xc,yc,zc;
float xg,yg,zg;
float nx,ny,nz;
char mess[20];
struct s_triangle *tete=NULL,*tri;
int i;
int color;
float dismin,dismax,diff,g[4];
static float arot[3][3]={{-1.,-1.,-1.},{-1.,-1.,-1.},{-1.,-1.,-1.}};
int pareil;


pareil=0;
if (tete!=NULL)
if (arot[0][0]==rot[0][0])
if (arot[1][0]==rot[1][0])
if (arot[2][0]==rot[2][0])
if (arot[0][1]==rot[0][1])
if (arot[1][1]==rot[1][1])
if (arot[2][1]==rot[2][1])
if (arot[0][2]==rot[0][2])
if (arot[1][2]==rot[1][2])
if (arot[2][2]==rot[2][2]) pareil=1;


if (pareil==0)
{
        tete=NULL;
        arot[0][0]=rot[0][0];
        arot[1][0]=rot[1][0];
        arot[2][0]=rot[2][0];
        arot[0][1]=rot[0][1];
        arot[1][1]=rot[1][1];
        arot[2][1]=rot[2][1];
        arot[0][2]=rot[0][2];
        arot[1][2]=rot[1][2];
        arot[2][2]=rot[2][2];
        /* oeuil */
        xe=1000.*maxcoord-0.5*(xmin+xmax);
        ye=1000.*maxcoord-0.5*(ymin+ymax);
        ze=1000.*maxcoord-0.5*(zmin+zmax);
        /* oeuil dans le repere local */
        xee=rot[0][0]*xe+rot[1][0]*ye+rot[2][0]*ze;
        yee=rot[0][1]*xe+rot[1][1]*ye+rot[2][1]*ze;
        zee=rot[0][2]*xe+rot[1][2]*ye+rot[2][2]*ze;
        iordre=(int *)calloc(nb_triangle+1,sizeof(int));
        distance=(float *)calloc(nb_triangle+1,sizeof(float));
        for (i=0;i<nb_triangle;i++)
                {
        triangle[i].dis=0.;
        triangle[i].suivant=NULL;
        xa=point[triangle[i].n1].x;
        ya=point[triangle[i].n1].y;
        za=point[triangle[i].n1].z;
        xb=point[triangle[i].n2].x;
        yb=point[triangle[i].n2].y;
        zb=point[triangle[i].n2].z;
        xc=point[triangle[i].n3].x;
        yc=point[triangle[i].n3].y;
        zc=point[triangle[i].n3].z;
        xg=0.333333*(xa+xb+xc);
        yg=0.333333*(ya+yb+yc);
        zg=0.333333*(za+zb+zc);
    triangle[i].dis=(xee-xg)*(xee-xg)+(yee-yg)*(yee-yg)+(zee-zg)*(zee-zg);
    if (i==0)
                {
            dismin=triangle[i].dis;
                        dismax=triangle[i].dis;
                        }
        else
                {
            if (triangle[i].dis<dismin) dismin=triangle[i].dis;
                        if (triangle[i].dis>dismax) dismax=triangle[i].dis;
            }
    nx=(yb-ya)*(zc-za)-(zb-za)*(yc-ya);
    ny=(xc-xa)*(zb-za)-(xb-xa)*(zc-za);
    nz=(xb-xa)*(yc-ya)-(yb-ya)*(xc-xa);
    triangle[i].angle=(float)((nx*(xee-xg)+ny*(yee-yg)+nz*(zee-zg))/sqrt((double)((nx*nx+ny*ny+nz*nz)*(triangle[i].dis))));
    if (triangle[i].angle>=1.) triangle[i].angle=0.99999999;
    if (triangle[i].angle<=(-1.)) triangle[i].angle=(-0.99999999);
        distance[i]=triangle[i].dis;
        iordre[i]=i+1;

    }
trirea(distance,&(nb_triangle),iordre);
diff=dismax-dismin;
dismin=dismin-0.025*diff;
dismax=dismax+0.025*diff;
Printer()->Canvas->Pen->Color=Cshade[0];
Printer()->Canvas->Brush->Color=Cshade[1];
for (i=nb_triangle-1;i>(-1);i--)
        {
    tri=&(triangle[iordre[i]-1]);
    if ((tri->select==0)&& (Fenp->Debug1->Checked)) continue;
        if (((tri->num_ent==numface) || (numface==(-1))&&(tri->dis-dismin*loin-(1.-loin)*dismax>0.)))
                {
        xa=point[tri->n1].x;
                ya=point[tri->n1].y;
                za=point[tri->n1].z;
                xb=point[tri->n2].x;
                yb=point[tri->n2].y;
                zb=point[tri->n2].z;
                xc=point[tri->n3].x;
                yc=point[tri->n3].y;
                zc=point[tri->n3].z;
            if (Normal1->Checked==false)
                {
            g[0]=0.33333*(xa+xb+xc);
                        g[1]=0.33333*(ya+yb+yc);
                        g[2]=0.33333*(za+zb+zc);
            xa=g[0]+homo*(xa-g[0]);ya=g[1]+homo*(ya-g[1]);za=g[2]+homo*(za-g[2]);
                        xb=g[0]+homo*(xb-g[0]);yb=g[1]+homo*(yb-g[1]);zb=g[2]+homo*(zb-g[2]);
                        xc=g[0]+homo*(xc-g[0]);yc=g[1]+homo*(yc-g[1]);zc=g[2]+homo*(zc-g[2]);
            }
        projection(xa,ya,za,&xpoint[0],&xpoint[1],2);
                projection(xb,yb,zb,&xpoint[2],&xpoint[3],2);
                projection(xc,yc,zc,&xpoint[4],&xpoint[5],2);
        ypoint[0].x=xpoint[0];
        ypoint[0].y=xpoint[1];
        ypoint[1].x=xpoint[2];
        ypoint[1].y=xpoint[3];
                ypoint[2].x=xpoint[4];
        ypoint[2].y=xpoint[5];
        color=(int)(-10*fabs((double)(tri->angle))+10);
        if (color>10) color=10;
        if (color<1) color=1;
                Printer()->Canvas->Brush->Color=Cshade[color];
                Printer()->Canvas->Polygon(ypoint,2);
        Printer()->Canvas->Font->Color=clFuchsia;
                if (Noeud1->Checked==true)
                {
                sprintf(mess,"%d",tri->n1);
                Printer()->Canvas->TextOut(xpoint[0]+2,xpoint[1]+2,mess);
                sprintf(mess,"%d",tri->n2);
                Printer()->Canvas->TextOut(xpoint[2]+2,xpoint[3]+2,mess);
                sprintf(mess,"%d",tri->n3);
                Printer()->Canvas->TextOut(xpoint[4]+2,xpoint[5]+2,mess);
                }
        if (tri->n1==numero_noeud)
                {
            Printer()->Canvas->Font->Color=clBlack;
                sprintf(mess,"%d",tri->n1);
                Printer()->Canvas->TextOut(xpoint[0]+2,xpoint[1]+2,mess);
                }
        if (tri->n2==numero_noeud)
                {
            Printer()->Canvas->Font->Color=clBlack;
                sprintf(mess,"%d",tri->n2);
                Printer()->Canvas->TextOut(xpoint[2]+2,xpoint[3]+2,mess);
                }
        if (tri->n3==numero_noeud)
                {
            Printer()->Canvas->Font->Color=clBlack;
                sprintf(mess,"%d",tri->n3);
                Printer()->Canvas->TextOut(xpoint[4]+2,xpoint[5]+2,mess);
                }
        if (tri->num!=numero_maille) Printer()->Canvas->Font->Color=clYellow;
        else
                {
                Printer()->Canvas->Font->Color=clBlack;
                sprintf(mess,"%d",tri->num);
                Printer()->Canvas->TextOut((int)((xpoint[0]+xpoint[2]+xpoint[4])*0.333333),(int)((xpoint[1]+xpoint[3]+xpoint[5])*0.33333333),mess);
                }
                if (Maille2->Checked==true)
                {
                sprintf(mess,"%d",tri->num);
                Printer()->Canvas->TextOut((int)((xpoint[0]+xpoint[2]+xpoint[4])*0.333333),(int)((xpoint[1]+xpoint[3]+xpoint[5])*0.33333333),mess);
                }
                }
        }
        free(iordre);
        free(distance);
        }
}
void __fastcall TFenp::print_axe(void)
{
 int xc,yc;

 Printer()->Canvas->Pen->Color=clBlue;
 projection(75.,0.,0.,&xc,&yc,0);
 Printer()->Canvas->MoveTo(100,100);
 Printer()->Canvas->LineTo(100+xc,100+yc);
 Printer()->Canvas->Pen->Color=clRed;
 projection(0.,75.,0.,&xc,&yc,0);
 Printer()->Canvas->MoveTo(100,100);
 Printer()->Canvas->LineTo(100+xc,100+yc);
 Printer()->Canvas->Pen->Color=clBlack;
 projection(0.,0.,75.,&xc,&yc,0);
 Printer()->Canvas->MoveTo(100,100);
 Printer()->Canvas->LineTo(100+xc,100+yc);
 }


Revision:	3
Committed:	Tue Jun 12 12:42:51 2007 UTC (17 years, 11 months ago)
File size:	19997 byte(s)
Log Message:
#	Content
1	//---------------------------------------------------------------------------
2	#include <vcl\vcl.h>
3	#include <vcl\printers.hpp>
4	#include <math.h>
5	#include <string.h>
6	#include <stdlib.h>
7	#pragma hdrstop
8
9	#include "fenetre.h"
10	#include "s_valeur.h"
11	//---------------------------------------------------------------------------
12
13	extern void trirea(float a,int la,int *ir);
14
15
16	void __fastcall TFenp::imprime(void)
17	{
18	int i,xc,yc;
19	float zoomx,zoomy;
20	int xcmax,ycmin,ycmax,xcmin;
21	int Height,Height2,Width;
22	char mess[500];
23	float zoomtmp,xcentretmp,ycentretmp;
24
25	Width=(int)(Printer()->PageWidth*rapport);
26	Height=(int)(Printer()->PageHeightrapport0.95);
27	Height2=(int)(Printer()->PageHeight*rapport);
28	/* normalisation sur la page */
29	for (i=0;i<nb_point;i++)
30	{
31	projection(point[i].x,point[i].y,point[i].z,&xc,&yc,0);
32	if (i==0)
33	{
34	xcmax=xc;
35	xcmin=xc;
36	ycmax=yc;
37	ycmin=yc;
38	}
39	if (xc>xcmax) xcmax=xc;
40	if (yc>ycmax) ycmax=yc;
41	if (xc<xcmin) xcmin=xc;
42	if (yc<ycmin) ycmin=yc;
43	}
44	xcmax++;xcmin--;
45	ycmax++;ycmin--;
46	if (xcmax!=xcmin) zoomx=Width*1./(xcmax-xcmin); else zoomx=1E30;
47	if (ycmax!=ycmin) zoomy=Height*1./(ycmax-ycmin); else zoomy=1E30;
48	if (zoomx>zoomy) zoom_p=zoomy; else zoom_p=zoomx;
49	zoom_p=zoom_p/1.15;
50	xcentre_p=(Width/2.-zoom_p*(xcmin+xcmax)/2.);
51	ycentre_p=(Height/2.-zoom_p*(ycmin+ycmax)/2.);
52	zoomtmp=zoom;
53	xcentretmp=xcentre;
54	ycentretmp=ycentre;
55	norme();
56	zoom_p=zoom_p*zoomtmp/zoom;
57	zoom=zoomtmp;
58	xcentre_p=xcentre_p+(xcentretmp-xcentre)/Fenp->Width*Width;
59	ycentre_p=ycentre_p+(ycentretmp-ycentre)/Fenp->Height*Height;
60	xcentre=xcentretmp;
61	ycentre=ycentretmp;
62
63
64	if (Geometrie1->Checked==true) strcpy(mess,"Geometrie de ");
65	if (Maillage1D1->Checked==true) strcpy(mess,"Maillage 1D de ");
66	if (Maillage2D1->Checked==true) strcpy(mess,"Maillage 2D de ");
67	if (Maillage3D1->Checked==true) strcpy(mess,"Maillage 3D de ");
68	if (Remaillage2D1->Checked==true) strcpy(mess,"Remaillage 2D de ");
69	if (Remaillage3D1->Checked==true) strcpy(mess,"Remaillage 3D de ");
70	strcat(mess,Fenp->Fichier);
71	strcpy(info,"Entrez la l�gende du dessin ?");
72	Valeur->Edit1->Text=mess;
73	Valeur->ShowModal();
74	strcpy(mess,Valeur->Edit1->Text.c_str());
75
76
77
78
79
80
81
82	/* debut d 'impression */
83	Printer()->Title=Fenp->Fichier;
84	Printer()->BeginDoc();
85	Printer()->Canvas->Pen->Color=clBlack;
86	Printer()->Canvas->MoveTo(1,1);
87	Printer()->Canvas->LineTo(1,Height2-2);
88	Printer()->Canvas->LineTo(Width-2,Height2-2);
89	Printer()->Canvas->LineTo(Width-2,1);
90	Printer()->Canvas->LineTo(1,1);
91	Printer()->Canvas->MoveTo(1,(int)(Height));
92	Printer()->Canvas->LineTo(Width-2,(int)(Height));
93
94	if (Geometrie1->Checked==true) print_geo();
95	if (Maillage1D1->Checked==true) print_1D();
96	if (Maillage2D1->Checked==true) print_2D();
97	if (Maillage3D1->Checked==true) print_3D();
98	if (Remaillage2D1->Checked==true) print_2D();
99	if (Remaillage3D1->Checked==true) print_3D();
100	if (Front3D1->Checked==true) print_2D();
101	print_axe();
102	Printer()->Canvas->Font->Color=clBlack;
103	Printer()->Canvas->Font->Size=(int)(16*rapport);
104	Printer()->Canvas->Brush->Color=clWhite;
105	Printer()->Canvas->Pen->Color=clBlack;
106	Printer()->Canvas->Rectangle(1,(int)(Height),Width-2,Height2-2);
107	Printer()->Canvas->TextOut((int)(0.05Width),(int)(Height+0.3(Height2-Height)),mess);
108	Printer()->EndDoc();
109	}
110
111
112
113	void __fastcall TFenp::print_geo(void)
114	{
115	int i,n1,n2;
116	float x1,x2,y1,y2,z1,z2;
117	int xc1,xc2,yc1,yc2;
118
119	Printer()->Canvas->Pen->Color=clBlack;
120	for (i=0;i<nb_segment;i++)
121	{
122	n1=segment[i].n1;
123	n2=segment[i].n2;
124	x1=point[n1].x;
125	y1=point[n1].y;
126	z1=point[n1].z;
127	x2=point[n2].x;
128	y2=point[n2].y;
129	z2=point[n2].z;
130	projection(x1,y1,z1,&xc1,&yc1,2);
131	projection(x2,y2,z2,&xc2,&yc2,2);
132	Printer()->Canvas->MoveTo(xc1,yc1);
133	Printer()->Canvas->LineTo(xc2,yc2);
134	}
135	}
136
137
138	void __fastcall TFenp::print_1D(void)
139	{
140	int i,n1,n2;
141	float x1,x2,y1,y2,z1,z2;
142	int xc1,xc2,yc1,yc2;
143	char mess[20];
144
145	Printer()->Canvas->Pen->Color=clRed;
146	Printer()->Canvas->Font->Color=clFuchsia;
147	for (i=0;i<nb_point;i++)
148	{
149	x1=point[i].x;
150	y1=point[i].y;
151	z1=point[i].z;
152	projection(x1,y1,z1,&xc1,&yc1,2);
153	Printer()->Canvas->MoveTo(xc1-2,yc1);
154	Printer()->Canvas->LineTo(xc1+2,yc1);
155	Printer()->Canvas->MoveTo(xc1,yc1-2);
156	Printer()->Canvas->LineTo(xc1,yc1+2);
157	if (Noeud1->Checked==true)
158	{
159	sprintf(mess,"%d",i);
160	Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
161	}
162	}
163	Printer()->Canvas->Font->Color=clYellow;
164	for (i=0;i<nb_segment;i++)
165	{
166	if ((segment[i].select==0)&& (Fenp->Debug1->Checked)) continue;
167	n1=segment[i].n1;
168	n2=segment[i].n2;
169	x1=point[n1].x;
170	y1=point[n1].y;
171	z1=point[n1].z;
172	x2=point[n2].x;
173	y2=point[n2].y;
174	z2=point[n2].z;
175	projection(x1,y1,z1,&xc1,&yc1,2);
176	projection(x2,y2,z2,&xc2,&yc2,2);
177	Printer()->Canvas->MoveTo(xc1,yc1);
178	Printer()->Canvas->LineTo(xc2,yc2);
179	if (Maille2->Checked==true)
180	{
181	sprintf(mess,"%d",i);
182	Printer()->Canvas->TextOut((int)((xc1+xc2)0.5),(int)((yc1+yc2)0.5),mess);
183	}
184	}
185	if (numero_noeud>=0)
186	{
187	Printer()->Canvas->Font->Color=clBlack;
188	x1=point[numero_noeud].x;
189	y1=point[numero_noeud].y;
190	z1=point[numero_noeud].z;
191	projection(x1,y1,z1,&xc1,&yc1,2);
192	sprintf(mess,"%d",numero_noeud);
193	Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
194	}
195	if (numero_maille>=0)
196	{
197	Printer()->Canvas->Font->Color=clBlack;
198	n1=segment[numero_maille].n1;
199	n2=segment[numero_maille].n2;
200	x1=point[n1].x;
201	y1=point[n1].y;
202	z1=point[n1].z;
203	x2=point[n2].x;
204	y2=point[n2].y;
205	z2=point[n2].z;
206	projection(x1,y1,z1,&xc1,&yc1,2);
207	projection(x2,y2,z2,&xc2,&yc2,2);
208	sprintf(mess,"%d",numero_maille);
209	Printer()->Canvas->TextOut((int)((xc1+xc2)0.5),(int)((yc1+yc2)0.5),mess);
210	}
211	}
212	//---------------------------------------------------------------------------
213	void __fastcall TFenp::print_2D(void)
214	{
215	int i,n1,n2,n3;
216	float x1,x2,y1,y2,z1,z2,x3,y3,z3;
217	float ab[3],ac[3],g[3],crit;
218	int xc1,xc2,yc1,yc2,xc3,yc3;
219	TRect r;
220	char mess[20];
221
222	if (fc==1)
223	{
224	print_reel();
225	return;
226	}
227	Printer()->Canvas->Pen->Color=clRed;
228	if (Eclatee1->Checked==true) Printer()->Canvas->Pen->Color=clOlive;
229	for (i=0;i<nb_triangle;i++)
230	{
231	if ((triangle[i].num_ent==numface) \|\| (numface==(-1)))
232	{
233	if ((triangle[i].select==0) && (Fenp->Debug1->Checked)) continue;
234	n1=triangle[i].n1;
235	n2=triangle[i].n2;
236	n3=triangle[i].n3;
237	x1=point[n1].x;
238	y1=point[n1].y;
239	z1=point[n1].z;
240	x2=point[n2].x;
241	y2=point[n2].y;
242	z2=point[n2].z;
243	x3=point[n3].x;
244	y3=point[n3].y;
245	z3=point[n3].z;
246	if (Normal1->Checked==false)
247	{
248	if (Qualite1->Checked==true)
249	{
250	crit=crit_cal2d(n1,n2,n3);
251	if (crit>0.5) Printer()->Canvas->Pen->Color=clLime;
252	else if (crit>0.2) Printer()->Canvas->Pen->Color=clBlue;
253	else if (crit>0.1) Printer()->Canvas->Pen->Color=clRed;
254	else Printer()->Canvas->Pen->Color=clMaroon;
255	}
256	ab[0]=x2-x1;ab[1]=y2-y1;ab[2]=z2-z1;
257	ac[0]=x3-x1;ac[1]=y3-y1;ac[2]=z3-z1;
258	g[0]=x1+0.33333*(ab[0]+ac[0]);
259	g[1]=y1+0.33333*(ab[1]+ac[1]);
260	g[2]=z1+0.33333*(ab[2]+ac[2]);
261	x1=g[0]+0.8(x1-g[0]);y1=g[1]+0.8(y1-g[1]);z1=g[2]+0.8*(z1-g[2]);
262	x2=g[0]+0.8(x2-g[0]);y2=g[1]+0.8(y2-g[1]);z2=g[2]+0.8*(z2-g[2]);
263	x3=g[0]+0.8(x3-g[0]);y3=g[1]+0.8(y3-g[1]);z3=g[2]+0.8*(z3-g[2]);
264	}
265	projection(x1,y1,z1,&xc1,&yc1,2);
266	projection(x2,y2,z2,&xc2,&yc2,2);
267	projection(x3,y3,z3,&xc3,&yc3,2);
268	Printer()->Canvas->MoveTo(xc1,yc1);
269	Printer()->Canvas->LineTo(xc2,yc2);
270	Printer()->Canvas->LineTo(xc3,yc3);
271	Printer()->Canvas->LineTo(xc1,yc1);
272	Printer()->Canvas->Font->Color=clFuchsia;
273	if (Noeud1->Checked==true)
274	{
275	sprintf(mess,"%d",n1);
276	projection(point[n1].x,point[n1].y,point[n1].z,&xc1,&yc1,2);
277	Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
278	sprintf(mess,"%d",n2);
279	projection(point[n2].x,point[n2].y,point[n2].z,&xc2,&yc2,2);
280	Printer()->Canvas->TextOut(xc2+2,yc2+2,mess);
281	sprintf(mess,"%d",n3);
282	projection(point[n3].x,point[n3].y,point[n3].z,&xc3,&yc3,2);
283	Printer()->Canvas->TextOut(xc3+2,yc3+2,mess);
284	}
285	Printer()->Canvas->Font->Color=clYellow;
286	if (Maille2->Checked==true)
287	{
288	sprintf(mess,"%d",i);
289	Printer()->Canvas->TextOut((int)((xc1+xc2+xc3)0.333333),(int)((yc1+yc2+yc3)0.33333333),mess);
290	}
291	}
292	}
293	if (numero_noeud>=0)
294	{
295	Printer()->Canvas->Font->Color=clBlack;
296	x1=point[numero_noeud].x;
297	y1=point[numero_noeud].y;
298	z1=point[numero_noeud].z;
299	projection(x1,y1,z1,&xc1,&yc1,2);
300	sprintf(mess,"%d",numero_noeud);
301	Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
302	}
303	if (numero_maille>=0)
304	{
305	Printer()->Canvas->Font->Color=clBlack;
306	n1=triangle[numero_maille].n1;
307	n2=triangle[numero_maille].n2;
308	n3=triangle[numero_maille].n3;
309	x1=point[n1].x;
310	y1=point[n1].y;
311	z1=point[n1].z;
312	x2=point[n2].x;
313	y2=point[n2].y;
314	z2=point[n2].z;
315	x3=point[n3].x;
316	y3=point[n3].y;
317	z3=point[n3].z;
318	projection(x1,y1,z1,&xc1,&yc1,2);
319	projection(x2,y2,z2,&xc2,&yc2,2);
320	projection(x3,y3,z3,&xc3,&yc3,2);
321	sprintf(mess,"%d",numero_maille);
322	Printer()->Canvas->TextOut((int)((xc1+xc2+xc3)0.333333),(int)((yc1+yc2+yc3)0.33333333),mess);
323	}
324
325	}
326	//---------------------------------------------------------------------------
327	void __fastcall TFenp::print_3D(void)
328	{
329	int i,n1,n2,n3,n4;
330	float x1,x2,y1,y2,z1,z2,x3,y3,z3,x4,y4,z4;
331	float ab[3],ac[3],ad[3],g[3],crit;
332	int xc1,xc2,yc1,yc2,xc3,yc3,xc4,yc4;
333	char mess[20];
334	TRect r;
335
336	if ((Activerlerendu3D1->Checked==true)&&(fc==1))
337	{
338	print_reel();
339	return;
340	}
341	Printer()->Canvas->Pen->Color=clRed;
342	if (Eclatee1->Checked==true) Printer()->Canvas->Pen->Color=clOlive;
343	for (i=0;i<nb_tetra;i++)
344	{
345	if ((tetra[i].select==0)&& (Fenp->Debug1->Checked)) continue;
346	n1=tetra[i].n1;
347	n2=tetra[i].n2;
348	n3=tetra[i].n3;
349	n4=tetra[i].n4;
350	x1=point[n1].x;
351	y1=point[n1].y;
352	z1=point[n1].z;
353	x2=point[n2].x;
354	y2=point[n2].y;
355	z2=point[n2].z;
356	x3=point[n3].x;
357	y3=point[n3].y;
358	z3=point[n3].z;
359	x4=point[n4].x;
360	y4=point[n4].y;
361	z4=point[n4].z;
362	if (Normal1->Checked==false)
363	{
364	if (Qualite1->Checked==true)
365	{
366	crit=crit_cal3d(n1,n2,n3,n4);
367	if (crit>0.5) Printer()->Canvas->Pen->Color=clLime;
368	else if (crit>0.2) Printer()->Canvas->Pen->Color=clBlue;
369	else if (crit>0.1) Printer()->Canvas->Pen->Color=clRed;
370	else Printer()->Canvas->Pen->Color=clMaroon;
371	}
372	ab[0]=x2-x1;ab[1]=y2-y1;ab[2]=z2-z1;
373	ac[0]=x3-x1;ac[1]=y3-y1;ac[2]=z3-z1;
374	ad[0]=x4-x1;ad[1]=y4-y1;ad[2]=z4-z1;
375	g[0]=x1+0.25*(ab[0]+ac[0]+ad[0]);
376	g[1]=y1+0.25*(ab[1]+ac[1]+ad[1]);
377	g[2]=z1+0.25*(ab[2]+ac[2]+ad[2]);
378	x1=g[0]+0.8(x1-g[0]);y1=g[1]+0.8(y1-g[1]);z1=g[2]+0.8*(z1-g[2]);
379	x2=g[0]+0.8(x2-g[0]);y2=g[1]+0.8(y2-g[1]);z2=g[2]+0.8*(z2-g[2]);
380	x3=g[0]+0.8(x3-g[0]);y3=g[1]+0.8(y3-g[1]);z3=g[2]+0.8*(z3-g[2]);
381	x4=g[0]+0.8(x4-g[0]);y4=g[1]+0.8(y4-g[1]);z4=g[2]+0.8*(z4-g[2]);
382	}
383	projection(x1,y1,z1,&xc1,&yc1,2);
384	projection(x2,y2,z2,&xc2,&yc2,2);
385	projection(x3,y3,z3,&xc3,&yc3,2);
386	projection(x4,y4,z4,&xc4,&yc4,2);
387	Printer()->Canvas->MoveTo(xc1,yc1);
388	Printer()->Canvas->LineTo(xc2,yc2);
389	Printer()->Canvas->LineTo(xc3,yc3);
390	Printer()->Canvas->LineTo(xc1,yc1);
391	Printer()->Canvas->MoveTo(xc4,yc4);
392	Printer()->Canvas->LineTo(xc1,yc1);
393	Printer()->Canvas->MoveTo(xc4,yc4);
394	Printer()->Canvas->LineTo(xc2,yc2);
395	Printer()->Canvas->MoveTo(xc4,yc4);
396	Printer()->Canvas->LineTo(xc3,yc3);
397	Printer()->Canvas->Font->Color=clFuchsia;
398	if (Noeud1->Checked==true)
399	{
400	sprintf(mess,"%d",n1);
401	Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
402	sprintf(mess,"%d",n2);
403	Printer()->Canvas->TextOut(xc2+2,yc2+2,mess);
404	sprintf(mess,"%d",n3);
405	Printer()->Canvas->TextOut(xc3+2,yc3+2,mess);
406	sprintf(mess,"%d",n4);
407	Printer()->Canvas->TextOut(xc4+2,yc4+2,mess);
408	}
409	Printer()->Canvas->Font->Color=clYellow;
410	if (Maille2->Checked==true)
411	{
412	sprintf(mess,"%d",i);
413	Printer()->Canvas->TextOut((int)((xc1+xc2+xc3+xc4)0.25),(int)((yc1+yc2+yc3+yc4)0.25),mess);
414	}
415	}
416	if (numero_noeud>=0)
417	{
418	Printer()->Canvas->Font->Color=clBlack;
419	x1=point[numero_noeud].x;
420	y1=point[numero_noeud].y;
421	z1=point[numero_noeud].z;
422	projection(x1,y1,z1,&xc1,&yc1,2);
423	sprintf(mess,"%d",numero_noeud);
424	Printer()->Canvas->TextOut(xc1+2,yc1+2,mess);
425	}
426	if (numero_maille>=0)
427	{
428	Printer()->Canvas->Font->Color=clBlack;
429	n1=tetra[numero_maille].n1;
430	n2=tetra[numero_maille].n2;
431	n3=tetra[numero_maille].n3;
432	n4=tetra[numero_maille].n4;
433	x1=point[n1].x;
434	y1=point[n1].y;
435	z1=point[n1].z;
436	x2=point[n2].x;
437	y2=point[n2].y;
438	z2=point[n2].z;
439	x3=point[n3].x;
440	y3=point[n3].y;
441	z3=point[n3].z;
442	x4=point[n4].x;
443	y4=point[n4].y;
444	z4=point[n4].z;
445	projection(x1,y1,z1,&xc1,&yc1,2);
446	projection(x2,y2,z2,&xc2,&yc2,2);
447	projection(x3,y3,z3,&xc3,&yc3,2);
448	projection(x4,y4,z4,&xc4,&yc4,2);
449	sprintf(mess,"%d",numero_maille);
450	Printer()->Canvas->TextOut((int)((xc1+xc2+xc3+xc4)0.25),(int)((yc1+yc2+yc3+yc4)0.25),mess);
451	}
452
453
454	}
455	//---------------------------------------------------------------------------
456	void __fastcall TFenp::print_reel(void)
457	{
458	int *iordre;
459	float *distance;
460	int xpoint[8];
461	TPoint ypoint[8];
462	float xe,ye,ze,xee,yee,zee;
463	float xa,ya,za;
464	float xb,yb,zb;
465	float xc,yc,zc;
466	float xg,yg,zg;
467	float nx,ny,nz;
468	char mess[20];
469	struct s_triangle tete=NULL,tri;
470	int i;
471	int color;
472	float dismin,dismax,diff,g[4];
473	static float arot[3][3]={{-1.,-1.,-1.},{-1.,-1.,-1.},{-1.,-1.,-1.}};
474	int pareil;
475
476
477
478
479	pareil=0;
480	if (tete!=NULL)
481	if (arot[0][0]==rot[0][0])
482	if (arot[1][0]==rot[1][0])
483	if (arot[2][0]==rot[2][0])
484	if (arot[0][1]==rot[0][1])
485	if (arot[1][1]==rot[1][1])
486	if (arot[2][1]==rot[2][1])
487	if (arot[0][2]==rot[0][2])
488	if (arot[1][2]==rot[1][2])
489	if (arot[2][2]==rot[2][2]) pareil=1;
490
491
492
493	if (pareil==0)
494	{
495	tete=NULL;
496	arot[0][0]=rot[0][0];
497	arot[1][0]=rot[1][0];
498	arot[2][0]=rot[2][0];
499	arot[0][1]=rot[0][1];
500	arot[1][1]=rot[1][1];
501	arot[2][1]=rot[2][1];
502	arot[0][2]=rot[0][2];
503	arot[1][2]=rot[1][2];
504	arot[2][2]=rot[2][2];
505	/* oeuil */
506	xe=1000.maxcoord-0.5(xmin+xmax);
507	ye=1000.maxcoord-0.5(ymin+ymax);
508	ze=1000.maxcoord-0.5(zmin+zmax);
509	/* oeuil dans le repere local */
510	xee=rot[0][0]xe+rot[1][0]ye+rot[2][0]*ze;
511	yee=rot[0][1]xe+rot[1][1]ye+rot[2][1]*ze;
512	zee=rot[0][2]xe+rot[1][2]ye+rot[2][2]*ze;
513	iordre=(int *)calloc(nb_triangle+1,sizeof(int));
514	distance=(float *)calloc(nb_triangle+1,sizeof(float));
515	for (i=0;i<nb_triangle;i++)
516	{
517	triangle[i].dis=0.;
518	triangle[i].suivant=NULL;
519	xa=point[triangle[i].n1].x;
520	ya=point[triangle[i].n1].y;
521	za=point[triangle[i].n1].z;
522	xb=point[triangle[i].n2].x;
523	yb=point[triangle[i].n2].y;
524	zb=point[triangle[i].n2].z;
525	xc=point[triangle[i].n3].x;
526	yc=point[triangle[i].n3].y;
527	zc=point[triangle[i].n3].z;
528	xg=0.333333*(xa+xb+xc);
529	yg=0.333333*(ya+yb+yc);
530	zg=0.333333*(za+zb+zc);
531	triangle[i].dis=(xee-xg)(xee-xg)+(yee-yg)(yee-yg)+(zee-zg)*(zee-zg);
532	if (i==0)
533	{
534	dismin=triangle[i].dis;
535	dismax=triangle[i].dis;
536	}
537	else
538	{
539	if (triangle[i].dis<dismin) dismin=triangle[i].dis;
540	if (triangle[i].dis>dismax) dismax=triangle[i].dis;
541	}
542	nx=(yb-ya)(zc-za)-(zb-za)(yc-ya);
543	ny=(xc-xa)(zb-za)-(xb-xa)(zc-za);
544	nz=(xb-xa)(yc-ya)-(yb-ya)(xc-xa);
545	triangle[i].angle=(float)((nx(xee-xg)+ny(yee-yg)+nz(zee-zg))/sqrt((double)((nxnx+nyny+nznz)*(triangle[i].dis))));
546	if (triangle[i].angle>=1.) triangle[i].angle=0.99999999;
547	if (triangle[i].angle<=(-1.)) triangle[i].angle=(-0.99999999);
548	distance[i]=triangle[i].dis;
549	iordre[i]=i+1;
550
551	}
552	trirea(distance,&(nb_triangle),iordre);
553	diff=dismax-dismin;
554	dismin=dismin-0.025*diff;
555	dismax=dismax+0.025*diff;
556	Printer()->Canvas->Pen->Color=Cshade[0];
557	Printer()->Canvas->Brush->Color=Cshade[1];
558	for (i=nb_triangle-1;i>(-1);i--)
559	{
560	tri=&(triangle[iordre[i]-1]);
561	if ((tri->select==0)&& (Fenp->Debug1->Checked)) continue;
562	if (((tri->num_ent==numface) \|\| (numface==(-1))&&(tri->dis-disminloin-(1.-loin)dismax>0.)))
563	{
564	xa=point[tri->n1].x;
565	ya=point[tri->n1].y;
566	za=point[tri->n1].z;
567	xb=point[tri->n2].x;
568	yb=point[tri->n2].y;
569	zb=point[tri->n2].z;
570	xc=point[tri->n3].x;
571	yc=point[tri->n3].y;
572	zc=point[tri->n3].z;
573	if (Normal1->Checked==false)
574	{
575	g[0]=0.33333*(xa+xb+xc);
576	g[1]=0.33333*(ya+yb+yc);
577	g[2]=0.33333*(za+zb+zc);
578	xa=g[0]+homo(xa-g[0]);ya=g[1]+homo(ya-g[1]);za=g[2]+homo*(za-g[2]);
579	xb=g[0]+homo(xb-g[0]);yb=g[1]+homo(yb-g[1]);zb=g[2]+homo*(zb-g[2]);
580	xc=g[0]+homo(xc-g[0]);yc=g[1]+homo(yc-g[1]);zc=g[2]+homo*(zc-g[2]);
581	}
582	projection(xa,ya,za,&xpoint[0],&xpoint[1],2);
583	projection(xb,yb,zb,&xpoint[2],&xpoint[3],2);
584	projection(xc,yc,zc,&xpoint[4],&xpoint[5],2);
585	ypoint[0].x=xpoint[0];
586	ypoint[0].y=xpoint[1];
587	ypoint[1].x=xpoint[2];
588	ypoint[1].y=xpoint[3];
589	ypoint[2].x=xpoint[4];
590	ypoint[2].y=xpoint[5];
591	color=(int)(-10*fabs((double)(tri->angle))+10);
592	if (color>10) color=10;
593	if (color<1) color=1;
594	Printer()->Canvas->Brush->Color=Cshade[color];
595	Printer()->Canvas->Polygon(ypoint,2);
596	Printer()->Canvas->Font->Color=clFuchsia;
597	if (Noeud1->Checked==true)
598	{
599	sprintf(mess,"%d",tri->n1);
600	Printer()->Canvas->TextOut(xpoint[0]+2,xpoint[1]+2,mess);
601	sprintf(mess,"%d",tri->n2);
602	Printer()->Canvas->TextOut(xpoint[2]+2,xpoint[3]+2,mess);
603	sprintf(mess,"%d",tri->n3);
604	Printer()->Canvas->TextOut(xpoint[4]+2,xpoint[5]+2,mess);
605	}
606	if (tri->n1==numero_noeud)
607	{
608	Printer()->Canvas->Font->Color=clBlack;
609	sprintf(mess,"%d",tri->n1);
610	Printer()->Canvas->TextOut(xpoint[0]+2,xpoint[1]+2,mess);
611	}
612	if (tri->n2==numero_noeud)
613	{
614	Printer()->Canvas->Font->Color=clBlack;
615	sprintf(mess,"%d",tri->n2);
616	Printer()->Canvas->TextOut(xpoint[2]+2,xpoint[3]+2,mess);
617	}
618	if (tri->n3==numero_noeud)
619	{
620	Printer()->Canvas->Font->Color=clBlack;
621	sprintf(mess,"%d",tri->n3);
622	Printer()->Canvas->TextOut(xpoint[4]+2,xpoint[5]+2,mess);
623	}
624	if (tri->num!=numero_maille) Printer()->Canvas->Font->Color=clYellow;
625	else
626	{
627	Printer()->Canvas->Font->Color=clBlack;
628	sprintf(mess,"%d",tri->num);
629	Printer()->Canvas->TextOut((int)((xpoint[0]+xpoint[2]+xpoint[4])0.333333),(int)((xpoint[1]+xpoint[3]+xpoint[5])0.33333333),mess);
630	}
631	if (Maille2->Checked==true)
632	{
633	sprintf(mess,"%d",tri->num);
634	Printer()->Canvas->TextOut((int)((xpoint[0]+xpoint[2]+xpoint[4])0.333333),(int)((xpoint[1]+xpoint[3]+xpoint[5])0.33333333),mess);
635	}
636	}
637	}
638	free(iordre);
639	free(distance);
640	}
641	}
642	void __fastcall TFenp::print_axe(void)
643	{
644	int xc,yc;
645
646	Printer()->Canvas->Pen->Color=clBlue;
647	projection(75.,0.,0.,&xc,&yc,0);
648	Printer()->Canvas->MoveTo(100,100);
649	Printer()->Canvas->LineTo(100+xc,100+yc);
650	Printer()->Canvas->Pen->Color=clRed;
651	projection(0.,75.,0.,&xc,&yc,0);
652	Printer()->Canvas->MoveTo(100,100);
653	Printer()->Canvas->LineTo(100+xc,100+yc);
654	Printer()->Canvas->Pen->Color=clBlack;
655	projection(0.,0.,75.,&xc,&yc,0);
656	Printer()->Canvas->MoveTo(100,100);
657	Printer()->Canvas->LineTo(100+xc,100+yc);
658	}
659
660
661