74 while ( fabs(dent(0).val_point(rr, M_PI/2., pp)) > 1.e-15 ||
75 fabs(dent(1).val_point(rr, M_PI/2., pp)) > 5.e-15) {
82 while ( dy > 1.e-15 ) {
88 rr =
sqrt( xxp*xxp + yyp*yyp ) ;
94 pp =
acos( xxp / rr ) ;
97 dval_y = dent(1).val_point(rr, M_PI/2., pp) ;
106 while( diff_dent_y > 1.e-15 ) {
109 ytmp = yyp + ss * nn * dy ;
111 rr =
sqrt( xxp*xxp + ytmp*ytmp ) ;
122 pp =
acos( xxp / rr ) ;
125 diff_dent_y = ss * dent(1).val_point(rr, M_PI/2., pp) ;
128 yyp += ss * (nn - 1) * dy ;
132 while ( dx > 1.e-15 ) {
138 rr =
sqrt( xxp*xxp + yyp*yyp ) ;
144 pp =
acos( xxp / rr ) ;
147 dval_x = dent(0).val_point(rr, M_PI/2., pp) ;
156 while( diff_dent_x > 1.e-15 ) {
159 xtmp = xxp + ss * mm * dx ;
161 rr =
sqrt( xtmp*xtmp + yyp*yyp ) ;
172 pp =
acos( xtmp / rr ) ;
175 diff_dent_x = ss * dent(0).val_point(rr, M_PI/2., pp) ;
178 xxp += ss * (mm - 1) * dx ;