312
Z=Z+DZ;
313
UW=UW+DUW;
314
UI=PHI+UW;
315
CHI=SNYDER(PHI,DSAT(PHI),DSAT(PHI+DPHI)-DSAT(PHI));
316
UN=CHI*UW+(1-CHI)*UI;
317
}
318
319
void Trigidice::GLOBAL(){
320
InitTCALC();
321
QW1 = Y*VI + ((Y-1)*(WSAT-WF))*VB;
322
GRADUW1 = FW - QW1/KWSAT;
323
QH1 = -RKH(WSAT)*GRADTB ;
324
FSI = VI/(VI+VB);
325
GRADTHETA1 = GRADTB;
326
327
GRADPHI1 = (Y-1)* GRADUW1 - Y*H* GRADTHETA1;
328
DPHI = resolution*PHIB*sqrt(GRADPHI1);
329
DZ = DPHI/GRADPHI1;
330
resolution = DZ*DPHI;
331
UW=0;
332
UI=PHI;
333
CHI=SNYDER(PHI,DSAT(PHI),DSAT(PHI+DPHI)-DSAT(PHI));
334
UN=CHI*UW+(1-CHI)*UI;
335
THETA=-PHI/(Y*H);
336
}
337
338
339
340
int Trigidice::prec(){
341
342
if (( fabs((sum/(ZI-ZN)-WF)/WD) <= precision)
343
|| (LAYERPROB) || (GRADPROB)){
344
return 0;
345
}
346
else{
347
WF=sum/(ZI-ZN);
348
iteration++;
349
return 1;
350
}
351
}
352
353
354
int Trigidice::check(){
355
356
if (counter > LCMAX){
357
LAYERPROB=1;
358
errors();
359
iteration = -999;
360
return 0;
361
}
362
if (GRADPROB){
363
errors();
364
iteration = -888;
365
return 0;
24