2003-07-10 15:06:41 +06:00
|
|
|
// SMESH MEFISTO2 : algorithm for meshing
|
|
|
|
|
//
|
|
|
|
|
// Copyright (C) 2003 Laboratoire J.-L. Lions UPMC Paris
|
|
|
|
|
//
|
|
|
|
|
// This library is free software; you can redistribute it and/or
|
|
|
|
|
// modify it under the terms of the GNU Lesser General Public
|
|
|
|
|
// License as published by the Free Software Foundation; either
|
|
|
|
|
// version 2.1 of the License.
|
|
|
|
|
//
|
|
|
|
|
// This library 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
|
|
|
|
|
// Lesser General Public License for more details.
|
|
|
|
|
//
|
|
|
|
|
// You should have received a copy of the GNU Lesser General Public
|
|
|
|
|
// License along with this library; if not, write to the Free Software
|
|
|
|
|
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
|
|
|
|
|
//
|
|
|
|
|
// See http://www.ann.jussieu.fr/~perronne or email Perronnet@ann.jussieu.fr
|
|
|
|
|
//
|
|
|
|
|
//
|
|
|
|
|
//
|
|
|
|
|
// File : aptrte.h
|
|
|
|
|
// Author: Alain PERRONNET
|
|
|
|
|
// Module : SMESH
|
|
|
|
|
|
2003-05-19 20:07:00 +06:00
|
|
|
#ifndef aptrte__h
|
|
|
|
|
#define aptrte__h
|
|
|
|
|
|
|
|
|
|
#include <limits.h> // limites min max int long real ...
|
|
|
|
|
#include <unistd.h> // gethostname, ...
|
|
|
|
|
#include <stdio.h>
|
|
|
|
|
#include <iostream.h> // pour cout cin ...
|
|
|
|
|
#include <iomanip.h> // pour le format des io setw, stx, setfill, ...
|
|
|
|
|
#include <string.h> // pour les fonctions sur les chaines de caracteres
|
|
|
|
|
#include <ctype.h>
|
|
|
|
|
#include <stdlib.h>
|
|
|
|
|
#include <math.h> // pour les fonctions mathematiques
|
|
|
|
|
#include <time.h>
|
|
|
|
|
|
|
|
|
|
#include <sys/types.h>
|
|
|
|
|
#include <sys/time.h>
|
|
|
|
|
|
|
|
|
|
void qualitetrte( R3 *mnpxyd,
|
|
|
|
|
Z & mosoar, Z & mxsoar, Z *mnsoar,
|
|
|
|
|
Z & moartr, Z & mxartr, Z *mnartr,
|
|
|
|
|
Z & nbtria, R & quamoy, R & quamin );
|
|
|
|
|
// +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
|
|
|
|
|
// but : calculer la qualite moyenne et minimale de la triangulation
|
|
|
|
|
// ----- actuelle definie par les tableaux nosoar et noartr
|
|
|
|
|
// entrees:
|
|
|
|
|
// --------
|
|
|
|
|
// mnpxyd : tableau des coordonnees 2d des points
|
|
|
|
|
// par point : x y distance_souhaitee
|
|
|
|
|
// mosoar : nombre maximal d'entiers par arete et
|
|
|
|
|
// indice dans nosoar de l'arete suivante dans le hachage
|
|
|
|
|
// mxsoar : nombre maximal d'aretes stockables dans le tableau nosoar
|
|
|
|
|
// attention: mxsoar>3*mxsomm obligatoire!
|
|
|
|
|
// nosoar : numero des 2 sommets , no ligne, 2 triangles de l'arete,
|
|
|
|
|
// chainage des aretes frontalieres, chainage du hachage des aretes
|
|
|
|
|
// hachage des aretes = nosoar(1)+nosoar(2)*2
|
|
|
|
|
// avec mxsoar>=3*mxsomm
|
|
|
|
|
// une arete i de nosoar est vide <=> nosoar(1,i)=0 et
|
|
|
|
|
// nosoar(2,arete vide)=l'arete vide qui precede
|
|
|
|
|
// nosoar(3,arete vide)=l'arete vide qui suit
|
|
|
|
|
// moartr : nombre maximal d'entiers par arete du tableau noartr
|
|
|
|
|
// mxartr : nombre maximal de triangles declarables
|
|
|
|
|
// noartr : les 3 aretes des triangles +-arete1, +-arete2, +-arete3
|
|
|
|
|
// arete1 = 0 si triangle vide => arete2 = triangle vide suivant
|
|
|
|
|
// sorties:
|
|
|
|
|
// --------
|
|
|
|
|
// nbtria : nombre de triangles internes au domaine
|
|
|
|
|
// quamoy : qualite moyenne des triangles actuels
|
|
|
|
|
// quamin : qualite minimale des triangles actuels
|
|
|
|
|
// +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
void aptrte( Z nutysu, R aretmx,
|
|
|
|
|
Z nblf, Z *nudslf, R2 *uvslf,
|
|
|
|
|
Z nbpti, R2 *uvpti,
|
|
|
|
|
Z & nbst, R2 * & uvst, Z & nbt, Z * & nust,
|
|
|
|
|
Z & ierr );
|
|
|
|
|
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
|
|
|
|
|
// but : appel de la triangulation par un arbre-4 recouvrant
|
|
|
|
|
// ----- de triangles equilateraux
|
|
|
|
|
// le contour du domaine plan est defini par des lignes fermees
|
|
|
|
|
// la premiere ligne etant l'enveloppe de toutes les autres
|
|
|
|
|
// la fonction areteideale_(s,d) donne la taille d'arete
|
|
|
|
|
// au point s dans la direction d (direction inactive pour l'instant)
|
|
|
|
|
// des lors toute arete issue d'un sommet s devrait avoir une longueur
|
|
|
|
|
// comprise entre 0.65 areteideale_(s,d) et 1.3 areteideale_(s,d)
|
|
|
|
|
//
|
|
|
|
|
//Attention:
|
|
|
|
|
// Les tableaux uvslf et uvpti sont supposes ne pas avoir de sommets identiques!
|
|
|
|
|
// De meme, un sommet d'une ligne fermee ne peut appartenir a une autre ligne fermee
|
|
|
|
|
//
|
|
|
|
|
// entrees:
|
|
|
|
|
// --------
|
|
|
|
|
// nutysu : numero de traitement de areteideale_() selon le type de surface
|
|
|
|
|
// 0 pas d'emploi de la fonction areteideale_() et aretmx est active
|
|
|
|
|
// 1 il existe une fonction areteideale_(s,d)
|
|
|
|
|
// dont seules les 2 premieres composantes de uv sont actives
|
|
|
|
|
// ... autres options a definir ...
|
|
|
|
|
// aretmx : longueur maximale des aretes de la future triangulation
|
|
|
|
|
// nblf : nombre de lignes fermees de la surface
|
|
|
|
|
// nudslf : numero du dernier sommet de chacune des nblf lignes fermees
|
|
|
|
|
// nudslf(0)=0 pour permettre la difference sans test
|
|
|
|
|
// Attention le dernier sommet de chaque ligne est raccorde au premier
|
|
|
|
|
// tous les sommets et les points internes ont des coordonnees
|
|
|
|
|
// UV differentes <=> Pas de point double!
|
|
|
|
|
// uvslf : uv des nudslf(nblf) sommets des lignes fermees
|
|
|
|
|
// nbpti : nombre de points internes futurs sommets de la triangulation
|
|
|
|
|
// uvpti : uv des points internes futurs sommets de la triangulation
|
|
|
|
|
//
|
|
|
|
|
// sorties:
|
|
|
|
|
// --------
|
|
|
|
|
// nbst : nombre de sommets de la triangulation finale
|
|
|
|
|
// uvst : coordonnees uv des nbst sommets de la triangulation
|
|
|
|
|
// nbt : nombre de triangles de la triangulation finale
|
|
|
|
|
// nust : 3 numeros dans uvst des sommets des nbt triangles
|
|
|
|
|
// ierr : 0 si pas d'erreur
|
|
|
|
|
// > 0 sinon
|
|
|
|
|
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
|
|
|
|
|
// auteur : Alain Perronnet Analyse Numerique Paris UPMC decembre 2001
|
|
|
|
|
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
|
|
|
|
|
|
|
|
|
|
extern "C" { void tempscpu_( double & tempsec ); }
|
|
|
|
|
//Retourne le temps CPU utilise en secondes
|
|
|
|
|
|
|
|
|
|
extern "C" { void deltacpu_( R & dtcpu ); }
|
|
|
|
|
//Retourne le temps CPU utilise en secondes depuis le precedent appel
|
|
|
|
|
|
|
|
|
|
//initialiser le tableau mnsoar pour le hachage des aretes
|
|
|
|
|
extern "C" {void insoar_( Z & mxsomm, Z & mosoar, Z & mxsoar, Z & n1soar, Z * mnsoar );}
|
|
|
|
|
|
|
|
|
|
//mettre a zero les nb entiers de tab
|
|
|
|
|
extern "C" {void azeroi_( Z & nb, Z * tab );}
|
|
|
|
|
|
|
|
|
|
extern "C" {void fasoar_( Z & ns1, Z & ns2, Z & nt1, Z & nt2, Z & nolign,
|
|
|
|
|
Z & mosoar, Z & mxsoar, Z & n1soar, Z * mnsoar, Z * mnarst,
|
|
|
|
|
Z & noar, Z & ierr );}
|
|
|
|
|
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
|
|
|
|
|
// but : former l'arete de sommet ns1-ns2 dans le hachage du tableau
|
|
|
|
|
// ----- nosoar des aretes de la triangulation
|
|
|
|
|
// entrees:
|
|
|
|
|
// --------
|
|
|
|
|
// ns1 ns2: numero pxyd des 2 sommets de l'arete
|
|
|
|
|
// nt1 : numero du triangle auquel appartient l'arete
|
|
|
|
|
// nt1=-1 si numero inconnu
|
|
|
|
|
// nt2 : numero de l'eventuel second triangle de l'arete si connu
|
|
|
|
|
// nt2=-1 si numero inconnu
|
|
|
|
|
// nolign : numero de la ligne fermee de l'arete
|
|
|
|
|
// =0 si l'arete n'est une arete de ligne
|
|
|
|
|
// ce numero est ajoute seulement si l'arete est creee
|
|
|
|
|
// mosoar : nombre maximal d'entiers par arete du tableau nosoar
|
|
|
|
|
// mxsoar : nombre maximal d'aretes stockables dans le tableau nosoar
|
|
|
|
|
// modifies:
|
|
|
|
|
// ---------
|
|
|
|
|
// n1soar : numero de la premiere arete vide dans le tableau nosoar
|
|
|
|
|
// une arete i de nosoar est vide <=> nosoar(1,i)=0
|
|
|
|
|
// chainage des aretes vides amont et aval
|
|
|
|
|
// l'arete vide qui precede=nosoar(4,i)
|
|
|
|
|
// l'arete vide qui suit =nosoar(5,i)
|
|
|
|
|
// nosoar : numero des 2 sommets, no ligne, 2 triangles de l'arete,
|
|
|
|
|
// chainage momentan'e d'aretes, chainage du hachage des aretes
|
|
|
|
|
// hachage des aretes = min( nosoar(1), nosoar(2) )
|
|
|
|
|
// noarst : noarst(np) numero d'une arete du sommet np
|
|
|
|
|
|
|
|
|
|
// ierr : si < 0 en entree pas d'affichage en cas d'erreur du type
|
|
|
|
|
// "arete appartenant a plus de 2 triangles et a creer!"
|
|
|
|
|
// si >=0 en entree affichage de ce type d'erreur
|
|
|
|
|
// sorties:
|
|
|
|
|
// --------
|
|
|
|
|
// noar : >0 numero de l'arete retrouvee ou ajoutee
|
|
|
|
|
// ierr : =0 si pas d'erreur
|
|
|
|
|
// =1 si le tableau nosoar est sature
|
|
|
|
|
// =2 si arete a creer et appartenant a 2 triangles distincts
|
|
|
|
|
// des triangles nt1 et nt2
|
|
|
|
|
// =3 si arete appartenant a 2 triangles distincts
|
|
|
|
|
// differents des triangles nt1 et nt2
|
|
|
|
|
// =4 si arete appartenant a 2 triangles distincts
|
|
|
|
|
// dont le second n'est pas le triangle nt2
|
|
|
|
|
//+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
|
|
|
|
|
|
|
|
|
|
//initialisation du tableau letree et ajout dans letree des sommets 1 a nbsomm
|
|
|
|
|
extern "C" {void teajte_( Z & mxsomm, Z & nbsomm, R3 * mnpxyd, R3 * comxmi,
|
|
|
|
|
R & aretmx, Z & mxtree, Z * letree,
|
|
|
|
|
Z & ierr );}
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
extern "C" {void tehote_( Z & nutysu, Z & nbarpi, Z & mxsomm, Z & nbsomm, R3 * mnpxyd,
|
|
|
|
|
R3 * comxmi, R & aretmx,
|
|
|
|
|
Z * letree, Z & mxqueu, Z * mnqueu,
|
|
|
|
|
Z & ierr );}
|
|
|
|
|
// homogeneisation de l'arbre des te a un saut de taille au plus
|
|
|
|
|
// prise en compte des tailles d'aretes souhaitees autour des sommets initiaux
|
|
|
|
|
|
|
|
|
|
extern "C" {void tetrte_( R3 * comxmi, R & aretmx, Z & nbarpi, Z & mxsomm, R3 * mnpxyd,
|
|
|
|
|
Z & mxqueu, Z * mnqueu, Z * mntree,
|
|
|
|
|
Z & mosoar, Z & mxsoar, Z & n1soar, Z * mnsoar,
|
|
|
|
|
Z & moartr, Z & mxartr, Z & n1artr, Z * mnartr, Z * mnarst,
|
|
|
|
|
Z & ierr );}
|
|
|
|
|
// trianguler les triangles equilateraux feuilles a partir de leurs 3 sommets
|
|
|
|
|
// et des points de la frontiere, des points internes imposes interieurs
|
|
|
|
|
|
|
|
|
|
extern "C" {void aisoar_( Z & mosoar, Z & mxsoar, Z * mnsoar, Z & na );}
|
|
|
|
|
// formation du chainage 6 des aretes internes a echanger eventuellement
|
|
|
|
|
|
|
|
|
|
extern "C" {void tedela_( R3 * mnpxyd, Z * mnarst,
|
|
|
|
|
Z & mosoar, Z & mxsoar, Z & n1soar, Z * mnsoar, Z & na,
|
|
|
|
|
Z & moartr, Z & mxartr, Z & n1artr, Z * mnartr, Z & n );}
|
|
|
|
|
// boucle sur les aretes internes (non sur une ligne de la frontiere)
|
|
|
|
|
// avec echange des 2 diagonales afin de rendre la triangulation delaunay
|
|
|
|
|
|
|
|
|
|
extern "C" {void terefr_( Z & nbarpi, R3 * mnpxyd,
|
|
|
|
|
Z & mosoar, Z & mxsoar, Z & n1soar, Z * mnsoar,
|
|
|
|
|
Z & moartr, Z & n1artr, Z * mnartr, Z * mnarst,
|
|
|
|
|
Z & mxarcf, Z * mnarc1, Z * mnarc2,
|
|
|
|
|
Z * mnarc3, Z * mnarc4,
|
|
|
|
|
Z & n, Z & ierr );}
|
|
|
|
|
// detection des aretes frontalieres initiales perdues
|
|
|
|
|
// triangulation frontale pour les restaurer
|
|
|
|
|
|
|
|
|
|
extern "C" {void tesuex_( Z & nblf, Z * nulftr,
|
|
|
|
|
Z & ndtri0, Z & nbsomm, R3 * mnpxyd, Z * mnslig,
|
|
|
|
|
Z & mosoar, Z & mxsoar, Z * mnsoar,
|
|
|
|
|
Z & moartr, Z & mxartr, Z & n1artr, Z * mnartr, Z * mnarst,
|
|
|
|
|
Z & nbtria, Z * mntrsu, Z & ierr );}
|
|
|
|
|
// suppression des triangles externes a la surface
|
|
|
|
|
|
|
|
|
|
extern "C" {void teamqt_( Z & nutysu,
|
|
|
|
|
Z * mnarst, Z & mosoar, Z & mxsoar, Z & n1soar, Z * mnsoar,
|
|
|
|
|
Z & moartr, Z & mxartr, Z & n1artr, Z * mnartr,
|
|
|
|
|
Z & mxarcf, Z * mntrcf, Z * mnstbo,
|
|
|
|
|
Z * n1arcf, Z * mnarcf, Z * mnarc1,
|
|
|
|
|
R3 * comxmi, Z & nbarpi, Z & nbsomm, Z & mxsomm,
|
|
|
|
|
R3 * mnpxyd, Z * mnslig,
|
|
|
|
|
Z & ierr );}
|
|
|
|
|
// amelioration de la qualite de la triangulation par
|
|
|
|
|
// barycentrage des sommets internes a la triangulation
|
|
|
|
|
// suppression des aretes trop longues ou trop courtes
|
|
|
|
|
// modification de la topologie des groupes de triangles
|
|
|
|
|
// mise en delaunay de la triangulation
|
|
|
|
|
|
|
|
|
|
extern "C" {void nusotr_( Z & nt, Z & mosoar, Z * mnsoar, Z & moartr, Z * mnartr,
|
|
|
|
|
Z * nosotr );}
|
|
|
|
|
//retrouver les numero des 3 sommets du triangle nt
|
|
|
|
|
|
|
|
|
|
extern "C" {void qutr2d_( R3 & p1, R3 & p2, R3 & p3, R & qualite );}
|
|
|
|
|
//calculer la qualite d'un triangle de R2 de sommets p1, p2, p3
|
|
|
|
|
|
|
|
|
|
extern "C" { R surtd2_( R3 & p1, R3 & p2, R3 & p3 ); }
|
|
|
|
|
//calcul de la surface d'un triangle defini par 3 points de r**2
|
|
|
|
|
|
|
|
|
|
#endif
|
