Changeset 3615 in Sophya for trunk/SophyaLib/NTools


Ignore:
Timestamp:
May 1, 2009, 1:34:31 PM (16 years ago)
Author:
cmv
Message:

Modifs relatives a l'introduction de RandomGeneratorInterface + delete de srandgen.c, cmv 01/05/2009

Location:
trunk/SophyaLib/NTools
Files:
5 edited

Legend:

Unmodified
Added
Removed

  • trunk/SophyaLib/NTools/Makefile

    r3608 r3615  
    3333  $(SOPHYAINCP)ndatablock.h \
    3434  $(SOPHYAINCP)utilarr.h \
     35  $(SOPHYAINCP)randinterf.h \
    3536  $(SOPHYAINCP)tmatrix.h \
    3637  $(SOPHYAINCP)tmatrix_tsnl.h \
    … …  
    6768  $(SOPHYAINCP)ndatablock.h \
    6869  $(SOPHYAINCP)utilarr.h \
     70  $(SOPHYAINCP)randinterf.h \
    6971  $(SOPHYAINCP)tmatrix_tsnl.h \
    7072  $(SOPHYAINCP)fioarr.h \
    … …  
    114116  $(SOPHYAINCP)ndatablock.h \
    115117  $(SOPHYAINCP)utilarr.h \
     118  $(SOPHYAINCP)randinterf.h \
    116119  $(SOPHYAINCP)tmatrix_tsnl.h \
    117120  $(SOPHYAINCP)tvector_tsnl.h \
    … …  
    120123 
    121124$(SOPHYAOBJP)dynccd.o: dynccd.cc $(SOPHYAINCP)sopnamsp.h \
    122   $(SOPHYAINCP)machdefs.h fmath.h perandom.h \
    123   $(SOPHYAINCP)tvector.h \
    124   $(SOPHYAINCP)tmatrix.h \
    125   $(SOPHYAINCP)machdefs.h \
    126   $(SOPHYAINCP)tarray.h \
    127   $(SOPHYAINCP)basarr.h \
    128   $(SOPHYAINCP)anydataobj.h \
     125  $(SOPHYAINCP)machdefs.h fmath.h \
     126  $(SOPHYAINCP)srandgen.h \
     127  $(SOPHYAINCP)machdefs.h \
     128  $(SOPHYAINCP)randinterf.h \
     129  $(SOPHYAINCP)anydataobj.h perandom.h \
     130  $(SOPHYAINCP)tvector.h \
     131  $(SOPHYAINCP)tmatrix.h \
     132  $(SOPHYAINCP)tarray.h \
     133  $(SOPHYAINCP)basarr.h \
    129134  $(SOPHYAINCP)mutyv.h \
    130135  $(SOPHYAINCP)timestamp.h \
    … …  
    146151  $(SOPHYAINCP)perrors.h \
    147152  $(SOPHYAINCP)fmath.h \
    148   $(SOPHYAINCP)tvector.h \
    149   $(SOPHYAINCP)srandgen.h classfunc.h cimage.h \
     153  $(SOPHYAINCP)tvector.h classfunc.h cimage.h \
    150154  $(SOPHYAINCP)tmatrix.h \
    151155  $(SOPHYAINCP)fioarr.h dynccd.h \
    … …  
    175179  $(SOPHYAINCP)ndatablock.h \
    176180  $(SOPHYAINCP)utilarr.h \
     181  $(SOPHYAINCP)randinterf.h \
    177182  $(SOPHYAINCP)tmatrix_tsnl.h \
    178183  $(SOPHYAINCP)tvector_tsnl.h generaldata.h \
    … …  
    207212  $(SOPHYAINCP)ndatablock.h \
    208213  $(SOPHYAINCP)utilarr.h \
     214  $(SOPHYAINCP)randinterf.h \
    209215  $(SOPHYAINCP)tmatrix_tsnl.h \
    210216  $(SOPHYAINCP)tvector_tsnl.h generaldata.h \
    … …  
    239245  $(SOPHYAINCP)ndatablock.h \
    240246  $(SOPHYAINCP)utilarr.h \
     247  $(SOPHYAINCP)randinterf.h \
    241248  $(SOPHYAINCP)tmatrix_tsnl.h \
    242249  $(SOPHYAINCP)tvector.h \
    … …  
    266273  $(SOPHYAINCP)ndatablock.h \
    267274  $(SOPHYAINCP)utilarr.h \
     275  $(SOPHYAINCP)randinterf.h \
    268276  $(SOPHYAINCP)tmatrix_tsnl.h \
    269277  $(SOPHYAINCP)tvector.h \
    … …  
    292300  $(SOPHYAINCP)ndatablock.h \
    293301  $(SOPHYAINCP)utilarr.h \
     302  $(SOPHYAINCP)randinterf.h \
    294303  $(SOPHYAINCP)tmatrix_tsnl.h \
    295304  $(SOPHYAINCP)tvector.h \
    … …  
    319328  $(SOPHYAINCP)ndatablock.h \
    320329  $(SOPHYAINCP)utilarr.h \
     330  $(SOPHYAINCP)randinterf.h \
    321331  $(SOPHYAINCP)tmatrix_tsnl.h \
    322332  $(SOPHYAINCP)tvector_tsnl.h \
    … …  
    346356  $(SOPHYAINCP)ndatablock.h \
    347357  $(SOPHYAINCP)utilarr.h \
     358  $(SOPHYAINCP)randinterf.h \
    348359  $(SOPHYAINCP)tmatrix_tsnl.h \
    349360  $(SOPHYAINCP)tvector_tsnl.h generaldata.h \
    … …  
    378389  $(SOPHYAINCP)ndatablock.h \
    379390  $(SOPHYAINCP)utilarr.h \
     391  $(SOPHYAINCP)randinterf.h \
    380392  $(SOPHYAINCP)tmatrix_tsnl.h \
    381393  $(SOPHYAINCP)tvector_tsnl.h generaldata.h \
    … …  
    418430  $(SOPHYAINCP)ndatablock.h \
    419431  $(SOPHYAINCP)utilarr.h \
     432  $(SOPHYAINCP)randinterf.h \
    420433  $(SOPHYAINCP)tmatrix_tsnl.h \
    421434  $(SOPHYAINCP)tvector_tsnl.h generaldata.h \
    … …  
    458471  $(SOPHYAINCP)ndatablock.h \
    459472  $(SOPHYAINCP)utilarr.h \
     473  $(SOPHYAINCP)randinterf.h \
    460474  $(SOPHYAINCP)tmatrix_tsnl.h \
    461475  $(SOPHYAINCP)tvector_tsnl.h \
    … …  
    488502  $(SOPHYAINCP)ndatablock.h \
    489503  $(SOPHYAINCP)utilarr.h \
     504  $(SOPHYAINCP)randinterf.h \
    490505  $(SOPHYAINCP)tmatrix_tsnl.h cimage.h \
    491506  $(SOPHYAINCP)fioarr.h \
    … …  
    512527  $(SOPHYAINCP)machdefs.h \
    513528  $(SOPHYAINCP)pexceptions.h \
    514   $(SOPHYAINCP)machdefs.h perandom.h \
    515   $(SOPHYAINCP)tvector.h \
    516   $(SOPHYAINCP)tmatrix.h \
    517   $(SOPHYAINCP)tarray.h \
    518   $(SOPHYAINCP)basarr.h \
    519   $(SOPHYAINCP)anydataobj.h \
     529  $(SOPHYAINCP)machdefs.h \
     530  $(SOPHYAINCP)srandgen.h \
     531  $(SOPHYAINCP)randinterf.h \
     532  $(SOPHYAINCP)anydataobj.h perandom.h \
     533  $(SOPHYAINCP)tvector.h \
     534  $(SOPHYAINCP)tmatrix.h \
     535  $(SOPHYAINCP)tarray.h \
     536  $(SOPHYAINCP)basarr.h \
    520537  $(SOPHYAINCP)mutyv.h \
    521538  $(SOPHYAINCP)timestamp.h \
    … …  
    537554  $(SOPHYAINCP)perrors.h \
    538555  $(SOPHYAINCP)fmath.h \
    539   $(SOPHYAINCP)tvector.h \
    540   $(SOPHYAINCP)srandgen.h classfunc.h pemath.h
     556  $(SOPHYAINCP)tvector.h classfunc.h
    541557        $(CXXCOMPILE) $(CXXTEMPFLG) -o $@  perandom.cc
    542558 
    … …  
    564580  $(SOPHYAINCP)ndatablock.h \
    565581  $(SOPHYAINCP)utilarr.h \
     582  $(SOPHYAINCP)randinterf.h \
    566583  $(SOPHYAINCP)tmatrix_tsnl.h \
    567584  $(SOPHYAINCP)tvector_tsnl.h \
    … …  
    594611  $(SOPHYAINCP)ndatablock.h \
    595612  $(SOPHYAINCP)utilarr.h \
     613  $(SOPHYAINCP)randinterf.h \
    596614  $(SOPHYAINCP)tmatrix_tsnl.h \
    597615  $(SOPHYAINCP)tvector_tsnl.h \
    … …  
    629647  $(SOPHYAINCP)ndatablock.h \
    630648  $(SOPHYAINCP)utilarr.h \
     649  $(SOPHYAINCP)randinterf.h \
    631650  $(SOPHYAINCP)tmatrix_tsnl.h \
    632651  $(SOPHYAINCP)tvector_tsnl.h generalfit.h \
    … …  
    680699  $(SOPHYAINCP)ndatablock.h \
    681700  $(SOPHYAINCP)utilarr.h \
     701  $(SOPHYAINCP)randinterf.h \
    682702  $(SOPHYAINCP)tmatrix.h \
    683703  $(SOPHYAINCP)tmatrix_tsnl.h \

  • trunk/SophyaLib/NTools/dynccd.cc

    r2882 r3615  
    55
    66#include "fmath.h"
     7#include "srandgen.h"
    78#include "perandom.h"
    89

  • trunk/SophyaLib/NTools/perandom.cc

    r3608 r3615  
    11#include "sopnamsp.h"
    22#include "machdefs.h"
     3#include <iostream>
    34#include "pexceptions.h"
     5#include "srandgen.h"
    46#include "perandom.h"
    5 #include "pemath.h"
    6 #include <iostream>
    77
    88
    … …  
    290290}
    291291
    292 
    293 
    294 /////////////////////////////////////////////////////////////////
    295 /*
    296 **** Remarques sur complex< r_8 > ComplexGaussRan(double sig) ****
    297 
    298 --- variables gaussiennes x,y independantes
    299 x gaussien: pdf f(x) = 1/(sqrt(2Pi) Sx) exp(-(x-Mx)^2/(2 Sx^2))
    300 y gaussien: pdf f(y) = 1/(sqrt(2Pi) Sy) exp(-(y-My)^2/(2 Sy^2))
    301 x,y independants --> pdf f(x,y) = f(x) f(y)
    302 On a:
    303   <x>   = Integrate[x*f(x)]   = Mx
    304   <x^2> = Integrate[x^2*f(x)] = Mx^2 + Sx^2
    305 
    306 --- On cherche la pdf g(r,t) du module et de la phase
    307   x = r cos(t) ,  y = r sin(t)
    308   r=sqrt(x^2+y^2 , t=atan2(y,x)
    309   (r,t) --> (x,y): le Jacobien = r
    310 
    311   g(r,t) = r f(x,y) = r f(x) f(y)
    312          = r/(2Pi Sx Sy) exp(-(x-Mx)^2/(2 Sx^2)) exp(-(y-My)^2/(2 Sy^2))
    313 
    314 - Le cas general est complique
    315   (cf D.Pelat cours DEA "bruits et signaux" section 4.5)
    316 
    317 - Cas ou "Mx = My = 0" et "Sx = Sy = S"
    318   c'est la pdf du module et de la phase d'un nombre complexe
    319      dont les parties reelles et imaginaires sont independantes
    320      et sont distribuees selon des gaussiennes de variance S^2
    321   g(r,t) = r/(2Pi S^2) exp(-r^2/(2 S^2))
    322   La distribution de "r" est donc:
    323     g(r) = Integrate[g(r,t),{t,0,2Pi}]
    324          = r/S^2 exp(-r^2/(2 S^2))
    325   La distribution de "t" est donc:
    326     g(t) = Integrate[g(r,t),{r,0,Infinity}]
    327          = 1 / 2Pi  (distribution uniforme sur [0,2Pi[)
    328   Les variables aleatoires r,t sont independantes:
    329     g(r,t) = g(r) g(t)
    330 On a:
    331   <r>   = Integrate[r*g(r)]   = sqrt(PI/2)*S
    332   <r^2> = Integrate[r^2*g(r)] = 2*S^2
    333   <r^3> = Integrate[r^3*g(r)] = 3*sqrt(PI/2)*S^3
    334   <r^4> = Integrate[r^4*g(r)] = 8*S^4
    335 
    336 - Attention:
    337 La variable complexe "c = x+iy = r*exp(i*t)" ainsi definie verifie:
    338               <|c|^2> = <c c*> = <x^2+y^2> = <r^2> = 2 S^2
    339 Si on veut generer une variable complexe gaussienne telle que
    340      <c c*> = s^2 alors il faut prendre S = s/sqrt(2) comme argument
    341 
    342 */

  • trunk/SophyaLib/NTools/perandom.h

    r3608 r3615  
    1212#include "tvector.h"
    1313#include "histos.h"
    14 #include "srandgen.h"
    1514#include "classfunc.h"
    16 #include <complex>
    1715
    1816namespace SOPHYA { 
    … …  
    5452};
    5553
    56 
    57 /*!
    58   \ingroup NTools
    59   \brief Returns a random complex number such that real and imaginary parts are gaussians with variance sig^2
    60 */
    61 inline complex< r_8 > ComplexGaussRan(void)
    62   {return complex< r_8 >(NorRand(),NorRand());}
    63 inline complex< r_8 > ComplexGaussRan(double sig)
    64   {return complex< r_8 >(GauRnd(0.,sig),GauRnd(0.,sig));}
    65 
    66 /*!
    67   \ingroup NTools
    68   \brief Returns the module of a random complex number generated by ComplexGaussRan
    69 */
    70 inline double ModComplexGaussRan(double sig=1.)
    71   {double r=-log(1.-drand01()); return sig*sqrt(2.*r);}
    72 
    7354} // namespace SOPHYA
    7455

  • trunk/SophyaLib/NTools/smakefile

    r3608 r3615  
    3333  $(SOPHYAINCP)ndatablock.h \
    3434  $(SOPHYAINCP)utilarr.h \
     35  $(SOPHYAINCP)randinterf.h \
    3536  $(SOPHYAINCP)tmatrix.h \
    3637  $(SOPHYAINCP)tmatrix_tsnl.h \
    … …  
    6768  $(SOPHYAINCP)ndatablock.h \
    6869  $(SOPHYAINCP)utilarr.h \
     70  $(SOPHYAINCP)randinterf.h \
    6971  $(SOPHYAINCP)tmatrix_tsnl.h \
    7072  $(SOPHYAINCP)fioarr.h \
    … …  
    114116  $(SOPHYAINCP)ndatablock.h \
    115117  $(SOPHYAINCP)utilarr.h \
     118  $(SOPHYAINCP)randinterf.h \
    116119  $(SOPHYAINCP)tmatrix_tsnl.h \
    117120  $(SOPHYAINCP)tvector_tsnl.h \
    … …  
    120123 
    121124$(SOPHYAOBJP)dynccd.o: dynccd.cc $(SOPHYAINCP)sopnamsp.h \
    122   $(SOPHYAINCP)machdefs.h fmath.h perandom.h \
    123   $(SOPHYAINCP)tvector.h \
    124   $(SOPHYAINCP)tmatrix.h \
    125   $(SOPHYAINCP)machdefs.h \
    126   $(SOPHYAINCP)tarray.h \
    127   $(SOPHYAINCP)basarr.h \
    128   $(SOPHYAINCP)anydataobj.h \
     125  $(SOPHYAINCP)machdefs.h fmath.h \
     126  $(SOPHYAINCP)srandgen.h \
     127  $(SOPHYAINCP)machdefs.h \
     128  $(SOPHYAINCP)randinterf.h \
     129  $(SOPHYAINCP)anydataobj.h perandom.h \
     130  $(SOPHYAINCP)tvector.h \
     131  $(SOPHYAINCP)tmatrix.h \
     132  $(SOPHYAINCP)tarray.h \
     133  $(SOPHYAINCP)basarr.h \
    129134  $(SOPHYAINCP)mutyv.h \
    130135  $(SOPHYAINCP)timestamp.h \
    … …  
    146151  $(SOPHYAINCP)perrors.h \
    147152  $(SOPHYAINCP)fmath.h \
    148   $(SOPHYAINCP)tvector.h \
    149   $(SOPHYAINCP)srandgen.h classfunc.h cimage.h \
     153  $(SOPHYAINCP)tvector.h classfunc.h cimage.h \
    150154  $(SOPHYAINCP)tmatrix.h \
    151155  $(SOPHYAINCP)fioarr.h dynccd.h \
    … …  
    175179  $(SOPHYAINCP)ndatablock.h \
    176180  $(SOPHYAINCP)utilarr.h \
     181  $(SOPHYAINCP)randinterf.h \
    177182  $(SOPHYAINCP)tmatrix_tsnl.h \
    178183  $(SOPHYAINCP)tvector_tsnl.h generaldata.h \
    … …  
    207212  $(SOPHYAINCP)ndatablock.h \
    208213  $(SOPHYAINCP)utilarr.h \
     214  $(SOPHYAINCP)randinterf.h \
    209215  $(SOPHYAINCP)tmatrix_tsnl.h \
    210216  $(SOPHYAINCP)tvector_tsnl.h generaldata.h \
    … …  
    239245  $(SOPHYAINCP)ndatablock.h \
    240246  $(SOPHYAINCP)utilarr.h \
     247  $(SOPHYAINCP)randinterf.h \
    241248  $(SOPHYAINCP)tmatrix_tsnl.h \
    242249  $(SOPHYAINCP)tvector.h \
    … …  
    266273  $(SOPHYAINCP)ndatablock.h \
    267274  $(SOPHYAINCP)utilarr.h \
     275  $(SOPHYAINCP)randinterf.h \
    268276  $(SOPHYAINCP)tmatrix_tsnl.h \
    269277  $(SOPHYAINCP)tvector.h \
    … …  
    292300  $(SOPHYAINCP)ndatablock.h \
    293301  $(SOPHYAINCP)utilarr.h \
     302  $(SOPHYAINCP)randinterf.h \
    294303  $(SOPHYAINCP)tmatrix_tsnl.h \
    295304  $(SOPHYAINCP)tvector.h \
    … …  
    319328  $(SOPHYAINCP)ndatablock.h \
    320329  $(SOPHYAINCP)utilarr.h \
     330  $(SOPHYAINCP)randinterf.h \
    321331  $(SOPHYAINCP)tmatrix_tsnl.h \
    322332  $(SOPHYAINCP)tvector_tsnl.h \
    … …  
    346356  $(SOPHYAINCP)ndatablock.h \
    347357  $(SOPHYAINCP)utilarr.h \
     358  $(SOPHYAINCP)randinterf.h \
    348359  $(SOPHYAINCP)tmatrix_tsnl.h \
    349360  $(SOPHYAINCP)tvector_tsnl.h generaldata.h \
    … …  
    378389  $(SOPHYAINCP)ndatablock.h \
    379390  $(SOPHYAINCP)utilarr.h \
     391  $(SOPHYAINCP)randinterf.h \
    380392  $(SOPHYAINCP)tmatrix_tsnl.h \
    381393  $(SOPHYAINCP)tvector_tsnl.h generaldata.h \
    … …  
    418430  $(SOPHYAINCP)ndatablock.h \
    419431  $(SOPHYAINCP)utilarr.h \
     432  $(SOPHYAINCP)randinterf.h \
    420433  $(SOPHYAINCP)tmatrix_tsnl.h \
    421434  $(SOPHYAINCP)tvector_tsnl.h generaldata.h \
    … …  
    458471  $(SOPHYAINCP)ndatablock.h \
    459472  $(SOPHYAINCP)utilarr.h \
     473  $(SOPHYAINCP)randinterf.h \
    460474  $(SOPHYAINCP)tmatrix_tsnl.h \
    461475  $(SOPHYAINCP)tvector_tsnl.h \
    … …  
    488502  $(SOPHYAINCP)ndatablock.h \
    489503  $(SOPHYAINCP)utilarr.h \
     504  $(SOPHYAINCP)randinterf.h \
    490505  $(SOPHYAINCP)tmatrix_tsnl.h cimage.h \
    491506  $(SOPHYAINCP)fioarr.h \
    … …  
    512527  $(SOPHYAINCP)machdefs.h \
    513528  $(SOPHYAINCP)pexceptions.h \
    514   $(SOPHYAINCP)machdefs.h perandom.h \
    515   $(SOPHYAINCP)tvector.h \
    516   $(SOPHYAINCP)tmatrix.h \
    517   $(SOPHYAINCP)tarray.h \
    518   $(SOPHYAINCP)basarr.h \
    519   $(SOPHYAINCP)anydataobj.h \
     529  $(SOPHYAINCP)machdefs.h \
     530  $(SOPHYAINCP)srandgen.h \
     531  $(SOPHYAINCP)randinterf.h \
     532  $(SOPHYAINCP)anydataobj.h perandom.h \
     533  $(SOPHYAINCP)tvector.h \
     534  $(SOPHYAINCP)tmatrix.h \
     535  $(SOPHYAINCP)tarray.h \
     536  $(SOPHYAINCP)basarr.h \
    520537  $(SOPHYAINCP)mutyv.h \
    521538  $(SOPHYAINCP)timestamp.h \
    … …  
    537554  $(SOPHYAINCP)perrors.h \
    538555  $(SOPHYAINCP)fmath.h \
    539   $(SOPHYAINCP)tvector.h \
    540   $(SOPHYAINCP)srandgen.h classfunc.h pemath.h
     556  $(SOPHYAINCP)tvector.h classfunc.h
    541557        $(CXXCOMPILE) $(CXXTEMPFLG) -o $@  perandom.cc
    542558 
    … …  
    564580  $(SOPHYAINCP)ndatablock.h \
    565581  $(SOPHYAINCP)utilarr.h \
     582  $(SOPHYAINCP)randinterf.h \
    566583  $(SOPHYAINCP)tmatrix_tsnl.h \
    567584  $(SOPHYAINCP)tvector_tsnl.h \
    … …  
    594611  $(SOPHYAINCP)ndatablock.h \
    595612  $(SOPHYAINCP)utilarr.h \
     613  $(SOPHYAINCP)randinterf.h \
    596614  $(SOPHYAINCP)tmatrix_tsnl.h \
    597615  $(SOPHYAINCP)tvector_tsnl.h \
    … …  
    629647  $(SOPHYAINCP)ndatablock.h \
    630648  $(SOPHYAINCP)utilarr.h \
     649  $(SOPHYAINCP)randinterf.h \
    631650  $(SOPHYAINCP)tmatrix_tsnl.h \
    632651  $(SOPHYAINCP)tvector_tsnl.h generalfit.h \
    … …  
    680699  $(SOPHYAINCP)ndatablock.h \
    681700  $(SOPHYAINCP)utilarr.h \
     701  $(SOPHYAINCP)randinterf.h \
    682702  $(SOPHYAINCP)tmatrix.h \
    683703  $(SOPHYAINCP)tmatrix_tsnl.h \
Note: See TracChangeset for help on using the changeset viewer.