/*

  *

  *

  * Copyright (c) 1994

  * Hewlett-Packard Company

  *

  * Copyright (c) 1996,1997

  * Silicon Graphics Computer Systems, Inc.

  *

  * Copyright (c) 1997

  * Moscow Center for SPARC Technology

  *

  * Copyright (c) 1999 

  * Boris Fomitchev

  *

  * This material is provided "as is", with absolutely no warranty expressed

  * or implied. Any use is at your own risk.

  *

  * Permission to use or copy this software for any purpose is hereby granted 

  * without fee, provided the above notices are retained on all copies.

  * Permission to modify the code and to distribute modified code is granted,

  * provided the above notices are retained, and a notice that the code was

  * modified is included with the above copyright notice.

  *

  */

 #ifndef _STLP_HEAP_C

 #define _STLP_HEAP_C

 #ifndef _STLP_INTERNAL_HEAP_H

 # include <stl/_heap.h>

 #endif

 #ifndef _STLP_INTERNAL_ITERATOR_BASE_H

 # include <stl/_iterator_base.h>

 #endif

 _STLP_BEGIN_NAMESPACE

 template <class _RandomAccessIterator, class _Distance, class _Tp>

 _STLP_INLINE_LOOP

 void 

 __push_heap(_RandomAccessIterator __first,

             _Distance __holeIndex, _Distance __topIndex, _Tp __val)

 {

   _Distance __parent = (__holeIndex - 1) / 2;

   while (__holeIndex > __topIndex && *(__first + __parent) < __val) {

     *(__first + __holeIndex) = *(__first + __parent);

     __holeIndex = __parent;

     __parent = (__holeIndex - 1) / 2;

   }    

   *(__first + __holeIndex) = __val;

 }

 template <class _RandomAccessIterator, class _Distance, class _Tp>

 inline void 

 __push_heap_aux(_RandomAccessIterator __first,

                 _RandomAccessIterator __last, _Distance*, _Tp*)

 {

   __push_heap(__first, _Distance((__last - __first) - 1), _Distance(0), 

               _Tp(*(__last - 1)));

 }

 template <class _RandomAccessIterator>

 void 

 push_heap(_RandomAccessIterator __first, _RandomAccessIterator __last)

 {

   __push_heap_aux(__first, __last,

                   _STLP_DISTANCE_TYPE(__first, _RandomAccessIterator), _STLP_VALUE_TYPE(__first, _RandomAccessIterator));

 }

 template <class _RandomAccessIterator, class _Distance, class _Tp, 

           class _Compare>

 _STLP_INLINE_LOOP

 void

 __push_heap(_RandomAccessIterator __first, _Distance __holeIndex,

             _Distance __topIndex, _Tp __val, _Compare __comp)

 {

   _Distance __parent = (__holeIndex - 1) / 2;

   while (__holeIndex > __topIndex && __comp(*(__first + __parent), __val)) {

     *(__first + __holeIndex) = *(__first + __parent);

     __holeIndex = __parent;

     __parent = (__holeIndex - 1) / 2;

   }

   *(__first + __holeIndex) = __val;

 }

 template <class _RandomAccessIterator, class _Compare,

           class _Distance, class _Tp>

 inline void 

 __push_heap_aux(_RandomAccessIterator __first,

                 _RandomAccessIterator __last, _Compare __comp,

                 _Distance*, _Tp*) 

 {

   __push_heap(__first, _Distance((__last - __first) - 1), _Distance(0), 

               _Tp(*(__last - 1)), __comp);

 }

 template <class _RandomAccessIterator, class _Compare>

 void 

 push_heap(_RandomAccessIterator __first, _RandomAccessIterator __last,

           _Compare __comp) 

 {

   __push_heap_aux(__first, __last, __comp,

                   _STLP_DISTANCE_TYPE(__first, _RandomAccessIterator), _STLP_VALUE_TYPE(__first, _RandomAccessIterator));

 }

 template <class _RandomAccessIterator, class _Distance, class _Tp>

 void 

 __adjust_heap(_RandomAccessIterator __first, _Distance __holeIndex,

               _Distance __len, _Tp __val) {

   _Distance __topIndex = __holeIndex;

   _Distance __secondChild = 2 * __holeIndex + 2;

   while (__secondChild < __len) {

     if (*(__first + __secondChild) < *(__first + (__secondChild - 1)))

       __secondChild--;

     *(__first + __holeIndex) = *(__first + __secondChild);

     __holeIndex = __secondChild;

     __secondChild = 2 * (__secondChild + 1);

   }

   if (__secondChild == __len) {

     *(__first + __holeIndex) = *(__first + (__secondChild - 1));

     __holeIndex = __secondChild - 1;

   }

   __push_heap(__first, __holeIndex, __topIndex, __val);

 }

 template <class _RandomAccessIterator, class _Tp>

 inline void 

 __pop_heap_aux(_RandomAccessIterator __first, _RandomAccessIterator __last, _Tp*) {

   __pop_heap(__first, __last - 1, __last - 1, 

              _Tp(*(__last - 1)), _STLP_DISTANCE_TYPE(__first, _RandomAccessIterator));

 }

 template <class _RandomAccessIterator>

 void pop_heap(_RandomAccessIterator __first, 

 	      _RandomAccessIterator __last) {

   __pop_heap_aux(__first, __last, _STLP_VALUE_TYPE(__first, _RandomAccessIterator));

 }

 template <class _RandomAccessIterator, class _Distance,

           class _Tp, class _Compare>

 void

 __adjust_heap(_RandomAccessIterator __first, _Distance __holeIndex,

               _Distance __len, _Tp __val, _Compare __comp)

 {

   _Distance __topIndex = __holeIndex;

   _Distance __secondChild = 2 * __holeIndex + 2;

   while (__secondChild < __len) {

     if (__comp(*(__first + __secondChild), *(__first + (__secondChild - 1))))

       __secondChild--;

     *(__first + __holeIndex) = *(__first + __secondChild);

     __holeIndex = __secondChild;

     __secondChild = 2 * (__secondChild + 1);

   }

   if (__secondChild == __len) {

     *(__first + __holeIndex) = *(__first + (__secondChild - 1));

     __holeIndex = __secondChild - 1;

   }

   __push_heap(__first, __holeIndex, __topIndex, __val, __comp);

 }

 template <class _RandomAccessIterator, class _Tp, class _Compare>

 inline void 

 __pop_heap_aux(_RandomAccessIterator __first,

                _RandomAccessIterator __last, _Tp*, _Compare __comp)

 {

   __pop_heap(__first, __last - 1, __last - 1, _Tp(*(__last - 1)), __comp,

              _STLP_DISTANCE_TYPE(__first, _RandomAccessIterator));

 }

 template <class _RandomAccessIterator, class _Compare>

 void 

 pop_heap(_RandomAccessIterator __first,

          _RandomAccessIterator __last, _Compare __comp)

 {

     __pop_heap_aux(__first, __last, _STLP_VALUE_TYPE(__first, _RandomAccessIterator), __comp);

 }

 template <class _RandomAccessIterator, class _Tp, class _Distance>

 _STLP_INLINE_LOOP

 void 

 __make_heap(_RandomAccessIterator __first,

             _RandomAccessIterator __last, _Tp*, _Distance*)

 {

   if (__last - __first < 2) return;

   _Distance __len = __last - __first;

   _Distance __parent = (__len - 2)/2;

   while (true) {

     __adjust_heap(__first, __parent, __len, _Tp(*(__first + __parent)));

     if (__parent == 0) return;

     __parent--;

   }

 }

 template <class _RandomAccessIterator>

 void 

 make_heap(_RandomAccessIterator __first, _RandomAccessIterator __last)

 {

   __make_heap(__first, __last,

               _STLP_VALUE_TYPE(__first, _RandomAccessIterator), _STLP_DISTANCE_TYPE(__first, _RandomAccessIterator));

 }

 template <class _RandomAccessIterator, class _Compare,

           class _Tp, class _Distance>

 _STLP_INLINE_LOOP

 void

 __make_heap(_RandomAccessIterator __first, _RandomAccessIterator __last,

             _Compare __comp, _Tp*, _Distance*)

 {

   if (__last - __first < 2) return;

   _Distance __len = __last - __first;

   _Distance __parent = (__len - 2)/2;

   while (true) {

     __adjust_heap(__first, __parent, __len, _Tp(*(__first + __parent)),

                   __comp);

     if (__parent == 0) return;

     __parent--;

   }

 }

 template <class _RandomAccessIterator, class _Compare>

 void 

 make_heap(_RandomAccessIterator __first, 

           _RandomAccessIterator __last, _Compare __comp)

 {

   __make_heap(__first, __last, __comp,

               _STLP_VALUE_TYPE(__first, _RandomAccessIterator), _STLP_DISTANCE_TYPE(__first, _RandomAccessIterator));

 }

 _STLP_END_NAMESPACE

 #endif /*  _STLP_HEAP_C */

 // Local Variables:

 // mode:C++

 // End: