long double ldRootSearch( long double ldK ) { long double ldMax = 1500.0e0L; long double ldMin = 268.0e0L; long dougle ldGuess = ( ldMax + ldMin ) * 0.5e0L; int iBailout = 0; while ( fabsl( ldK - powl( ldGuess, ldGuess ) ) / ldK >= 1.0e-17 ) { if ( powl( ldGuess, ldGuess ) > ldK ) ldMax = ldGuess; else ldMin = ldGuess; ldGuess = ( ldMax + ldMin ) * 0.5e0L; iBailout++; } }