Function initialize

src/geode/geometry/internal/predicates.cpp:1804–1833

        }

        void initialize()
        {
            // Taken from Jonathan Shewchuk's exactinit.
            double half;
            double check, lastcheck;
            int every_other;

            every_other = 1;
            half = 0.5;
            expansion_epsilon_ = 1.0;
            expansion_splitter_ = 1.0;
            check = 1.0;
            // Repeatedly divide `epsilon' by two until it is too small to add
            // to one without causing roundoff.  (Also check if the sum is equal
            // to the previous sum, for machines that round up instead of using
            // exact rounding.  Not that this library will work on such machines
            // anyway.
            do
            {
                lastcheck = check;
                expansion_epsilon_ *= half;
                if( every_other )
                {
                    expansion_splitter_ *= 2.0;
                }
                every_other = !every_other;
                check = 1.0 + expansion_epsilon_;
            } while( ( check != 1.0 ) && ( check != lastcheck ) );
            expansion_splitter_ += 1.0;
        }
    } // namespace PCK
} // namespace GEO