App Store: Кубатурник леса ГОСТ 2708-75
Описание
ФУНКЦИИ ПРИЛОЖЕНИЯ:
► объем круглых бревен по ГОСТ 2708-75
► объем оцилиндрованных бревен
► выбор любого диаметра и любой длины
► удобный ввод количества бревен (не нужно нажимать 40 раз)
► подсчет количества и объема бревен по каждому диаметру
► подсчет суммарного объема, количества и цены для каждой группы диаметров
► создание, редактирование, удаление и экспорт расчетов
► для каждого расчета: выбор валюты, клиента, даты отгрузки
► настройка цен для разных групп диаметров
► удобный выбор диаметров на панели: можно настроить выбор только четных или нечётных диаметров, а также диаметров с заданными номерами
► учет брака и балансов (по необходимости)
► выбор округления для оцилиндрованных бревен
► сохранение расчетов в приложении
► поиск и сортировка расчетов по разным параметрам: по типу, имени, клиенту, дате отгрузки или создания, объему
► экспорт данных в Excel или Google Таблицы
► скриншот итогов
► отправка файла или скриншота по почте, whatsApp или другим мессенджерам
► сохранение файла на устройстве или в iCloud
► передача файла через AirDrop
► добавление, изменение и удаление клиентов
► удобная помощь и оперативная техподдержка по почте
► NEW!!! Создание шаблонов для расчетов
► NEW!!! Удаление всех записей в расчете
► NEW!!! Замена длин для всех бревен в расчете
► NEW!!! Выбор сортировки бревен в Итогах: сортировка по Диаметру или по Длине бревен
► NEW!!! Выбор отображения бревен в Расчете: показывать сначала Диаметр, а потом Длину, либо наоборот: сначала Длину, а потом Диаметр
► NEW!!! Объединение нескольких расчетов в один расчет
► NEW!!! Выбор и Удаление нескольких расчетов
► NEW!!! Задать для всех новых расчетов диаметры и цены по умолчанию.
Объем по ГОСТ 2708-75:
Берется из стандарта «ГОСТ 2708-75. ЛЕСОМАТЕРИАЛЫ КРУГЛЫЕ. ТАБЛИЦЫ ОБЪЕМОВ». В приложении представлены объемы бревен:
— длиной от 1,0 до 9,5 м и толщиной от 3 до 120 см
— длиной от 0,5 до 0,9 м и толщиной от 6 до 15 см
— длиной от 10,0 до 13,5 м и толщиной от 8 до 38 см
Объем оцилиндрованных бревен:
Считается по формуле πr² • H, где:
π — 3,14
r² — радиус оцилиндр-го бревна в квадрате
r = D/2, где D — диаметр бревна
Н — длина оцилиндрованного бревна
Версия 4.2.0
ВНИМАНИЕ! Из-за санкций App Store больше не принимает платежи по карте. Теперь совершить покупку можно двумя способами.
► Способ 1: Настройте оплату через Мобильный телефон. Для этого откройте приложение Настройки -> Нажмите на Apple ID (Ваше имя) -> Оплата и доставка -> Добавить способ Оплаты -> Мобильный телефон. Убедитесь, что на счете Вашего телефона достаточно средств для покупки.
► Способ 2: Пополните Ваш счет Apple ID. Для этого откройте App Store → Профиль (человечек в правом верхнем углу) → Пополнить счет Apple ID. Если у Вас возникнут вопросы, напишите нам на почту [email protected], мы поможем!
Что нового в этой версии:
1) Добавлена функция «Задать диаметры и цены по умолчанию для всех новых расчетов». Для включения данной функции нажмите вкладку «Настройки приложения» в правом нижнем углу экрана, затем в Настройке «Цена по умолчанию за м3» введите нужные диаметры и цену для них.
2) Теперь в расчете данные диаметра, длины и количества выравниваются по центру.
Оценки и отзывы
Оценок: 34
Дмитрий
Отличный кубатурник, пользовался на андроиде, но таскать с собой второй телефон только из-за кубатурника надоело. Купил на айфон и не жалею)
Уровень на 5+
Кубатурник бревна отличнейший все ясно быстро и удобно ! Добавьте в нем и кубатуру Доски ! Цены вам не будет!)))
Большое спасибо за Ваш отзыв! Очень рады, что приложение удобное. Мы выпустили новое приложение «Кубатурник доски леса Pro ГОСТ», в нем очень много функций по расчету кубатуры, площади, погонных метров доски. Можно выбирать любые единицы измерения, удобно вводить цену, ширину, толщину, длину по умолчанию, также можно считать пильные и зачетные размеры по заданному коэффициенту потери, заданным допускам на усушку или по ГОСТ.
Будем рады, если Вы попробуете наше новое приложение и напишите, все ли Вас устраивает и нужны ли какие-либо дополнительные функции:)
Не доволен
Постоянно логает, уже в сохранённых списках меняет количество и диаметры самостоятельно.
1 июня 2020:
UPD: Добрый день! Версия 1.3.0. со всеми исправлениями уже доступна в App Store)
Мы будем очень признательны, если Вы напишите нам сюда или на почту (в нашем приложении Вкладка помощь -> кнопка Задать вопрос), если после обновления у Вас все же будут возникать проблемы в работе приложения.
Если же после обновления приложение будет работать корректно и без лагов, мы будем очень признательны, если Вы поставите нам более высокую оценку:)
Большое спасибо за внимание!
С уважением, команда разработчиков приложения «Кубатурник леса».
——
30 мая 2020:
Добрый день! Большое спасибо за Ваш ценный отзыв.
Мы исправили ошибку «уже в сохранённых списках меняет количество и диаметры самостоятельно». Данная ошибка возникала, когда Вы меняли количество или диаметр бревна через текстовое поле в одном расчете, а потом открывали другой расчет, и в нем также отображались измененные данные. Теперь эта проблема устранена и все работает исправно, пересохранений данных нет.
Также мы исправили вылетание программы в некоторых случаях: при изменении диаметра, количества и длины бревен в расчете, а также при изменении диаметров и цен в Настройках расчета. Теперь приложение работает стабильно, без лагов и вылетов.
Разработчик Alisa Potapova указал, что в соответствии с политикой конфиденциальности приложения данные могут обрабатываться так, как описано ниже. Подробные сведения доступны в политике конфиденциальности разработчика.
Сбор данных не ведется
Разработчик не ведет сбор данных в этом приложении.
Конфиденциальные данные могут использоваться по-разному в зависимости от вашего возраста, задействованных функций или других факторов. Подробнее
Информация
- Провайдер
- Alisa Potapova
- Размер
- 26 МБ
- Категория
Утилиты
- Возраст
- 4+
- Copyright
- © Alisa Potapova 2022
- Цена
- 179,00 ₽
Поддержка приложения
Политика конфиденциальности
Поддерживается
Другие приложения этого разработчика
Вам может понравиться
Таблица кубатурник пиломатериала, сколько штук бруса и доски в кубе
Таблица кубатурник доски, сколько досок в кубе
Таблица пиломатериала — кубатурник — сколько штук доски обрезной в кубе.
Размер, мм | Кол-во штук в одном м3 | Кол-во погонных метров в одном м3 | Объём одной доски, м3 | Площадь одной доски, кв.м. | Вес одной доски, кг (влажн. 20%) |
---|---|---|---|---|---|
20х100х6000 | 83,3 | 500 | 0,012 | 0,6 | 9 |
20х100х3000 | 167,6 | 500 | 0,006 | 0,3 | 4,5 |
20х150х6000 | 55,6 | 333,3 | 0,018 | 0,9 | 13,5 |
20х150х3000 | 111,1 | 333,3 | 0,009 | 0,45 | 6,75 |
22х100х6000 | 75,8 | 454,5 | 0,0132 | 0,6 | 9,9 |
22х100х3000 | 151,5 | 454,5 | 0,0066 | 0,3 | 4,95 |
22х125х6000 | 60,6 | 363,6 | 0,0165 | 0,75 | 12,375 |
22х125х3000 | 121,2 | 363,6 | 0,00825 | 0,0375 | 6,19 |
22х150х6000 | 50,5 | 303 | 0,0198 | 0,09 | 14,85 |
22х150х3000 | 101 | 303 | 0,0099 | 0,045 | 7,425 |
22х175х6000 | 43,3 | 259,7 | 0,0231 | 1,05 | 17,325 |
22х175х3000 | 86,6 | 259,7 | 0,01155 | 0,525 | 8,66 |
22х200х6000 | 37,9 | 227,3 | 0,0264 | 1,2 | 19,8 |
22х200х3000 | 75,8 | 227,3 | 0,0132 | 0,6 | 9,9 |
22х225х6000 | 33,7 | 202 | 0,0297 | 1,35 | 22,275 |
22х225х3000 | 37,3 | 202 | 0,01485 | 0,675 | 11,138 |
22х250х6000 | 30,3 | 181,8 | 0,033 | 1,5 | 24,75 |
22х250х3000 | 60,6 | 181,8 | 0,0165 | 0,75 | 12,375 |
25х100х6000 | 66,7 | 400 | 0,015 | 0,6 | 11,25 |
25х100х3000 | 133,3 | 400 | 0,0075 | 0,3 | 5,625 |
25х100х2000 | 200 | 400 | 0,005 | 0,2 | 3,75 |
25х125х6000 | 53,3 | 320 | 0,01875 | 0,75 | 14,06 |
25х125х3000 | 106,7 | 320 | 0,009375 | 0,0375 | 7,03 |
25х150х6000 | 44,4 | 266,7 | 0,0225 | 0,9 | 16,875 |
25х150х3000 | 88,9 | 266,7 | 0,01125 | 0,45 | 8,44 |
25х150х2000 | 133,3 | 266,7 | 0,0075 | 0,3 | 5,625 |
25х175х6000 | 38,1 | 228,6 | 0,02625 | 1,05 | 19,69 |
25х175х3000 | 76,2 | 228,6 | 0,012125 | 0,525 | 9,094 |
25х200х6000 | 33,3 | 200 | 0,03 | 1,2 | 22,5 |
25х200х3000 | 66,7 | 200 | 0,015 | 0,6 | 11,25 |
25х225х6000 | 29,6 | 177,8 | 0,03375 | 1,35 | 25,31 |
25х225х3000 | 59,3 | 177,8 | 0,016875 | 0,675 | 12,656 |
25х250х6000 | 26,7 | 160 | 0,0375 | 1,5 | 28,125 |
25х250х3000 | 53,3 | 160 | 0,01875 | 0,75 | 14,06 |
32х100х6000 | 52,1 | 312,5 | 0,0192 | 0,6 | 14,4 |
32х100х3000 | 104,2 | 312,5 | 0,0096 | 0,3 | 7,2 |
32х125х6000 | 41,7 | 250 | 0,024 | 0,75 | 18 |
32х125х3000 | 83,3 | 250 | 0,012 | 0,0375 | 9 |
32х150х6000 | 34,7 | 208,3 | 0,0288 | 0,9 | 21,6 |
32х150х3000 | 69,4 | 208,3 | 0,0144 | 0,45 | 10,8 |
32х175х6000 | 29,8 | 178,6 | 0,0336 | 1,05 | 25,2 |
32х175х3000 | 59,5 | 178,6 | 0,0168 | 0,525 | 12,6 |
32х200х6000 | 26 | 156,3 | 0,0384 | 1,2 | 28,8 |
32х200х3000 | 52,1 | 156,3 | 0,0192 | 0,6 | 14,4 |
32х225х6000 | 23,1 | 138,9 | 0,0432 | 1,35 | 32,4 |
32х225х3000 | 46,3 | 138,9 | 0,0216 | 0,675 | 16,2 |
32х250х6000 | 20,8 | 125 | 0,048 | 1,5 | 36 |
32х250х3000 | 41,7 | 125 | 0,024 | 0,75 | 18 |
40х100х6000 | 41,7 | 250 | 0,024 | 0,6 | 18 |
40х100х3000 | 83,3 | 250 | 0,012 | 0,3 | 9 |
40х125х6000 | 33,3 | 200 | 0,03 | 0,75 | 22,5 |
40х125х3000 | 66,7 | 200 | 0,015 | 0,0375 | 11,25 |
40х150х6000 | 27,8 | 166,7 | 0,036 | 0,9 | 27 |
40х150х3000 | 55,6 | 166,7 | 0,018 | 0,45 | 13,5 |
40х175х6000 | 23,8 | 142,9 | 0,042 | 1,05 | 31,5 |
40х175х3000 | 47,6 | 142,9 | 0,021 | 0,525 | 15,75 |
40х200х6000 | 20,8 | 125 | 0,048 | 1,2 | 36 |
40х200х3000 | 41,7 | 125 | 0,024 | 0,6 | 18 |
40х225х6000 | 18,5 | 111,1 | 0,054 | 1,35 | 40,5 |
40х225х3000 | 37 | 111,1 | 0,027 | 0,675 | 20,25 |
40х250х6000 | 16,7 | 100 | 0,06 | 1,5 | 45 |
40х250х3000 | 33,3 | 100 | 0,03 | 0,75 | 22,5 |
50х100х6000 | 33,3 | 200 | 0,03 | 0,6 | 22,5 |
50х100х3000 | 66,7 | 200 | 0,015 | 0,3 | 11,25 |
50х125х6000 | 26,7 | 160 | 0,0375 | 0,75 | 28,125 |
50х125х3000 | 53,3 | 160 | 0,01875 | 0,0375 | 14,06 |
50х150х6000 | 22,2 | 133,3 | 0,045 | 0,9 | 33,75 |
50х150х3000 | 44,4 | 133,3 | 0,0225 | 0,45 | 16,875 |
50х175х6000 | 19 | 114,3 | 0,0525 | 1,05 | 39,375 |
50х175х3000 | 38,1 | 114,3 | 0,02625 | 0,525 | 19,688 |
50х200х6000 | 16,7 | 100 | 0,06 | 1,2 | 45 |
50х200х3000 | 33,3 | 100 | 0,03 | 0,6 | 22,5 |
50х225х6000 | 14,8 | 88,9 | 0,0675 | 1,35 | 50,625 |
50х225х3000 | 29,6 | 88,9 | 0,03375 | 0,675 | 25,31 |
50х250х6000 | 13,3 | 80 | 0,075 | 1,5 | 56,25 |
50х250х3000 | 26,7 | 80 | 0,0375 | 0,75 | 28,125 |
60х125х6000 | 22,2 | 133,3 | 0,045 | 0,75 | 33,75 |
60х125х3000 | 44,4 | 133,3 | 0,0225 | 0,0375 | 28,125 |
60х150х6000 | 18,5 | 111,1 | 0,054 | 0,9 | 40,5 |
60х150х3000 | 37 | 111,1 | 0,027 | 0,45 | 20,25 |
60х175х6000 | 15,9 | 95,2 | 0,063 | 1,05 | 47,25 |
60х200х6000 | 13,9 | 83,3 | 0,072 | 1,2 | 54 |
60х225х6000 | 12,3 | 74,1 | 0,081 | 1,35 | 60,75 |
60х250х6000 | 11,1 | 66,7 | 0,09 | 1,5 | 67,5 |
60х250х3000 | 22,2 | 66,7 | 0,045 | 0,75 | 33,75 |
75х175х6000 | 12,7 | 76,2 | 0,07875 | 1,05 | 59,06 |
75х175х3000 | 25,4 | 76,2 | 0,0394 | 0,525 | 29,55 |
75х200х6000 | 11,1 | 66,7 | 0,09 | 1,2 | 67,5 |
75х200х3000 | 22,2 | 66,7 | 0,045 | 0,6 | 33,75 |
75х225х6000 | 9,9 | 59,3 | 0,101 | 1,35 | 75,75 |
75х225х3000 | 19,7 | 59,3 | 0,051 | 0,675 | 38,25 |
75х250х6000 | 8,9 | 53,3 | 0,1123 | 1,5 | 84,225 |
Таблица кубатурник бруса, сколько бруса в кубе
Таблица расчета кубатуры бруса — кубатурник — сколько штук бруса в кубе.
Размер, мм | Кол-во штук в одном м3 | Кол-во погонных метров в одном м3 | Объём одной штуки, м3 | Вес одной штуки, кг (влажн. 20%) | 50х50х6000 | 66,67 | 400 | 0,015 | 11,25 | 50х50х3000 | 133,33 | 400 | 0,0075 | 5,625 |
---|---|---|---|---|
50х100х6000 | 33,33 | 200 | 0,03 | 22,5 |
50х100х3000 | 66,67 | 200 | 0,015 | 11,25 |
60х60х6000 | 46,3 | 277,78 | 0,0216 | 16,2 |
60х60х3000 | 92,6 | 277,78 | 0,0108 | 8,1 |
60х100х6000 | 27,78 | 166,67 | 0,036 | 27 |
60х100х3000 | 55,55 | 166,67 | 0,018 | 13,5 |
75х75х6000 | 29,63 | 177,78 | 0,03375 | 25,31 |
75х75х3000 | 59,26 | 177,78 | 0,0169 | 12,675 |
75х100х6000 | 22,22 | 133,33 | 0,045 | 33,75 |
75х100х3000 | 44,44 | 133,33 | 0,0225 | 16,875 |
75х150х6000 | 14,8 | 88,89 | 0,0675 | 50,625 |
100х100х6000 | 16,67 | 100 | 0,06 | 45 |
100х100х3000 | 33,33 | 100 | 0,03 | 22,5 |
100х150х6000 | 11,11 | 66,67 | 0,09 | 67,5 |
100х200х6000 | 8,33 | 50 | 0,12 | 90 |
100х250х6000 | 6,67 | 40 | 0,15 | 112,5 |
125х125х6000 | 10,67 | 64 | 0,09375 | 70,31 |
150х150х6000 | 7,41 | 44,44 | 0,135 | 101,25 |
200х200х6000 | 4,17 | 25 | 0,24 | 180 |
200х250х6000 | 3,33 | 20 | 0,3 | 225 |
250х250х6000 | 2,67 | 16 | 0,375 | 281,25 |
ЗАДАТЬ ВОПРОС
Остались вопросы? Напишите нам, мы ответим Вам в ближайшее время!
Библиотека Кубы — arXiv Vanity
Т. Хан
Институт физики им. Макса Планка
Föhringer Ring 6, D–80805 Мюнхен, Германия
Abstract
Концепции и реализация библиотеки Cuba для многомерных
разъясняется численное интегрирование.
1 Обзор
Куба [2] — это библиотека для многомерных числовых
интеграция. Он имеет четыре алгоритма интеграции с интерфейсами
для Фортрана, C/C++ и Mathematica. Все четыре могут интегрировать вектор
подынтегральные выражения и их вызов очень похожи, чтобы упростить
переключать подпрограммы для сравнения. Обобщаются основные характеристики
ниже.
Обычный:
Вегас
Доступны основные методы интеграции:
∙ проба Соболя (квази Монте-Карло)
∙ Образец Mersenne Twister (псевдо Монте-Карло)
Снижение отклонения:
выборка по важности
Обычный:
учтивый
Доступны основные методы интеграции:
∙ проба Соболя (квази Монте-Карло)
∙ Образец Mersenne Twister (псевдо Монте-Карло)
Снижение отклонения:
выборка по важности в сочетании с глобально адаптивным подразделением
Обычный:
Дивонн
Доступны основные методы интеграции:
∙ Проба Коробова (решеточный метод)
∙ проба Соболя (квази Монте-Карло)
∙ Образец Mersenne Twister (псевдо Монте-Карло)
∙ кубатурные правила (детерминированный метод)
Снижение отклонения:
стратифицированная выборка с использованием методов
численная оптимизация
Обычный:
Куре
Доступен базовый метод интеграции:
∙ кубатурные правила (детерминированный метод)
Снижение отклонения:
глобально адаптивное подразделение
Прежде чем объяснять модные словечки, встречающиеся в этом списке, возможно,
стоит ответить на два часто задаваемых вопроса.
Во-первых, численное интегрирование быстро усложняется с
увеличивающееся измерение, независимо от того, сколько трюков встроено в
интегратор. Чтобы проиллюстрировать это, представьте, что вы вычисляете объем
д-тусклый. сферы Sd путем интегрирования ее характеристической функции χ=θ(1−∥x∥2) внутри окружающего гиперкуба Cd=[−1,1]d. Таким образом, с точки зрения Монте-Карло отношение r=VolSd/VolCd можно рассматривать как вероятность того, что универсальная
интегратор найдет шар вообще. Цифры четко видны
то, что часто называют «проклятием размерности»:
д251050р.785.164.00251,5×10−28 | (1) |
Вторая, Куба (и, если уж на то пошло, самая многомерная
интеграторы) могут вычислять только интегралы Римана вида
Если:=∫10ddxf(→x). | (2) |
Большинство вопросов касаются границ, хотя это просто
преобразовать конечную область интегрирования в единичный гиперкуб:
∫b1a1⋯∫bdadddxf(→x)=∫10ddyJf(→x), | (3) | ||
J=d∏i=1(bi-ai),xi=ai+(bi-ai)yi. |
2 концепции
2.1 Детерминистический анализ по сравнению с методом Монте-Карло 90 136
Куба содержит как детерминированные методы, так и методы интегрирования Монте-Карло.
Детерминированный подход основан на кубатурных правилах ,
If≈Cnf:=n∑i=1wif(→xi) | (4) |
со специально подобранными узлами →xi и весами wi. Оценка Монте-Карло, хотя и очень похожая по форме, концептуально сильно отличается, поскольку эта формула обозначает Начиная с ортогонального базиса функций {b1,…,bm} – для nd+n неизвестных →xi и wi. Они представляют грозный, После того, как будет доступна оценка ошибки для интеграла, глобальная Интегрируем всю область: Итот±Этот. , в то время как Etot>max(εrelItot,εabs) Найдите область r с наибольшей ошибкой. Разделить пополам (или иным образом разрезать) r. Интегрировать каждый подрегион r отдельно. Itot=∑Ii, конец, а Здесь уместно сделать замечание о двух точностях, εrel и Выборка по важности вводит весовую функцию в интеграл: с двумя требованиями: Работа по стратифицированной выборке по субрегионам выборки. Рассмотрим в общей сложности , тогда как в первом случае ее можно записать как Даже в этом простом примере последняя дисперсия в лучшем случае равна Методы квази-Монте-Карло основаны на неравенстве Коксмы–Главки. Кроме выбора другого подынтегрального выражения мало что можно сделать где ν(r) подсчитывает →xi, попадающие в r. Слово Куба предлагает на выбор квазислучайные последовательности Соболя [5] или Рисунок 1: Сравнение последовательностей. Решетчатые методы требуют периодического подынтегрального выражения, обычно получаемого →z выбирается (путем обширных компьютерных поисков) для выбивания как где L⊥={→k∈Zd:→k⋅→z=0(modn)} — обратная решетка. Vegas — это классический алгоритм Лепажа Монте-Карло [8] . Оно использует В реализации Кубы Вегас может запоминать свою сетку для последующего Внимательный — гибрид Вегаса и Скупца [9] , Монте-Карло. Алгоритм работы следующий: пока не будет достигнута требуемая точность, Сетка Вегаса хранится между подразделениями, т. е. регионом, который является Divonne — это значительно расширенная версия алгоритма CERNlib D151. Фаза 1: Разделение Регион интеграции разбит на субрегионы (приблизительно) Минимум и максимум каждого субрегиона ищутся с использованием методов из \SetScale.7 Фаза 2: Отбор проб Субрегионы, определенные на этапе 1, берутся независимо друг от друга с помощью Этап 3: уточнение Регионы, результаты которых на этапах 1 и 2 не совпадают в пределах их В реализацию Cuba добавлены две важные функции: Пользователь может указывать экстремумы для сложных подынтегральных выражений. Для подынтегральных выражений, выборка которых невозможна слишком близко к границе, Cuhre — детерминированный алгоритм. Он использует кубатуру Генца – Малика. Cuhre был повторно реализован на Кубе в основном для последовательного Сбалансированное сравнение алгоритмов интеграции почти В этом контексте следует отметить, что статистическая ошибка С учетом этих предостережений на следующем графике сравнивается производительность Отдельного упоминания заслуживает интерфейс Mathematica, так как он не \SetScale.75 После загрузки соответствующего исполняемого файла MathLink, например. с Куба включает «универсальный интерфейс», который еще больше упрощает Пользователь просто должен выбрать метод = 1,2,3,4 для переключения между Средство просмотра разделов Cuba отображает раздел, занятый интеграцией Уровень детализации 3 должен быть выбран в подпрограмме интеграции и , который затем отобразит 1–2 плоскость разбиения.
Оценка ошибки выполняется, например. по нулевым правилам Nm (m Если≈Mnf:=1nn∑i=1f(→xi), (5)
статистическое среднее по независимым и одинаково распределенным
случайные выборки →xi. В этом случае стандартное отклонение
дает вероятностную оценку ошибки интегрирования: σ(Mnf)=√Mnf2−M2nf. (6) 2.2 Построение правил кубатуры
обычно мономы — с которыми можно (надеюсь) аппроксимировать большинство f
достаточно хорошо, требуется, чтобы каждое bi интегрировалось точно по формуле
Cn: Ibi\tiny!=Cnbi. Это m уравнений моментов n∑k=1wkbi(→xk)=∫10ddxbi(→x) (7)
в общем случае нелинейная система уравнений. Дополнительные предположения, например. симметрии обычно необходимы для решения этой системы. на Кубе работают
правила Генца – Малика [3] , построенные из симметричного
мономиальная основа. 2.3 Глобально адаптивное подразделение
адаптивность легко реализовать:
Eобщ=√∑E2i.
εабс. Наивно навязывают относительную точность:
результат должен быть точным, скажем, до одной тысячной, т. е. εотн=10−3. Однако для интегральных значений, стремящихся к нулю,
эта цель становится все труднее и труднее достичь, и чтобы не тратить
чрезмерное количество времени в таких случаях, абсолютная точность
можно задать εабс, где обычно εабс≪εотн. 2.4 Выборка по важности
Если=∫10ddxw(→x)f(→x)w(→x), (8) w(→x)>0,Iw=1,
а) нужно иметь возможность выбирать из распределения w(→x),
б) f/w должен быть «гладким» в том смысле, что
σw(f/w)<σ(f), например w и f должны иметь
такая же пиковая структура.
Известно, что идеальным выбором является w(→x)=|f(→x)|/If, который
имеет σw(f/w)=0, но малопригоден, так как требует априорного
знание интегрального значения. 2.5 Стратифицированная выборка
выборка n точек в области r=ra+rb по сравнению с выборкой n/2 точек
в ра и н/2 в рб. В последнем случае дисперсия равна. σ2fn=σ2af+σ2bf2n+(Iaf−Ibf)24n. (10)
прежний, и только если интегральные значения идентичны. Оптимальный
можно показать, что уменьшение дисперсии происходит при na/nb=σaf/σbf [4] . Таким образом, рецепт состоит в том, чтобы разделить
область интегрирования на части с одинаковой дисперсией, а затем произвести выборку всех
части с одинаковым количеством очков. 2.6 Методы квази-Монте-Карло 90–136
который устанавливает верхнюю границу ошибки формулы интегрирования
Qnf=1n∑ni=1f(→xi), |Qnf−If|⩽V(f)D∗(→x1,…,→xn). (11)
о V(f), «вариации в смысле Харди и Краузе».
невязка D∗ последовательности →x1,…,→xn
определяется как D*=supr∈[0,1]d∣∣∣ν(r)n−Volr∣∣∣, (12)
«равнораспределенный» действительно обычно означает, что ν(r) пропорциональна
к Волру. Квазислучайные последовательности могут быть построены с помощью
существенно более низкое несоответствие, чем (псевдо-)случайные числа. Монте
Алгоритм Карло, основанный на этих последовательностях, обычно достигает сходимости.
ставки O(logd−1n/n), а не обычные O(1/√n).
псевдослучайные последовательности Mersenne Twister [6] для всех
Алгоритмы Монте-Карло. Рисунок 1 показывает, что
квазислучайные числа покрывают плоскость гораздо более однородно, чем
псевдослучайные числа. 2.7 Методы решетки
применение преобразование периодизации (Кубинская Дивонна использует
х→|2x−1|). Отбор проб производится на интегрированной решетке L
натянутый на тщательно подобранный целочисленный вектор →z: Lnf=1nn−1∑i=0f({in→z}), (13) {x}=дробная часть x.
как можно больше «брегговских отражений» низкого порядка в члене ошибки
(см. , например, [7] ): Lnf-If =∑→k∈Zd~f(→k)Lne2πi→k⋅→x−~f(→0) =∑→k∈L⊥,→k≠→0~f(→k), (14) 3 Реализация
3.1 Вегас
выборка по важности для уменьшения дисперсии, для которой итеративно
строит кусочно-постоянную весовую функцию, представленную на
прямоугольная сетка. Каждая итерация состоит из шага выборки, за которым следует
путем измельчения сетки.
вызовы, и он может периодически сохранять свое внутреннее состояние, так что
расчет может быть возобновлен, например. после аварии. 3.2 Вежливый
алгоритм, сочетающий выборку важности в стиле Вегаса с глобальным
адаптивное подразделение.
разделить пополам область с наибольшей ошибкой вдоль оси, в которой
флуктуации подынтегральной функции уменьшаются больше всего. Пропорционально количеству
новые образцы в каждой половине для его колебания.
результат n−1 подразделений имел n итераций Вегаса, выполненных на
это. С другой стороны, Suave несколько требователен к памяти, так как ему нужно
сохранить образцы для последующего использования. 3.3 Дивонн
[10] . По сути, это алгоритм Монте-Карло, но
Кубатурные правила также встроены для сравнения. Снижение дисперсии
стратифицированная выборка, которой помогают методы численного
оптимизация. Divonne имеет трехфазный алгоритм:
равный спред s, определяемый как s(r)=Volr2(max→x∈rf(→x)−min→x∈rf(→x)). (15)
численная оптимизация (по сути, квазиньютоновский поиск). Затем,
«разделители» перемещаются (см. рисунок), чтобы найти оптимальное
расщепление. Эта последняя процедура может быть искусно переведена в
решение линейной системы и, следовательно, довольно быстро (подробности см.
[10] ).
\CBox(5,5)(95,95)СинийПастельСиний
\SetWidth2
\SetColorRed
\Линия(0,20)(100,20)
\Линия(0,65)(100,65)
\Линия(30,0)(30,100)
\Линия(85,0)(85,100)
\SetWidth0
\SetColorBlack
\ДлиннаяСтрелка(1,18)(1,14)
\ДлиннаяСтрелка(1,22)(1,26)
\ДлиннаяСтрелка(1,63)(1,59)
\ДлиннаяСтрелка(1,67)(1,71)
\ДлиннаяСтрелка(28,1)(24,1)
\ДлиннаяСтрелка(32,1)(36,1)
\ДлиннаяСтрелка(83,1)(79,1)
\ДлиннаяСтрелка(87,1)(91,1)
одинаковое количество баллов у каждого. Последний экстраполируется из результатов
этапа 1.
ошибки подразделяются или снова отбираются. Этот этап является дополнением к
исходный алгоритм, так как было обнаружено, что достаточно часто ошибка
оценка, или даже интегральное значение, были отключены, потому что характеристики
подынтегральная функция не была найдена в фазе 1.
может быть задано «безопасное расстояние», в пределах которого значения будут
экстраполируется из двух внутренних точек. 3.4 Кухре
правила [3] в глобально адаптивной схеме подразделения.
Алгоритм таков: пока не будет достигнута требуемая точность, разделите пополам
область с наибольшей ошибкой по оси с наибольшей четвертой
разница.
интерфейс, он такой же, как исходная подпрограмма DCUHRE
[11] . 4 Сравнение
безнадежный. Производительность сильно зависит от подынтегральной функции, и есть
всегда случаи, и не только академические, в которых одна рутинная
превосходит другие, или наоборот, в котором одна процедура просто дает
неправильные результаты. Это, конечно, главная причина, по которой существует четыре
независимые и легко взаимозаменяемые алгоритмы в библиотеке Cuba.
оценка, приведенная алгоритмами Монте-Карло, просто констатирует
интервал одной сигмы, или, другими словами, вероятность того, что центральный
значение, лежащее в данном интервале, составляет (всего) около 68%.
четыре подпрограммы Кубы на реальной интеграции фазового пространства. Результаты
полученные четырьмя методами (здесь не показаны), действительно согласуются в пределах
запрашиваемая точность. 5 Интерфейс Mathematica
библиотека в прямом смысле. Это скорее четыре исполняемых файла, которые
общаться с Mathematica через MathLink API:
\Text(50,85)[b]Математика
\CBox(0,0)(100,80)КрасныйПастельКрасный
\Text(50,65)[]\BlackVegas[f, …]\Text(50,30)[]\Blackintegrand f\Text(50,15)[]\Black(скомпилированная функция)\SetOffset(135,0)
\Текст(50,85)[b]С
\CBox(0,0)(100,80)СинийПастельСиний
\Text(50,65)[]\Blackvoid Vegas(…)\Text(50,25)[]\Образцы Blackrequest\SetOffset(105,0)
\SetWidth2.5
\ДлиннаяСтрелка(-36,65)(36,65)
\Text(0,71)[b]MathLink\Text(0,34)[b]{→x1,→x2,…}\LongArrow(36,30)(-36,30)
\ДлиннаяСтрелка(-36,20)(36,20)
\Text(0,16)[t]{f1,f2,…}
Install[«Vegas»], соответствующую процедуру можно использовать практически
как родной NIntegrate для Mathematica. Подынтегральная функция
полностью оценивается в Mathematica, что означает, что можно делать вещи
как Cuhre[Zeta[x y], {x,2,3}, {y,4,5}]
6 Дополнительные инструменты
6.1 Селектор
вызов: подпрограмма Куба(метод, ndim, ncomp,
подынтегральная функция, интеграл, ошибка, проба)
Вегас, Суаве, Дивонн, Кюре. Все параметры индивидуальны
подпрограммы «скрыты» (определены внутри подпрограммы), т. е. это не
готовый продукт, но должен быть адаптирован пользователем. 6.2 Средство просмотра разделов
алгоритм. Иногда это полезно для визуализации того, где подынтегральная функция
лежат характерные области. Это действительно полезно только в малом
умеренные габариты однако.
вывод передается через утилиту partview, например myprogram | часть 1 2
На рис. 2 показан снимок экрана.
Рисунок 2: Скриншот программы просмотра разделов Кубы.
7 Резюме
Cuba – библиотека для многомерного численного интегрирования, написанная
в C. Он содержит четыре независимых алгоритма: Vegas, Suave, Divonne,
и Cuhre, которые имеют аналогичные заклинания и могут быть легко заменены на
сравнение. Все подпрограммы могут интегрировать векторные интегралы и иметь
Интерфейс Fortran, C/C++ и Mathematica. Дополнительные инструменты
включены, такие как универсальный вызов и средство просмотра разделов. Куба
доступен на http://www.feynarts.de/cuba под лицензией
LGPL и прост в сборке (autoconf).
Каталожные номера
- [1]
- [2]
Т. Хан, комп. физ. коммун. 168 (2005) 78 [hep-ph/0404043].
- [3]
А. Генц, А. Малик, SIAM J. Numer. Анальный. 20 (1983) 580.
- [4]
В.Х. Пресс, С.А. Теукольский, В.Т. Феттерлинг, Б.П. Фланнери, Числовой
рецепты на Фортране, 2-е издание, Cambridge University Press,
1992. - [5]
П. Братли, Б.Л. Фокс, ACM Trans. Мат. Программное обеспечение 14 (1988) 88.
Алгоритм ТОМС 659. - [6]
М. Мацумото, Т. Нисимура, ACM Trans. Моделирование Комп. Моделирование 8 (1998) 3 и
http://www.math.sci.hiroshima-u.ac.jp/~m-mat/MT/emt.html. - [7]
Х. Л. Кенг, В. Юань, Приложения теории чисел к численному анализу,
Спрингер-Верлаг, 1981. - [8]
Г.П. Лепаж, Ж. Комп. физ. 27 (1978) 192.
Г.П. Лепаж, Отчет CLNS-80/447, Корнельский университет, Итака, Нью-Йорк, 1980. - [9]
В.Х. Пресс, Г.Р. Фаррар, Комп. в физ. 4 (1990) 190.
- [10]
Дж. Х. Фридман, М.Х. Райт, ACM Trans. Мат. Программное обеспечение 7 (1981) 76.
Дж.Х. Фридман, М. Х. Райт, отчет SLAC CGTM-193-РЕВ, ЦГТМ-193, 1981 г. - [11]
Дж. Бернтсен, Т. Эспелид, А. Генц, ACM Trans. Мат. Программное обеспечение 17 (1991) 437 и 452.
Алгоритм ТОМС 698.
Стиль изделия Белый / XS | РевЗилла Артикул № 1564991 | МФР. Товар № 4518913360 | ДоступностьНет в наличии: введите свой адрес электронной почты ниже, чтобы получить уведомление, когда этот товар появится в наличии. |
Стиль продукта Белый / SM | РевЗилла Артикул № 1564992 | МФР. Товар № 1325472 | ДоступностьНет в наличии: введите свой адрес электронной почты ниже, чтобы получить уведомление, когда этот товар появится в наличии. |
Стиль продукта Белый / MD | РевЗилла Артикул № 1564993 | МФР. Товар № 1325509 | ДоступностьНет в наличии: введите свой адрес электронной почты ниже, чтобы получить уведомление, когда этот товар появится в наличии. |
Стиль продукта Белый / LG | РевЗилла Артикул № 1564994 | МФР. Товар № 1325549 | ДоступностьНет в наличии: введите свой адрес электронной почты ниже, чтобы получить уведомление, когда этот товар появится в наличии. |
Стиль изделия Белый / XL | РевЗилла Артикул № 1564995 | МФР. Товар № 1325442 | ДоступностьНет в наличии: введите свой адрес электронной почты ниже, чтобы получить уведомление, когда этот товар появится в наличии. |
Стиль изделия Белый / 2XL | РевЗилла Артикул № 1564996 | МФР. Товар № 1325514 | ДоступностьНет в наличии: введите свой адрес электронной почты ниже, чтобы получить уведомление, когда этот товар появится в наличии. |
Стиль продукта Желтый / XS | РевЗилла Артикул № 1564997 | МФР. Товар № 4518923360 | ДоступностьНет в наличии: введите свой адрес электронной почты ниже, чтобы получить уведомление, когда этот товар появится в наличии. |
Стиль продукта Желтый / SM | РевЗилла Артикул № 1564998 | МФР. Товар № 1325415 | ДоступностьОсталось только 1: Доставка в течение 24 часов |
Стиль продукта Желтый / MD | РевЗилла Артикул № 1564999 | МФР. Товар № 1325499 | ДоступностьНет в наличии: введите свой адрес электронной почты ниже, чтобы получить уведомление, когда этот товар появится в наличии. |
Стиль продукта Желтый / LG | РевЗилла Артикул № 1565000 | МФР. Товар № 1325502 | ДоступностьНет в наличии: введите свой адрес электронной почты ниже, чтобы получить уведомление, когда этот товар появится в наличии. |
Стиль продукта Желтый / XL | РевЗилла Артикул № 1565001 | МФР. Товар №1325440 | ДоступностьНет в наличии: введите свой адрес электронной почты ниже, чтобы получить уведомление, когда этот товар появится в наличии. |
Стиль продукта Желтый / 2XL | РевЗилла Артикул № 1565002 | МФР. |