| 390 | } |
| 391 | |
| 392 | double ChargeF(int i1,int i2,int i3,Complex u[3], double const & eps) |
| 393 | { |
| 394 | double ch = |
| 395 | intdphase(i1<i2,u[0],u[1],eps) |
| 396 | +intdphase(i2<i3,u[1],u[2],eps) |
| 397 | +intdphase(i3<i1,u[2],u[0],eps); |
| 398 | double chn = round(ch); |
| 399 | if(verbosity>99) cout << u[0]<< " " << u[1] << " " << u[2] << " ch = " << ch << endl; |
| 400 | |
| 401 | return ch ; |
| 402 | } |
| 403 | long uZero(pf3c const & fu, pf3r const & fuc, double const &eps) |
| 404 | { |
| 405 | const double epscharge = 1e-3; // |