1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
|
/*
* Copyright (c) 2000-2001 Philippe Grandclement
*
* This file is part of LORENE.
*
* LORENE is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* LORENE 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 General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with LORENE; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
*/
char map_poisson_vect_C[] = "$Header: /cvsroot/Lorene/C++/Source/Map/map_poisson_vect.C,v 1.5 2014/10/13 08:53:06 j_novak Exp $" ;
/*
* $Id: map_poisson_vect.C,v 1.5 2014/10/13 08:53:06 j_novak Exp $
* $Log: map_poisson_vect.C,v $
* Revision 1.5 2014/10/13 08:53:06 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.4 2003/12/19 16:21:43 j_novak
* Shadow hunt
*
* Revision 1.3 2002/07/09 16:46:23 p_grandclement
* The Param in the case of an affine mapping is now 0x0 and not deleted
* (I wonder why it was working before)
*
* Revision 1.2 2002/05/07 07:10:45 e_gourgoulhon
* Compatibilty with xlC compiler on IBM SP2:
* suppressed the parenthesis around argument of instruction new:
* e.g. aa = new (Tbl*[nzone]) ---> aa = new Tbl*[nzone]
* result = new (Param) ---> result = new Param
*
* Revision 1.1.1.1 2001/11/20 15:19:27 e_gourgoulhon
* LORENE
*
* Revision 2.3 2000/03/07 16:53:51 eric
* niter est desormais passe en parametres.
*
* Revision 2.2 2000/02/15 10:25:49 phil
* suppression des fonctions poisson_vect et poisson_vect_oohara
*
* Revision 2.1 2000/02/09 10:01:46 phil
* ajout version oohara
*
* Revision 2.0 2000/01/21 12:59:10 phil
* *** empty log message ***
*
*
* $Header: /cvsroot/Lorene/C++/Source/Map/map_poisson_vect.C,v 1.5 2014/10/13 08:53:06 j_novak Exp $
*
*/
#include "map.h"
#include "tenseur.h"
#include "param.h"
namespace Lorene {
Param* Map_af::donne_para_poisson_vect(Param&, int) const {
return 0x0 ;
}
Param* Map_et::donne_para_poisson_vect(Param& para, int i) const {
assert ((i>=0) && (i<4)) ;
Param* result ;
result = new Param ;
result->add_int(para.get_int()) ; //nbre max iterations
result->add_double(para.get_double(0), 0) ; // relaxation
result->add_double(para.get_double(1), 1) ; // precision
if (i!=3)
result->add_cmp_mod(para.get_tenseur_mod().set(i)) ; //source au pas precedent.
else
result->add_cmp_mod(para.get_cmp_mod()) ; //la scalaire...
result->add_int_mod(para.get_int_mod()) ;
return result ;
}
}
|
