source: trunk/examples/berrada/MFLBC_RD/bet.h @ 385

Last change on this file since 385 was 385, checked in by lnalod, 14 years ago

Add of the Berrada applications MFLBC and MFLBC_RD. The second is also useful for the test of the parallelization algorithm.
  • Property svn:eol-style set to native
File size: 1.3 KB
Line 
1
2/************************ Berrada Mohamed 2-2005 ***************************
3                 module classe bet.h  -  description
4  rappel :
5  .YSi   correspond ? la i?me sortie
6  .YJiIn correspond ? la d?riv?e de la i?me sortie par rapport ni?me entr?e
7 ***************************************************************************/
8
9  //===========================================================================
10  //                          methodes forward
11  forward (YREAL xdiacr,YREAL xdiacc,YREAL xgagr,YREAL xgagc)
12  {
13        complex<double>         xdiac (xdiacr,xdiacc);
14        complex<double>         xgag (xgagr,xgagc);
15
16        complex<double> var1(0.,0.);
17
18//
19        if(Yj != 0){
20                var1=xdiac;
21                if(Yj != 1) var1=xdiac-zig*xgag;
22                }
23//
24        YS1 = real(var1);
25        YS2 = imag(var1);
26
27
28  }
29
30  //===========================================================================
31  //                         methodes  backward
32  backward (YREAL xdiacr,YREAL xdiacc,YREAL xgagr,YREAL xgagc)
33  {
34        complex<double>         xdiac (xdiacr,xdiacc);
35        complex<double>         xgag (xgagr,xgagc);
36
37        if(Yj != 0){
38                Ycx_jac( J, 1, 1);
39                if(Yj != 1) {
40                        Ycx_jac( -zig*J, 1, 2);
41                        }
42                }
43
44//
45   }
46
47  //===========================================================================
48//************************** FIN DU MODULE bet ********************************
49
50
51
52
53
54
55
56
57
58
59
Note: See TracBrowser for help on using the repository browser.