Révision 20

Final-Parameters_sm_fsn_fss_fsv_freevasculature_pinnedup/pulvinus-optimize.py (revision 20)
94 94
###############################################################
95 95
###############################################################
96 96
# fitting the observed expansions        K=lambda+2*mu/3 (d=3 formula)
97
vasculature bottom can slide
97 98
###############################################################
98 99
###############################################################
99 100
Ok: Normal End
......
110 111
Results in the results file:
111 112

  
112 113
/*
113
sm;fsn;fss;fsv;sideangle;am;an;as;av;sam;san;sas;sav;tangentangle
114
sm;fsn;fss;fsv;sideangle;am;an;as;av;sam;san;sas;sav;tangentangle;cwidth
114 115
*/
115 116
-0.546671;0.972723;1.29778;0.533744;0.64322;1.67263;1.60879;1.79043;1.29461;0.964469;0.967521;0.935743;0.938534;0.784537
116 117

  
117 118

  
119
###############################################################
120
###############################################################
121
# fitting the observed expansions        K=lambda+2*mu/3 (d=3 formula)
122
vasculature bottom completely free
123
pinned up
124
###############################################################
125
###############################################################
118 126

  
127
Ok: Normal End
128
[1.67263 1.60879 1.79043 1.29462]
129
Out[2]: 
130
     fun: 7.656738339961133e-09
131
 message: 'Optimization terminated successfully.'
132
    nfev: 6005
133
     nit: 99
134
 success: True
135
       x: array([0.58126843, 0.85541866, 1.03900279, 0.52902973])
136

  
137

  
138

  
139
Results in the results file:
140

  
141
/*
142
sm;fsn;fss;fsv;sideangle;am;an;as;av;sam;san;sas;sav;tangentangle;cwidth
143
*/
144
-0.581268;0.855419;1.039;0.52903;0.116064;1.67263;1.60879;1.79043;1.29462;0.960895;0.986948;0.948169;0.948457;0.302494;1.54084
145

  
146

  
147

  
119 148
'''
Final-Parameters_sm_fsn_fss_fsv_freevasculature_pinnedup/geometry1_pinnedup.cpp (revision 20)
55 55
border vide3(t=0,1){x=-xM+t*(cwidth/2+xM);y=yM+t*(-ly-yM);label=53;};//ok
56 56

  
57 57
// ------------------- define the mesh
58

  
58
/*
59 59
plot(a1(nvertex*vthickness)
60 60
					//+a2(nvertex*2*rnect)
61 61
					+a21(nvertex*rnect)+a22(1)+a23(nvertex*rnect)
......
68 68
                     + d1(nvertex*ly)
69 69
                     );
70 70

  
71
*/
71 72

  
72 73

  
73

  
74 74
mesh Th = buildmesh (a1(nvertex*vthickness)
75 75
					//+a2(nvertex*2*rnect)
76 76
					+a21(nvertex*rnect)+a22(1)+a23(nvertex*rnect)
......
84 84
                     );
85 85

  
86 86

  
87
plot(Th);
87
//plot(Th);
88 88

  
89 89

  
90 90
/*************************************************************

Formats disponibles : Unified diff