Changeset 565 in Sophya for trunk/SophyaLib/Samba/utilgeom.cc
- Timestamp:
- Nov 10, 1999, 3:01:46 PM (26 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaLib/Samba/utilgeom.cc
r508 r565 2 2 #include <iostream.h> 3 3 #include "utilgeom.h" 4 4 //++ 5 // Titre Some utilitary functions for geometry (utilgeom.h)... 6 //-- 7 //++ 5 8 int sign(double d) 9 // 10 //-- 6 11 { 7 12 return (d >= 0) - (d < 0); 8 13 } 9 14 //++ 10 15 double absd(double d) 16 // 17 //-- 11 18 { 12 19 return sqrt(d*d); 13 20 } 14 21 //++ 15 22 double mod(double d, double periode) 23 // 24 //-- 16 25 { 17 26 if( d >= 0 ) return d-floor(d/periode)*periode; 18 27 else return d-ceil(d/periode)*periode+periode; 19 28 } 20 29 //++ 21 30 void swap(double& d1, double& d2) 31 // 32 //-- 22 33 { 23 34 double temp; … … 26 37 d1=temp; 27 38 } 28 39 //++ 29 40 double min(double d1, double d2) 41 // 42 //-- 30 43 { 31 44 if( d1 >= d2 ) return d2; 32 45 else return d1; 33 46 } 34 47 //++ 35 48 double max(double d1, double d2) 49 // 50 //-- 36 51 { 37 52 return -min(-d1,-d2); 38 53 } 39 54 //++ 40 55 int arrondi(double d) 56 // 57 //-- 41 58 { 42 59 return (int)(((d-floor(d)) >= 0.5)*ceil(d)+((d-floor(d)) < 0.5)*floor(d)); 43 60 } 44 61 //++ 45 62 long rangijd(int nc, int i, int j, int d) 63 // 64 //-- 46 65 { 47 66 if( i < j ) return 2*i*(2*nc-i-1)+4*(j-i-1)+d; … … 55 74 else return -1; 56 75 } 57 76 //++ 58 77 long rangdiff(int nc, int i, int j, int d) 78 // 79 //-- 59 80 { 60 81 if( i == j ) … … 68 89 return indm*(2*nc-indm-1)+2*(indM-indm-1)+d; 69 90 } 70 91 //++ 71 92 long rangik(int NtotEch, int i, int k) 93 // 94 //-- 72 95 { 73 96 return NtotEch*i+k; 74 97 } 75 98 //++ 76 99 long ranghk(int NtotEch, int h, int k) 100 // 101 //-- 77 102 { 78 103 return NtotEch*h+k; 79 104 } 80 105 //++ 81 106 double scangle(double sinus, double cosinus) 107 // 108 //-- 82 109 { 83 110 double eps=1e-10;
Note:
See TracChangeset
for help on using the changeset viewer.