Можно ли ускорить выполнение программы VBA (Страница 3 из 3)

<quote>
Владимир Линейцев пишет:
Тяжело ли адаптировать свои программы на Net?
</quote>
Совсем нет. Только я бы рекомендовал совсем уходить от com api, а пользовать net api.

Странно, Microsoft в своем MS Office сохраняет и развивает VBA (VBA 7), а Autodesk отказывается.

Согласно алгоритма Александра Ривилиса начал реализовывать способ двоичного поиска на VBA. По шагам пишу так:
1. нашел в сортированных массивах индексы координат расчетной точки P через двоичный поиск;
2. уменьшая или увеличивая индексы на 1 выполнил 4 проверки по коодинатам и определил нижние и верхние  границы координат в сортированных массивах по X и по Y, в центре которого находится точка P;
3. в двух последовательных циклах определил список точек лежащих внутри границ координат по оси X и по оси Y
4. т.к. получились повторы точек (одна и таже точка находилась в двух разных сорт массивах), то выполнил сортировку индексов точек массива в пункте 3
5. удалил повторения (хотя как раз они и будут находится в квадрате RxR) и точки за пределами круга R

На 25000 точек время выбора составило 8 мин 40 сек. Это пока только чистый выбор без всяких навесков. С навесками напишу попозже.
Как ускорить мой алгоритм? Где у меня не правильно, если сравнить с алгоритмом Александра Ривилиса. Или дело в медленности самого VBA?

Александр, выложите пожалуйста кусочек кода с выбором точек.

По поводу навесков и всего остального. С навесками получилось 9 мин 06 сек. Экономия времени на лицо - почти в 7 раз. Раньше этот процесс занимал 62 мин. Это здорово. Но все же...может можно еще быстрее?

Это общий код без сортировки координатных массивов. Массив TrackAll(m)(j) многослойный. Всего слоев KolTrack + 1.

    For m = 0 To KolTrack + 1
        For j = 0 To UBound()
            ' Поиск id точки mj в сортированных массивах координат
            Xsort = Replace(Format(CDbl(TrackAll(m)(j)(1)), "0.000"), ",", "")
            Ysort = Replace(Format(CDbl(TrackAll(m)(j)(2)), "0.000"), ",", "")
            mjFormat = Format(m + j / 1000000, "000.000000")
            Xmj = CDec(Xsort & mjFormat)  ' значение в сортированном массиве по X
            Ymj = CDec(Ysort & mjFormat)  ' значение в сортированном массиве по Y
            idPntX = Find2Id(xSortFull, Xmj)
            idPntY = Find2Id(ySortFull, Ymj)
            ' Определение id точек-соседей в квадрате Rn x Rn вокруг точки mj в сортированных массивах координат
            Rn = 255
            idXmin = 0: idXmax = UBound(xSortFull)
            If idPntX > 0 Then
                idXmin = idPntX - 1
                Do While Xmj - xSortFull(idXmin) < Rn * 1000000
                    idXmin = idXmin - 1
                    If idXmin < 0 Then Exit Do
                Loop
                idXmin = idXmin + 1
            End If
            If idPntX < UBound(xSortFull) Then
                idXmax = idPntX + 1
                Do While xSortFull(idXmax) - Xmj < Rn * 1000000
                    idXmax = idXmax + 1
                    If idXmax > UBound(xSortFull) Then Exit Do
                Loop
                idXmax = idXmax - 1
            End If
            idYmin = 0: idYmax = UBound(ySortFull)
            If idPntY > 0 Then
                idYmin = idPntY - 1
                Do While Ymj - ySortFull(idYmin) < Rn * 1000000
                    idYmin = idYmin - 1
                    If idYmin < 0 Then Exit Do
                Loop
                idYmin = idYmin + 1
            End If
            If idPntY < UBound(ySortFull) Then
                idYmax = idPntY + 1
                Do While ySortFull(idYmax) - Ymj < Rn * 1000000
                    idYmax = idYmax + 1
                    If idYmax > UBound(ySortFull) Then Exit Do
                Loop
                idYmax = idYmax - 1
            End If
            ' Суммирование точек в границах -Rn .. +Rn по обоим координатам в один массив
            ReDim IdPntAll(0)
            For idX = idXmin To idXmax
                idmj = xSortFull(idX) - 1000 * Fix(xSortFull(idX) / 1000)
                ReDim Preserve IdPntAll(idX - idXmin): IdPntAll(UBound(IdPntAll)) = Abs(idmj)
            Next idX
            a = idXmax - idXmin + 1
            For idY = idYmin To idYmax
                idmj = ySortFull(idY) - 1000 * Fix(ySortFull(idY) / 1000)
                ReDim Preserve IdPntAll(a + idY - idYmin): IdPntAll(UBound(IdPntAll)) = Abs(idmj)
            Next idY
            ' Быстрая сортировка точек по пути m
            SList = IdPntAll                 ' формирование массива данных быстрой сортировки SList
            Stack1 = SList                   ' присвоение размерности массиву Stack1
            Call QuickSort(0, UBound(SList)) ' Do the Quicksort - процедура быстрой сортировки (линейный массив)
            IdPntAll = SList
            ' Удаление повторов точек и точек R>Rn
            Temp = IdPntAll
            ReDim IdPntAll(0): i = 0
            For k = 0 To UBound(Temp)
                If k = 0 Then a = Temp(k) - 1 Else a = Temp(k - 1)
                If Temp(k) <> a Then
                    m1 = Int(Temp(k)): j1 = Int(1000000 * (Temp(k) - m1))
                    dX = TrackAll(m1)(j1)(1) - TrackAll(m)(j)(1)
                    dY = TrackAll(m1)(j1)(2) - TrackAll(m)(j)(2)
                    If dX * dX + dY * dY < Rn * Rn And m1 <= KolTrack Then ' проверка внутри круга
                        ReDim Preserve IdPntAll(i): IdPntAll(UBound(IdPntAll)) = Temp(k)
                        i = i + 1
                    End If
                End If
            Next k

        Next j
    Next m

Обратил внимание, что ты активно используешь ReDim Preserve внутри вложенного цикла. Думаю, что и на ней могут быть существенные потери времени. По этому поводу есть статья The ReDim Preserve Performance Trap с подзаголовком ReDim Preserve as a Performance Killer

Тут ты не прав. Ты можешь его создать сразу таким размером, чтобы помещались все точки.

Сделай его величиной равной количеству всех точек (у тебя в данном случае кажется их 25000). Когда подсчитаешь сколько точек действительно в него попадает и внесешь их в массив - уменьшишь размер массива до этого значения (хотя можно и не уменьшать).

Ты меня совсем запутал. Допустим у тебя всего точек N (это все, которые в принципе могут быть соседями). Создаешь массив iArr размером N. Затем вычисляешь реальных соседей и заполняешь последовательно массив iArr от 0 до M-1 (M - количество реальных соседей). После этого ты можешь (при желании) уменьшить массив до M. Т.е.: ReDim iArr(0 to M-1). Как-то так.