Changeset 565 in Sophya for trunk/SophyaLib/Samba/utilgeom.cc


Ignore:
Timestamp:
Nov 10, 1999, 3:01:46 PM (26 years ago)
Author:
ansari
Message:

ajout doc GLM

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/SophyaLib/Samba/utilgeom.cc

    r508 r565  
    22#include <iostream.h>
    33#include "utilgeom.h"
    4 
     4//++
     5// Titre        Some utilitary functions for geometry (utilgeom.h)...
     6//--
     7//++
    58int sign(double d)
     9//
     10//--
    611{
    712  return (d >= 0) - (d < 0);
    813}
    9 
     14//++
    1015double absd(double d)
     16//
     17//--
    1118{
    1219  return sqrt(d*d);
    1320}
    14 
     21//++
    1522double mod(double d, double periode)
     23//
     24//--
    1625{
    1726  if( d >= 0 ) return d-floor(d/periode)*periode;
    1827  else return d-ceil(d/periode)*periode+periode;
    1928}
    20 
     29//++
    2130void swap(double& d1, double& d2)
     31//
     32//--
    2233{
    2334  double temp;
     
    2637  d1=temp;
    2738}
    28 
     39//++
    2940double min(double d1, double d2)
     41//
     42//--
    3043{
    3144  if( d1 >= d2 ) return d2;
    3245  else return d1;
    3346}
    34 
     47//++
    3548double max(double d1, double d2)
     49//
     50//--
    3651{
    3752  return -min(-d1,-d2);
    3853}
    39 
     54//++
    4055int arrondi(double d)
     56//
     57//--
    4158{
    4259  return (int)(((d-floor(d)) >= 0.5)*ceil(d)+((d-floor(d)) < 0.5)*floor(d));
    4360}
    44 
     61//++
    4562long rangijd(int nc, int i, int j, int d)
     63//
     64//--
    4665{
    4766  if( i < j ) return 2*i*(2*nc-i-1)+4*(j-i-1)+d;
     
    5574  else return -1;
    5675}
    57 
     76//++
    5877long rangdiff(int nc, int i, int j, int d)
     78//
     79//--
    5980{
    6081  if( i == j )
     
    6889  return indm*(2*nc-indm-1)+2*(indM-indm-1)+d;
    6990}
    70 
     91//++
    7192long rangik(int NtotEch, int i, int k)
     93//
     94//--
    7295{
    7396  return NtotEch*i+k;
    7497}
    75 
     98//++
    7699long ranghk(int NtotEch, int h, int k)
     100//
     101//--
    77102{
    78103  return NtotEch*h+k;
    79104}
    80 
     105//++
    81106double scangle(double sinus, double cosinus)
     107//
     108//--
    82109{
    83110  double eps=1e-10;
Note: See TracChangeset for help on using the changeset viewer.