apptools.c

Go to the documentation of this file.

/* =============================================================================
**  This file is part of the mmg software package for the tetrahedral
**  mesh modification.
**  Copyright (c) Bx INP/CNRS/Inria/UBordeaux/UPMC, 2004-
**
**  mmg 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 3 of the License, or
**  (at your option) any later version.
**
**  mmg 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 and of the GNU General Public License along with mmg (in
**  files COPYING.LESSER and COPYING). If not, see
**  <http://www.gnu.org/licenses/>. Please read their terms carefully and
**  use this copy of the mmg distribution only if you accept them.
** =============================================================================
*/
 
#include "mmgcommon_private.h"
 
static inline
int MMG5_Alloc_inode( MMG5_pMesh mesh, MMG5_iNode **node ) {
 
  MMG5_ADD_MEM(mesh,sizeof(MMG5_iNode),"boundary reference node",
                return 0;);
 
  MMG5_SAFE_MALLOC(*node,1,MMG5_iNode,return 0);
 
  return 1;
}
 
int MMG5_Add_inode( MMG5_pMesh mesh, MMG5_iNode **liLi, int val ) {
  MMG5_iNode  *newNode, *cur;
 
  cur = *liLi;
 
  /* Travel through the linked list and search if the value val exist or, if
   * not, where to insert it */
  if ( cur ) {
    if ( val < (*liLi)->val ) {
      /* Add a value at the list head */
      if ( !MMG5_Alloc_inode(mesh,&newNode) ) return -1;
 
      newNode->val = val;
      newNode->nxt = (*liLi);
 
      (*liLi) = newNode;
 
      return 1;
 
    }
    else if (val == (*liLi)->val ) return 0;
 
    while ( cur->nxt && ( val >= (cur->nxt)->val) )
      cur = cur->nxt;
 
    if ( val == cur->val ) return 0;
 
    if ( !MMG5_Alloc_inode(mesh,&newNode) ) return -1;
 
    newNode->val = val;
    newNode->nxt = cur->nxt;
    cur->nxt = newNode;
  }
  else {
    if ( !MMG5_Alloc_inode(mesh,&newNode) ) return -1;
 
    newNode->val = val;
    newNode->nxt = NULL;
 
    *liLi = newNode;
  }
 
  return 1;
}
 
void MMG5_Free_ilinkedList( MMG5_pMesh mesh, MMG5_iNode *liLi ) {
  MMG5_iNode *cur,*nxt;
 
  cur = liLi;
  while (cur) {
    nxt = cur;
    cur = cur->nxt;
 
    MMG5_DEL_MEM(mesh,nxt);
  }
}
 
int MMG5_countLocalParamAtTri( MMG5_pMesh mesh,MMG5_iNode **bdryRefs) {
  MMG5_int    k;
  int         npar,ier;
 
  (*bdryRefs) = NULL;
 
  k = mesh->nt? mesh->tria[1].ref : 0;
 
  /* Try to alloc the first node */
  ier = MMG5_Add_inode( mesh, bdryRefs, k );
  if ( ier < 0 ) {
    fprintf(stderr,"\n  ## Error: %s: unable to allocate the first boundary"
           " reference node.\n",__func__);
    return 0;
  }
  else {
    assert(ier);
    npar = 1;
  }
 
  for ( k=1; k<=mesh->nt; ++k ) {
    ier = MMG5_Add_inode( mesh, bdryRefs, mesh->tria[k].ref );
 
    if ( ier < 0 ) {
      printf("  ## Warning: %s: unable to list the tria references."
             " Uncomplete parameters file.\n",__func__ );
      break;
    }
    else if ( ier ) ++npar;
  }
 
  return npar;
}
 
int MMG5_writeLocalParamAtTri( MMG5_pMesh mesh, MMG5_iNode *bdryRefs,
                                FILE *out ) {
  MMG5_iNode *cur;
 
  cur = bdryRefs;
  while( cur ) {
    fprintf(out,"%"MMG5_PRId" Triangle %e %e %e \n",cur->val,
            mesh->info.hmin, mesh->info.hmax,mesh->info.hausd);
    cur = cur->nxt;
  }
 
  MMG5_Free_ilinkedList(mesh,bdryRefs);
 
  return 1;
}