Не так давно мы выпустили профессиональный Квадрокоптер DEXP H-1
Но продаем его как “Радиоуправляемая игрушка квадрокоптер DEXP H-1”! Почему так спросите вы? Ну это все рекламный маркетинг виноват!!!!Коптер настолько новый и настолько профессиональный, что даже Яндекс в поиске по запросу “DEXP H-1” выдает товар на сайте DNS “Крепление DEXP h2”. Опять маркетологи DEXP не договорились и выпускают совершенно разный товар с одинаковым названием.
Пока не понятно для чего был запущен проморолик на youtube, с комментарием “у некоторых заинтересованных пользователей возникли вопросы”. У каких пользователей появились вопросы не понятно, на сайте DNS, основного продавца продукции DEXP, нет ни единого комментария по квадрокоптеру.
Не разбирайте данный АКБ, дабы не разочароваться в китайской емкости аккумуляторов
Очень жаль, что нет характеристик коптера DEXP H-1, так как на самом деле игрушка выглядит очень приятно и создает иллюзию действительно качественного продукта.
Пока об устройстве известно лишь, что максимальный радиус 500 метров (маловато для профессионала), имеет модуль GPS, автоматическая балансировка, полностью разборная модель, аккумулятор 3000 mAh, время работы от одного заряда 40 минут (явно завышено).начинка закреплена на металлическом кронштейне
В действительности коптер не дотягивает до профессионального дальностью полета, например DJI INSPIRE 1 улетает от пульта до 2 километров, нет встроенной камеры. Нет информации о таких необходимых вещах как: скорость и высота полета, подробности гироскопа, крепление камеры (питание, крепление), наличие карты памяти или сим карты, для чего ему нужен GPS (если сигнал не пишется, не передается), да что уж там даже вес квадрокоптера не указан, ни то что какой груз может поднять. На официальном сайте DEXP в обще нет никакой информации о характеристиках, сплошные картинки, ждем подробностей от покупателей.
ПДУ копия пультов на ОЕМ китайских коптерах: Wltoys V303, walkira qr x350
Как с таким подходом к продажам можно вообще говорить о продажах? Рекламный отдел DEXP видимо считает, что начав снимать ролик про квадрокоптер и выкладывая их на youtube, у них попрут продажи.
… ролик который даже с сайта DEXP нельзя посмотреть….А вы купите себе профессиональный квадрокоптер DEXP H-1? Игрушка или профессионал?
Ремонт Dexp Ursus Z210
Диагностика бесплатно (даже при отказе от ремонта)
Замена аудиокодека
Замена датчика приближения
Замена контроллера питания
Замена переключателя вибро
Замена разъема зарядки
Замена разъема карты памяти
Замена разъема наушников
Ремонт кнопок
Ремонт после воды
Быстро разряжается
Замена аккумулятора
Замена динамика
Замена дисплея
Замена задней крышки
Замена камеры
Замена кнопки включения
Замена контроллера зарядки
Замена корпусаЗамена крышки
Замена микрофона
Замена стекла
Замена шлейфа
Замена экрана
Кнопка Домой (Home)
Меня не слышат
Не видит карту памяти
Не видит Sim карту
Не видит наушники
Не гаснет экран
Не заряжается
Не работает wifi
Не работает виброзвонок
Не работает сенсор
Не регулируется громкость
Не слышно собеседника
Нет изображения
Нет подсветки
ПерезагружаетсяПолосит изображение
Разбито стекло
Ремонт антенны
Ремонт экрана
Сломан корпус
Сломана кнопка
Сломан разъем зарядки
Указанные цены на услуги действительны только при оформлении предварительной записи с сайта
Гарантийные сроки указаны на запчасти категории «ОРИГИНАЛ».
- Диагностика 0 ₽ (даже при отказе от ремонта)
- Замена экрана 840 ₽
- Замена стекла 840 ₽
- Замена разъема зарядки 720 ₽
- Замена аккумулятора 720 ₽
- Замена задней крышки 720 ₽
- Ремонт после воды 1190 ₽
- Ремонт кнопок 720 ₽
- Замена динамика 720 ₽
- Замена микрофона 720 ₽
- Замена камеры 840 ₽
- Замена кнопки включения 720 ₽
- Замена контроллера зарядки 1190 ₽
- Замена контроллера питания 1190 ₽
- Замена разъема наушников 840 ₽
- Замена корпуса 1190 ₽
- Замена шлейфа840 ₽
- Замена платы 1190 ₽
- Прошивка 990 ₽
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
| |||||
uz (ранее Torg) Ташкент
C dseed: исходное семя (dseed обновляется автоматически)
C nr: размерность массива двойной точности r
C r: реальный * 8 r (nr)
C GGUBS генерирует nr независимых однородных случайных величин
C GGNML генерирует nr независимых нормальных переменных ПОДПРОГРАММА GGUBS (DSEED, NR, R)
C ++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++ C
C GGUBS ГЕНЕРИРУЕТ NR ОДНОЙ ТОЧНОСТИ СЛУЧАЙНЫХ ВАРИАНТОВ ЕДИНИЦЫ C
C ON (0,1) ПО ЛИНЕЙНОЙ КОНГРУЕНЦИАЛЬНОЙ СХЕМЕ.
ЭТА ПРОЦЕДУРА C
C ЗАВИСИТ ОТ РАЗМЕРА СЛОВА МАШИНЫ. C
C ++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++ C
C ПОСЛЕДНЕЕ ИЗМЕНЕНИЕ C
C $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$, $$$$$$$$$$$$$$ C
C НА ВХОДЕ C
C DSEED ДВОЙНАЯ ТОЧНОСТЬ C
C СЕМЕНА ДЛЯ ГЕНЕРАТОРА C
C NR INTEGER C
C КОЛИЧЕСТВО ВАРИАНТОВ ОБРАЗОВАНИЯ C
C ПРИ ВОЗВРАЩЕНИИ C
C R REAL (NR) C
C ЕДИНАЯ ТОЧНАЯ МАССА, СОДЕРЖАЩАЯ ВАРИАТЫ C
C ЗВОНКИ GGUBS C
C DMOD C
C $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$, $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $$$$$$$$$$$$$$ C
INTEGER NR
НАСТОЯЩИЙ * 8 R (NR)
ДВОЙНАЯ ТОЧНОСТЬ DSEED
C
C МЕСТНОЕ
C
ЦЕЛОЕ I
ДВОЙНАЯ ТОЧНОСТЬ D2P31M, D2P31, DMULTX
C
C КОНСТАНТЫ МАШИНЫ
С D2P31M = (2 ** 31) – 1
C D2P31 = (2 ** 31) (ИЛИ НАСТРОЕННОЕ ЗНАЧЕНИЕ)
C
ДАННЫЕ D2P31M / 2147483647.D0 /
ДАННЫЕ D2P31 / 2147483711.
D0 /
ДАННЫЕ DMULTX / 16807.0D + 00 /
C
DO 5 I = 1, NR
DSEED = DMOD (DMULTX * DSEED, D2P31M)
R (I) = DSEED / D2P31
5 ПРОДОЛЖИТЬ
C
C КОНЕЦ GGUBS
C
ВОЗВРАЩЕНИЕ
КОНЕЦ
C ================================================= ============== C
C ================================================= ============== C
ПОДРОБНЫЙ GGNML (DSEED, NR, R)
C ++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++ C
C GGNML ГЕНЕРИРУЕТ NR ОДНОЙ ТОЧНОСТИ N (0,1) СЛУЧАЙНЫХ ВАРИАНТОВ C
C ++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++ C
C ПОСЛЕДНЕЕ ИЗМЕНЕНИЕ C
C $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$, $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $$$$$$$$$$$$$$ C
C НА ВХОДЕ C
C DSEED ДВОЙНАЯ ТОЧНОСТЬ C
C СЕМЕНА ДЛЯ ГЕНЕРАТОРА C
C NR INTEGER C
C КОЛИЧЕСТВО ВАРИАНТОВ ОБРАЗОВАНИЯ C
C ПРИ ВОЗВРАЩЕНИИ C
C R REAL (NR) C
C ЕДИНАЯ ТОЧНАЯ МАССА, СОДЕРЖАЩАЯ ВАРИАТЫ C
C $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$, $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $$$$$$$$$$$$$$ C
INTEGER NR
НАСТОЯЩИЙ * 8 R (NR)
ДВОЙНАЯ ТОЧНОСТЬ DSEED
C МЕСТНОЕ
INTEGER IER
C
C ПОЛУЧИТЬ NR СЛУЧАЙНЫЕ ЧИСЛА
C УНИФОРМА (0,1)
C
ВЫЗОВ GGUBS (DSEED, NR, R)
C
C ТРАНСФОРМАЦИЯ КАЖДОГО УНИФОРМА DEVIATE
C
DO 5 I = 1, NR
ПОЗВОНИТЬ MDNRIS (R (I), R (I), IER)
5 ПРОДОЛЖИТЬ
C
C КОНЕЦ GGNML
C
ВОЗВРАЩЕНИЕ
КОНЕЦ ВСПОМОГАТЕЛЬНЫЙ MDNRIS (P, Y, IER)
C ++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++ C
С С
C ++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++ C
C ПОСЛЕДНЕЕ ИЗМЕНЕНИЕ C
C $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$, $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $$$$$$$$$$$$$$ C
С С
C $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$, $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $$$$$$$$$$$$$$ C
C ТЕХНИЧЕСКИЕ ХАРАКТЕРИСТИКИ ДЛЯ АРГУМЕНТОВ
НАСТОЯЩИЙ * 8 P, Y
INTEGER IER
C ТЕХНИЧЕСКИЕ ХАРАКТЕРИСТИКИ ДЛЯ ЛОКАЛЬНЫХ ПЕРЕМЕННЫХ
НАСТОЯЩИЙ * 8 EPS, G0, G1, G2, G3, H0, h2, h3, A, W, WI, SN, SD
REAL * 8 SIGMA, SQRT2, X, XINF
ДАННЫЕ XINF / 1.
7014E + 38 /
ДАННЫЕ SQRT2 / 1.414214 /
ДАННЫЕ EPS / 1.1921E-07 /
ДАННЫЕ G0 / .1851159E-3 /, G1 / -. 2028152E-2 /
ДАННЫЕ G2 / -. 1498384 /, G3 / .1078639E-1 /
ДАННЫЕ H0 / .9952975E-1 /, h2 / .5211733 /
ДАННЫЕ h3 / -.6888301E-1 /
C ПЕРВОЕ ИСПОЛНИТЕЛЬНОЕ ЗАЯВЛЕНИЕ
IER = 0
ЕСЛИ (P .GT. 0.0. AND. P .LT. 1.0) ПЕРЕХОДИТЕ К ЭТАПУ 5
IER = 129 если (p .lt. 0,0D + 00), то
сигма = -1.0D + 00
еще
сигма = 1.0D + 00
endif C СИГМА = ЗНАК (1.0; P)
C написать (6,666) p, сигма
Формат C666 (‘Знак # 1:’, 2f15.8)
Пауза Y = SIGMA * XINF
ПЕРЕЙДИТЕ К 20
5 ЕСЛИ (P.LE.EPS) ПЕРЕЙДИТЕ К 10
Х = 1,0 – (P + P)
ПОЗВОНИТЬ MERFI (X, Y, IER)
Y = -SQRT2 * Y
ПЕРЕЙДИТЕ К 20
C P СЛИШКОМ МАЛЕНЬКИЙ, ВЫЧИТАЙТЕ Y НАПРЯМУЮ
10 А = П + П
W = КОРЕНЬ (-dLOG (A + (A-A * A)))
C ИСПОЛЬЗУЙТЕ РАЦИОНАЛЬНУЮ ФУНКЦИЮ В 1./ Вт
WI = 1 / Вт
SN = ((G3 * WI + G2) * WI + G1) * WI
SD = ((WI + h3) * WI + h2) * WI + H0
Y = W + W * (G0 + SN / SD)
Y = -Y * SQRT2
C КОНЕЦ MDNRIS
20 ВОЗВРАТ
КОНЕЦ
C ================================================= ============== C
ВСПОМОГАТЕЛЬНЫЙ MERFI (P, Y, IER)
C ++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++ C
С С
C ++++++++++++++++++++++++++++++++++++++++++++++++ ++++++++++++++ C
C ПОСЛЕДНЕЕ ИЗМЕНЕНИЕ C
C $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$, $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $$$$$$$$$$$$$$ C
С С
C $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$, $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $ $$$$$$$$$$$$$$ C
C ТЕХНИЧЕСКИЕ ХАРАКТЕРИСТИКИ ДЛЯ АРГУМЕНТОВ
НАСТОЯЩИЙ * 8 P, Y
INTEGER IER
C ТЕХНИЧЕСКИЕ ХАРАКТЕРИСТИКИ ДЛЯ ЛОКАЛЬНЫХ ПЕРЕМЕННЫХ
REAL * 8 A, B, X, Z, W, WI, SN, SD, F, Z2, RINFM, A1, A2, A3, B0, B1,
* B2, B3, C0, C1, C2, C3, D0, D1, D2, E0, E1, E2, E3, F0, F1,
* F2, G0, G1, G2, G3, H0, h2, h3, СИГМА
ДАННЫЕ A1 / -.
5751703 /, A2 / -1.896513 /, A3 / -. 5496261E-1 /
ДАННЫЕ B0 / -. 1137730 /, B1 / -3.293474 /, B2 / -2.374996 /
ДАННЫЕ B3 / -1.187515 /
ДАННЫЕ C0 / -. 1146666 /, C1 / -. 1314774 /, C2 / -. 2368201 /
ДАННЫЕ C3 / .5073975E-1 /
ДАННЫЕ D0 / -44.27977 /, D1 / 21.98546 /, D2 / -7.586103 /
ДАННЫЕ E0 / -.5668422E-1 /, E1 / .3937021 /, E2 / -. 3166501 /
ДАННЫЕ E3 / .6208963E-1 /
ДАННЫЕ F0 / -6.266786 /, F1 / 4.666263 /, F2 / -2.962883 /
ДАННЫЕ G0 / .1851159E-3 /, G1 / -. 2028152E-2 /
ДАННЫЕ G2 / -. 1498384 /, G3 / .1078639E-1 /
ДАННЫЕ H0 / .9952975E-1 /, h2 / .5211733 /
ДАННЫЕ h3 / -.6888301E-1 /
ДАННЫЕ RINFM / 1.7014E + 38 /
C ПЕРВОЕ ИСПОЛНИТЕЛЬНОЕ ЗАЯВЛЕНИЕ
IER = 0
X = P если (x. lt. 0,0D + 00), то
сигма = -1.0D + 00
еще
сигма = 1.0D + 00
endif C СИГМА = ЗНАК (1.0; X)
C написать (6,666) x, сигма
Формат C666 (‘Знак # 2:’, 2f15.8)
Пауза C ТЕСТ НА НЕДОСТАТОЧНЫЙ АРГУМЕНТ
ЕСЛИ (.
NOT. (X.GT.-1. .AND. X.LT.1.)) ПЕРЕЙДИТЕ К ЭТАПУ 30
Z = ABS (X)
ЕСЛИ (Z.LE. .85) ПЕРЕЙДИТЕ К 20
А = 1.-Z
B = Z
C СНИЖЕННЫЙ АРГУМЕНТ В (.85,1.),
C ПОЛУЧИТЬ ПРЕОБРАЗОВАННУЮ ПЕРЕМЕННУЮ
5 Вт = КОРЕНЬ (-dLOG (A + A * B))
ЕСЛИ (W.LT.2.5) ПЕРЕЙДИТЕ К 15
ЕСЛИ (W.LT.4.) ПЕРЕЙДИТЕ К ЭТАПУ 10
C W БОЛЬШЕ, ЧЕМ 4., ПРИМ. F BY A
РАЦИОНАЛЬНАЯ ФУНКЦИЯ В 1./W
WI = 1 / Вт
SN = ((G3 * WI + G2) * WI + G1) * WI
SD = ((WI + h3) * WI + h2) * WI + H0
F = W + W * (G0 + SN / SD)
ПЕРЕЙДИТЕ К 25
C W МЕЖДУ 2.5 И 4., ПРИМ. F
C РАЦИОНАЛЬНОЙ ФУНКЦИЕЙ В W
10 SN = ((E3 * W + E2) * W + E1) * W
SD = ((W + F2) * W + F1) * W + F0
F = W + W * (E0 + SN / SD)
ПЕРЕЙДИТЕ К 25
C W МЕЖДУ 1.13222 И 2.5, ПРИМ.
C F ПО РАЦИОНАЛЬНОЙ ФУНКЦИИ В W
15 СН = ((C3 * W + C2) * W + C1) * W
SD = ((W + D2) * W + D1) * W + D0
F = W + W * (C0 + SN / SD)
ПЕРЕЙДИТЕ К 25
C Z МЕЖДУ 0.И .85, ПРИМ. F
C РАЦИОНАЛЬНОЙ ФУНКЦИЕЙ В Z
20 Z2 = Z * Z
F = Z + Z * (B0 + A1 * Z2 / (B1 + Z2 + A2 / (B2 + Z2 + A3 / (B3 + Z2))))
C СОЗДАЙТЕ РЕШЕНИЕ ПО MULT.
F BY
C ПОДХОДЯЩИЙ ЗНАК
25 Y = СИГМА * F
IER = 0
ПЕРЕЙДИТЕ К 40
C ОШИБКА ВЫХОД.УСТАНОВИТЬ РЕШЕНИЕ НА ПЛЮС
C (ИЛИ МИНУС) БЕСКОНЕЧНОСТЬ
30 IER = 129
Y = СИГМА * RINFM
C КОНЕЦ MERFI
40 ВОЗВРАТ
КОНЕЦ
C ================================================= ============== C ВСПОМОГАТЕЛЬНЫЙ DDNOR (Y, GAUSS)
НЕПРАВИЛЬНОЕ РЕАЛЬНОЕ * 8 (A-H, O-Z)
НАСТОЯЩИЙ * 8 P (6), Q (5), A (9), B (8), C (5), D (4)
ДАННЫЕ P (1) / – 6.58749161529837803157D-04 /,
1 П (2) / – 1.60837851487422766278D-02 /,
2 П (3) / – 1.25781726111229246204D-01 /,
3 П (4) / – 3.60344899949804439429D-01 /,
4 П (5) / – 3.05326634961232344035D-01 /,
5 П (6) / – 1.63153871373020978498D-02 /
ДАННЫЕ Q (1) /2.33520497626869185443D-03/,
1 Q (2) /6.05183413124413191178D-02/,
2 Q (3) / 5.27905102951428412248D-01/,
3 Q (4) /1.87295284992346047209D00/,
4 Q (5) / 2.56852019228982242072D00 /
ДАННЫЕ A (1) /1.23033935479799725272D03/,
1 А (2) /2.
05107837782607146532D03/,
2 А (3) /1.71204761263407058314D03/,
3 А (4) / 8.81952221241769090411D02/,
4 А (5) / 2.98635138197400131132D02/,
5 А (6) /6.61191906371416294775D01/,
6 А (7) /8.88314979438837594118D00/,
7 А (8) /5.64188496988670089180D-01/,
8 А (9) /2.15311535474403846343D-08/
ДАННЫЕ B (1) /1.23033935480374942043D03/,
1 В (2) / 3.43936767414372163696D03 /,
2 В (3) /4.3626190
24715820D03/, 3 В (4) /3.29079923573345962678D03/, 4 В (5) /1.6213895745666
74D03/, 5 В (6) /5.37181101862009857509D02/, 6 В (7) / 1.17693950891312499305D02/, 7 В (8) /1.57449261107098347253D01/ ДАННЫЕ C (1) /3.209377589138469472562D03/, 1 С (2) /3.774852376853020208137D02/, 2 С (3) /1.138641541510501556495D02/, 3 С (4) /3.161123743870565596947D00/, 4 С (5) / 1.857777061846031526730D-01 / ДАННЫЕ D (1) /2.844236833439170622273D03/, 1 Д (2) /1.282616526077372275645D03/, 2 D (3) /2.
440246379344441733056D02/,
3 D (4) /2.360129095234412093499D01/
ДАННЫЕ ORPI / .5641895835477562869483D0 /,
1 ROOT2 / .70710678118654752440083D0 /
C В ЭТОЙ ПОДПРОГРАММЕ ИСПОЛЬЗУЕТСЯ МЕТОД CODY ДЛЯ ОЦЕНКИ КУМУЛЯТИВНОЙ
C НОРМАЛЬНОЕ РАСПРЕДЕЛЕНИЕ. ВЕРОЯТНО ТОЧНО ДО 19 ИЛИ 20
ЗНАЧИТЕЛЬНЫЕ ЦИФРЫ. ЭТО БЫЛ НАПИСАН ДЖЕЙМСОМ МАКИННОМ.
C 1977 г., НА ОСНОВАНИИ СТАТЬИ КОДИРОВКИ, УКАЗАННОЙ В ДОКУМЕНТАЦИИ
C ДЛЯ IMSL SUBROUTINE MDNOR.ISW = 1
ЕСЛИ (Y.LT.-16.D0) Y = -16.D0
ЕСЛИ (Y.GT.16.D0) Y = 16.D0
X = -Y * ROOT2
ЕСЛИ (X.GT.0.D0) ПЕРЕЙДИТЕ К 1
ЕСЛИ (X.LT.0.D0) ПЕРЕХОДИТЕ К ЭТАПУ 2
GAUSS = .5D0
ВОЗВРАЩЕНИЕ
2 ПРОДОЛЖИТЬ
Х = – Х
ISW = -1
1 ПРОДОЛЖИТЬ
ЕСЛИ (X.LT..477D0) ПЕРЕЙДИТЕ К 10
ЕСЛИ (X.LE.4.D0) ПЕРЕЙДИТЕ К 20
C ОЦЕНИТЬ ERFC ДЛЯ X.GT.4.0
Х2 = Х * Х
XM2 = 1.D0 / X2
XM4 = XM2 * XM2
XM6 = XM4 * XM2
XM8 = XM4 * XM4
XM10 = XM6 * XM4
ВЕРХ = P (1) + P (2) * XM2 + P (3) * XM4 + P (4) * XM6 + P (5) * XM8 + P (6) * XM10
БОТ = Q (1) + Q (2) * XM2 + Q (3) * XM4 + Q (4) * XM6 + Q (5) * XM8 + XM10
ДЕРЬМО = ORPI + TOP / (BOT * X2)
ERFC = DEXP (-X2) * ДЕРЬМО / X
C
ЕСЛИ (ISW.
EQ.-1) ERFC = 2.D0 – ERFC
GAUSS = ERFC * .5D0
ВОЗВРАЩЕНИЕ
20 ПРОДОЛЖИТЬ
C ОЦЕНКА ERFC ДЛЯ .477.LT.X.LE.4.0
Х2 = Х * Х
Х3 = Х2 * Х
Х4 = Х2 * Х2
Х5 = Х3 * Х2
Х6 = Х3 * Х3
Х7 = Х3 * Х4
Х8 = Х4 * Х4
ВЕРХ = A (1) + A (2) * X + A (3) * X2 + A (4) * X3 + A (5) * X4 + A (6) * X5 +
1 А (7) * Х6 + А (8) * Х7 + А (9) * Х8
БОТ = B (1) + B (2) * X + B (3) * X2 + B (4) * X3 + B (5) * X4 + B (6) * X5 +
1 В (7) * Х6 + В (8) * Х7 + Х8
ERFC = DEXP (-X2) * TOP / BOT
C
ЕСЛИ (ISW.EQ.-1) ERFC = 2.D0 – ERFC
GAUSS = ERFC * .5D0
ВОЗВРАЩЕНИЕ
10 ПРОДОЛЖИТЬ
C ОЦЕНИТЬ ERF ДЛЯ X.LT..477
Х2 = Х * Х
Х4 = Х2 * Х2
Х6 = Х4 * Х2
Х8 = Х4 * Х4
ВЕРХ = C (1) + C (2) * X2 + C (3) * X4 + C (4) * X6 + C (5) * X8
БОТ = D (1) + D (2) * X2 + D (3) * X4 + D (4) * X6 + X8
ERF = X * TOP / BOT
C
ERF = ERF * ISW
ERFC = 1.D0 – ERF
GAUSS = ERFC * .
5D0
ВОЗВРАЩЕНИЕ
КОНЕЦСЕК.gov | Превышен порог скорости запросов
Чтобы обеспечить равный доступ для всех пользователей, SEC оставляет за собой право ограничивать запросы, исходящие от необъявленных автоматизированных инструментов. Ваш запрос был идентифицирован как часть сети автоматизированных инструментов за пределами допустимой политики и будет обрабатываться до тех пор, пока не будут приняты меры по объявлению вашего трафика.
Пожалуйста, объявите свой трафик, обновив свой пользовательский агент, чтобы включить в него информацию о компании.
Для лучших практик по эффективной загрузке информации из SEC.gov, включая последние документы EDGAR, посетите sec.gov/developer. Вы также можете подписаться на рассылку обновлений по электронной почте о программе открытых данных SEC, включая передовые методы, которые делают загрузку данных более эффективной, и улучшения SEC.gov, которые могут повлиять на процессы загрузки по сценарию. Для получения дополнительной информации обращайтесь по адресу opendata@sec.
gov.
Для получения дополнительной информации см. Политику конфиденциальности и безопасности веб-сайта SEC. Благодарим вас за интерес к Комиссии по ценным бумагам и биржам США.
Идентификатор ссылки: 0.5dfd733e.1632650270.b68de029
Дополнительная информация
Политика безопасности в Интернете
Используя этот сайт, вы соглашаетесь на мониторинг и аудит безопасности. В целях безопасности и обеспечения того, чтобы общедоступная услуга оставалась доступной для пользователей, эта правительственная компьютерная система использует программы для мониторинга сетевого трафика для выявления несанкционированных попыток загрузки или изменения информации или иного причинения ущерба, включая попытки отказать пользователям в обслуживании.
Несанкционированные попытки загрузить информацию и / или изменить информацию в любой части этого сайта строго запрещены и подлежат судебному преследованию в соответствии с Законом о компьютерном мошенничестве и злоупотреблениях 1986 года и Законом о защите национальной информационной инфраструктуры 1996 года (см.
Раздел 18 U.S.C. §§ 1001 и 1030).
Чтобы обеспечить хорошую работу нашего веб-сайта для всех пользователей, SEC отслеживает частоту запросов на контент SEC.gov, чтобы гарантировать, что автоматический поиск не влияет на возможность доступа других лиц к контенту SEC.gov. Мы оставляем за собой право блокировать IP-адреса, которые отправляют чрезмерные запросы. Текущие правила ограничивают пользователей до 10 запросов в секунду, независимо от количества машин, используемых для отправки запросов.
Если пользователь или приложение отправляет более 10 запросов в секунду, дальнейшие запросы с IP-адреса (-ов) могут быть ограничены на короткий период.Как только количество запросов упадет ниже порогового значения на 10 минут, пользователь может возобновить доступ к контенту на SEC.gov. Эта практика SEC предназначена для ограничения чрезмерного автоматического поиска на SEC.gov и не предназначена и не ожидается, чтобы повлиять на людей, просматривающих веб-сайт SEC.gov.
Обратите внимание, что эта политика может измениться, поскольку SEC управляет SEC.
gov, чтобы гарантировать, что веб-сайт работает эффективно и остается доступным для всех пользователей.
Примечание: Мы не предлагаем техническую поддержку для разработки или отладки процессов загрузки по сценарию.
использовать portlib
неявное действительное * 8 (a-h, o-z)
размер h (20), r (20,20)
open (5, file = “out”)
write (0, *) ‘Дать m (количество переменных)’
99 прочтений *, м
если (m.eq.0) остановить
ч = 1.
г = 0,2
х = времяf ()
сделать 1 я = 1, м
г (я, я) = 1
1 продолжить
вызов dmv (m, 20, h, r, prob, 0.000001d0, ier, err, 15d0)
t = время f () – x
написать (0,2) m, проблема, t
написать (5,2) m, проблема, t
2 формата (i3, f10.7, ф11.2)
перейти к 99
конец ПОДПРОГРАММА DMV (M, K, H, R, PROB, EPS, IER, ERR, CUT)
НЕЯВНАЯ ДВОЙНАЯ ТОЧНОСТЬ (A-H, O-Z)
C
C
C DMV РАСЧЕТАЕТ МНОГООБРАЗНЫЙ НОРМАЛЬНЫЙ ИНТЕГРАЛ.
C
C ИНТЕГРАЛЬНЫЕ НИЖНИЕ ПРЕДЕЛЫ ЯВЛЯЮТСЯ -БЕСКОНЕЧНЫМИ ДЛЯ ВСЕХ ПЕРЕМЕННЫХ. C ВЕРХНИЕ ПРЕДЕЛЫ ПРЕДОСТАВЛЯЮТСЯ ВЕКТОРОМ H (ВХОД).
C МАТРИЦА КОРРЕЛЯЦИИ R (ВХОД).
C
ПАРАМЕТРЫ C:
C
C M – ЧИСЛО ПЕРЕМЕННЫХ. M10 (ВЫХОД).C
C ВЫРЕЗАТЬ – ЕСЛИ X0, И
C ЗНАЧЕНИЕ 2 ЕСЛИ БЕСКОНЕЧНОСТЬ ИСПОЛЬЗУЕТСЯ КАК ВЕРХНИЙ ПРЕДЕЛ И
ПЕРЕМЕННАЯ C УДАЛЯЕТСЯ ИЗ ИНТЕГРАЦИИ.
C LIST1 – СПИСОК ПЕРЕМЕННЫХ, УЧАСТВУЮЩИХ В
C ИНТЕГРАЦИЯ (НЕ ИМЕЮЩАЯ БЕСКОНЕЧНОСТЬ КАК ВЕРХНИЙ ПРЕДЕЛ).
C
ДО 10 I = 1, М
СПИСОК (I) = 0
ЕСЛИ (H (I) .GT.0.) СПИСОК (I) = 1
10 ПРОДОЛЖИТЬ
20 IPOS = 0
L = 0
ДО 30 I = 1, М
ЕСЛИ (СПИСОК (I) .EQ.2) ПЕРЕЙДИТЕ К 30
L = L + 1
СПИСОК1 (L) = I
30 ПРОДОЛЖИТЬ
ЕСЛИ (Л.EQ.0) ПЕРЕЙДИТЕ К 50
ДО 40 I = 1, L
L1 = СПИСОК1 (I)
h2 (I) = H (L1)
ЕСЛИ (СПИСОК (L1) .EQ.1) h2 (I) = – h2 (I)
ЕСЛИ (СПИСОК (L1) .EQ.1) IPOS = 1-IPOS
DO 40 Дж = 1, L
L2 = СПИСОК1 (J)
R1 (I, J) = R (L1, L2)
ЕСЛИ (СПИСОК (L1) .EQ.1) R1 (I, J) = – R1 (I, J)
ЕСЛИ (СПИСОК (L2) .
EQ.1) R1 (I, J) = – R1 (I, J)
40 ПРОДОЛЖИТЬ
50 S = РЕШИТЬ (L, K, h2, R1, ВЫРЕЗАТЬ, IER)
ЕСЛИ (IER.GT.0) ВОЗВРАТ
ЕСЛИ (IPOS.EQ.1) S = -S
DMV1 = DMV1 + S
ДО 60 I = 1, М
ЕСЛИ (СПИСОК (I).EQ.0) ПЕРЕЙДИТЕ К 60
СПИСОК (I) = СПИСОК (I) +1
ЕСЛИ (СПИСОК (I) .EQ.2) ПЕРЕЙДИТЕ К 20
СПИСОК (I) = 1
60 ПРОДОЛЖИТЬ
ВОЗВРАЩЕНИЕ
КОНЕЦ
РЕШЕНИЕ ФУНКЦИИ С ДВОЙНОЙ ТОЧНОСТЬЮ (M, K, H, R, CUT, IER)
C
C SOLVE РАСЧЕТАЕТ МНОГОМАРНУЮ ВЕРОЯТНОСТЬ ДАННОГО ВЕКТОРА H
C И ДАННЫЙ К.
C
C IER ЯВЛЯЕТСЯ ИНДИКАТОРОМ ОШИБКИ.
C IER = 0 ЯВЛЯЕТСЯ ОБЫЧНЫМ ПРЕКРАЩЕНИЕМ.
C IER = 98 УКАЗЫВАЕТ, ЧТО M ЗА ПРЕДЕЛАМИ ДИАПАЗОНА.
C IER = 99 означает, что R является единственным или его определяющим фактором является отрицательный.C ВСЕ ДРУГИЕ ПАРАМЕТРЫ ОПИСАНЫ ВЫШЕ В DMV1 ИЛИ DMV.
C SOLVE МОЖНО ВЫЗВАТЬ ПРЯМО. ИНТЕГРАЦИЯ
C ПАРАМЕТРЫ, ИЗВЛЕЧЕННЫЕ ИЗ [7], ИМЕЮТ ТОЧНОСТЬ 15 ЗНАЧЕНИЙ.
C
НЕЯВНАЯ ДВОЙНАЯ ТОЧНОСТЬ (A-H, O-Z)
РАЗМЕР H (20), h2 (20), h3 (20), СПИСОК (20), R (20,20), R1 (20,20)
КОЭФ (10,10), X (10,10), Y (20), D (2)
C
C КОЭФФИЦИЕНТЫ ВЕКТОРОВ И X СОДЕРЖАТ КВАДРАТУРНЫЕ КОЭФФИЦИЕНТЫ ГАУССА.
C И ТОЧКИ СОБСТВЕННО И ПОЛУЧАЮТСЯ ИЗ [6].C
ДАННЫЕ COEF / 10 * 0.D0, .6405291796843790D0, .245697745768379D0,8 * 0.D0
*,. 446029770466658D0, .396468266998335D0,4.37288879877644D-2,7 * 0.D0,
* .325302999756919D0, .421107101852062D0, .133442500357520D0,
* .637432348625728D-2,6 * 0.D0, .248406152028443D0, .392331066652399D0,
* .211418193076057D0, .332466603513439D-1, .824853344515628D-3,5 * 0.D0,
* .196849675488598D0, .349154201525395D0, .257259520584421D0,
* .0760131375840057D0, .685191862513596D-2,9.84716452019267D-5,
* 4 * 0.D0, .160609965149261D0, .306319808158099D0, .275527141784905D0,
* .120630193130784D0, .0218922863438067D0, .00123644672831056D0,
* 1.10841575911059D-5,3 * 0.D0, .13410918845336D0, .26833075447264D0,
* .275953397988422D0, .15744828261879D0, .0448141099174625D0,
* .00536793575602526D0,2.02063649132407D-4,1.19259692659532D-6,
* 2 * 0.D0, .114088970242118D0, .235940791223685D0, .266425473630253D0,
* .
183251679101663D0, .0713440493066916D0,.0139814184155604D0,
* .00116385272078519D0, .305670214897831D-4,1.23790511337496D-7,
* 0.D0,0.0985520975191087D0, .208678066608185D0, .252051688403761D0,
* .198684340038387D0, .097198422760062D0, .0270244164355446D0,
* .00380464962249537D0,2.28886243044656D-4,4.34534479844469D-6,
* 1.24773714817825D-8 /
ДАННЫЕ X / 10 * 0.D0, .300193931060839D0,1.25242104533372D0,8 * 0.D0,
* .190554149798192D0, .848251867544577D0,1.79977657841573D0,7 * 0.D0,
*.133776446996068D0, .6243246
190D0,1.34253782564499D0,
* 2.26266447701036D0,6 * 0.D0, .100242151968216D0, .482813966046201D0,
* 1.06094982152572D0,1.77972941852026D0,2.66976035608766D0,5 * 0.D0
*,. 0786006594130979D0, .386739410270631D0, .866429471682044D0,
* 1.46569804966352D0,2.172707796939D0,3.03682016932287D0,4 * 0.D0,
* .0637164846067008D0, .318192018888619D0, .724198989258373D0,
* 1.23803559921509D0,1.83852822027095D0,2.
53148815132768D0,
* 3.37345643012458D0,3 * 0.D0, .0529786439318514D0, .267398372167767D0,
* .616302884182402D0,1.06424631211623D0,1.58885586227006D0
*, 2.18392115309586D0,2.86313388370808D0,3.6860071627244D0,2 * 0.D0,
* .0449390308011934D0, .228605305560535D0, .532195844331646D0,
* .927280745338081D0,1.39292385519588D0,1.91884309919743D0,
* 2.50624783400574D0,3.17269213348124D0,3.97889886978978D0,0.D0,
* .0387385243257289D0, .198233304013083D0, .465201111814767D0
* ,.816861885592273D0,1.23454132402818D0,1.70679814968913D0,
* 2.22994008892494D0,2.80910374689875D0,3.46387241949586D0,
* 4.25536180636608D0 /
PI = 4.0D0 * DATAN (1.D0)
IER = 0
ЕСЛИ (M.LT.0.OR.M.GT.20) IER = 98
ЕСЛИ (IER.GT.0) ВОЗВРАТ
ЕСЛИ (K.LT.2.OR.K.GT.10) K = 10
ЕСЛИ (M.GT.0) ПЕРЕЙДИТЕ К 10
РЕШИТЬ = 1.
IER = 0
ВОЗВРАЩЕНИЕ
10 ДО 20 I = 1, М
DO 20 Дж = 1, М
20 R1 (I, J) = R (I, J)
ПОЗВОНИТЬ DPOFA (R1,20, M, IER)
ЕСЛИ (IER.
GT.0) IER = 99
ЕСЛИ (IER.EQ.99) ВОЗВРАТ
ПОЗВОНИТЬ DPODI (R1,20, M, D, 11)
ЕСЛИ ((D (1) .LE.0.). OR. (D (2) .LE.-30)) IER = 99
ЕСЛИ (IER.EQ.99) ВОЗВРАТ
C
C ОБРАТИТЕ ВНИМАНИЕ, ЧТО ДРУГАЯ ПОДПРОГРАММА МОЖЕТ ДАВАТЬ ОПРЕДЕЛЕНИЕ НАПРЯМУЮ
C В ОБНАРУЖЕНИИ, ПОЭТОМУ СЛЕДУЮЩЕЕ ЗАЯВЛЕНИЕ МОЖЕТ БЫТЬ НЕОБХОДИМО.
C R1 СОДЕРЖИТ ОБРАТКУ R. LOOP 30 CATERS ВОЗМОЖНОСТИ
C ЧТО ИНВЕРСИЯ СОХРАНЯЕТСЯ ТОЛЬКО В НИЖНЕМ ТРЕУГОЛЬНИКЕ.
C
DET = D (1) * 10.D0 ** D (2)
ДО 30 I = 1, М
DO 30 Дж = I, M
30 R1 (J, I) = R1 (I, J)
ПРОД = 1.D0
ДО 40 I = 1, М
h2 (I) = DSQRT (2./R1 (I, I))
h3 (I) = H (I) / h2 (I)
ПРОД = ПРОД * ПИ * R1 (I, I)
40 СПИСОК (I) = 1
ПРОД = 1.D0 / DSQRT (DET * PROD)
ДО 50 I = 1, М
DO 50 Дж = 1, М
50 R1 (I, J) = R1 (I, J) * h2 (I) * h2 (Дж) * 0.5D0
РЕШИТЬ = 0.
60 СУММ = 0.
ДО 70 I = 1, М
L = СПИСОК (I)
СУММ = СУММ + X (L, K) * X (L, K)
70 Y (I) = X (L, K) -h3 (I)
ДО 80 I = 1, М
DO 80 Дж = 1, М
80 СУММ = СУММ-Y (I) * Y (J) * R1 (I, J)
ЕСЛИ (СУММ.
LT.-CUT) ПЕРЕЙДИТЕ К 100
СУММ = DEXP (СУММ)
ДО 90 I = 1, М
90 СУММ = СУММ * КОЭФ (СПИСОК (I), K)
РЕШИТЬ = РЕШИТЬ + СУММ
100 ДО 110 I = 1, М
СПИСОК (I) = СПИСОК (I) +1
ЕСЛИ (СПИСОК (I) .LE.K) ПЕРЕХОДИТЕ К ЭТАПУ 60
СПИСОК (I) = 1
110 ПРОДОЛЖИТЬ
РЕШИТЬ = РЕШИТЬ * ПРОД
ВОЗВРАЩЕНИЕ
КОНЕЦ
подпрограмма dpofa (a, lda, n, info)
целое lda, n, информация
двойная точность a (lda, 1)
c
c dpofa множитель симметричный положительно определенный двойной точности
c матрица.c
c dpofa обычно вызывается dpoco, но его можно назвать
c напрямую с экономией времени, если rcond не нужен.
c (время для dpoco) = (1 + 18 / n) * (время для dpofa).
c
c при входе
c
c двойная точность (lda, n)
c симметричная матрица, подлежащая факторизации. только
c используются диагональ и верхний треугольник.
c
c lda целое число
c ведущее измерение массива a.c
c n целое число
c порядок матрицы a.
c
c по возвращении
c
c a верхнетреугольная матрица r, так что a = trans (r) * r
c, где trans (r) – транспонирование.
c строгий нижний треугольник не изменился.
c, если info .ne. 0 факторизация не завершена.
c
c информация целое число
c = 0 для нормального возврата.
c = k сигнализирует о состоянии ошибки.ведущий минор
c порядка k не является положительно определенным.
c
c linpack. эта версия от 14.08.78.
c Клив Молер, университет Нью-Мексико, национальная лаборатория аргонна.
c
c подпрограммы и функции
c
c blas ddot
c fortran dsqrt
c
c внутренние переменные
c
ddot двойной точности, т
двойная точность s
целое число j, jm1, k
c начать блок с … выход на 40
c
c
сделать 30 j = 1, n
info = j
s = 0.0d0
jm1 = j – 1
если (jm1 .lt. 1) перейти к 20
сделать 10 k = 1, jm1
t = a (k, j) – ddot (k-1, a (1, k), 1, a (1, j), 1)
т = т / а (к, к)
а (к, j) = т
s = s + t * t
10 продолжить
20 продолжить
s = a (j, j) – s
c .
….. выход
если (s. le. 0,0d0) перейти к 40
а (j, j) = dsqrt (s)
30 продолжить
info = 0
40 продолжить
возвращение
конец
подпрограмма dpodi (a, lda, n, det, job)
целое lda, n, работа
двойная точность a (lda, 1)
двойная точность det (2)
c
c dpodi вычисляет определитель и инверсию некоторого
c симметричная положительно определенная матрица двойной точности (см. ниже)
c с использованием коэффициентов, вычисленных dpoco, dpofa или dqrdc.c
c при входе
c
c двойная точность (lda, n)
c вывод a из dpoco или dpofa
c или вывод x из dqrdc.
c
c lda целое число
c ведущее измерение массива a.
c
c n целое число
c порядок матрицы a.
c
c целое число
c = 11 как определитель, так и обратный.
c = 01 только инверсия.Только определитель c = 10.
c
c по возвращении
c
c a если dpoco или dpofa использовались для факторизации a, то
c dpodi производит верхнюю половину инверсии (a).
c, если dqrdc использовался для разложения x, тогда
c dpodi создает верхнюю половину инверсии (trans (x) * x)
c, где trans (x) – транспонирование.
c элементы a ниже диагонали не изменяются.
c, если цифра единиц задания равна нулю, а не изменяется.c
c det двойной точности (2)
c определитель a или trans (x) * x, если требуется.
c иначе не упоминается.
детерминант c = det (1) * 10.0 ** det (2)
c с 1.0 .le. det (1) .lt. 10.0
c или det (1) .eq. 0,0.
c
c состояние ошибки
c
c деление на ноль произойдет, если входной коэффициент содержит
c ноль по диагонали и требуется обратное.c этого не произойдет, если подпрограммы вызываются правильно
c и если dpoco или dpofa установили info .eq. 0.
c
c linpack. эта версия от 14.08.78.
c Клив Молер, университет Нью-Мексико, национальная лаборатория аргонна.
c
c подпрограммы и функции
c
c blas daxpy, dscal
c fortran мод
c
c внутренние переменные
c
двойная точность t
двойная точность s
целое число i, j, jm1, k, kp1
c
c вычислить определитель
c
если (работа / 10.экв. 0) перейти к 70
det (1) = 1.0d0
det (2) = 0,0d0
s = 10,0d0
сделать 50 я = 1, п
det (1) = a (i, i) ** 2 * det (1)
c .
.. выход
if (det (1) .eq. 0.0d0) перейти к 60
10 if (det (1) .ge. 1.0d0) перейти к 20
det (1) = s * det (1)
Дет (2) = Дет (2) – 1.0d0
перейти к 10
20 продолжить
30, если (det (1) .lt. S) перейти к 40
det (1) = det (1) / с
Дет (2) = Дет (2) + 1.0d0
перейти к 30
40 продолжить
50 продолжить
60 продолжить
70 продолжить
c
c вычислить обратное (r)
c
if (mod (job, 10) .eq. 0) перейти к 140
сделать 100 k = 1, n
а (к, к) = 1.0d0 / а (к, к)
т = -а (к, к)
вызовите dscal (k-1, t, a (1, k), 1)
kp1 = k + 1
если (n .lt. kp1) перейти к 90
do 80 j = kp1, n
т = а (к, j)
а (к, j) = 0.0d0
вызовите daxpy (k, t, a (1, k), 1, a (1, j), 1)
80 продолжить
90 продолжить
100 продолжить
c
c форма inverse (r) * trans (inverse (r))
c
сделать 130 j = 1, n
jm1 = j – 1
если (jm1 .
lt. 1) перейти к 120
сделать 110 k = 1, jm1
т = а (к, j)
вызовите daxpy (k, t, a (1, j), 1, a (1, k), 1)
110 продолжить
120 продолжить
т = а (j, j)
вызовите dscal (j, t, a (1, j), 1)
130 продолжить
140 продолжить
возвращение
конец
подпрограмма daxpy (n, da, dx, incx, dy, incy)
c
c константа, умноженная на вектор плюс вектор.c использует развернутые циклы для приращений, равных единице.
c джек донгарра, линпак, 11.03.78.
c
двойной точности dx (1), dy (1), da
целое число i, incx, incy, ixiy, m, mp1, n
c
если (n.le.0) return
if (da .eq. 0.0d0) возврат
if (incx.eq.1.and.incy.eq.1) перейдите к 20
c
c код для неравных или равных приращений
c не равно 1
c
ix = 1
iy = 1
если (incx.lt.0) ix = (-n + 1) * incx + 1
если (incy.lt.0) iy = (-n + 1) * incy + 1
сделать 10 я = 1, п
dy (iy) = dy (iy) + da * dx (ix)
ix = ix + incx
iy = iy + incy
10 продолжить
возвращение
c
c код для обоих приращений, равных 1
c
c
c цикл очистки
c
20 м = mod (n, 4)
если (m.
экв. 0) перейти к 40
do 30 i = 1, м
dy (i) = dy (i) + da * dx (i)
30 продолжить
if (n .lt. 4) return
40 мп1 = м + 1
сделать 50 я = mp1, n, 4
dy (i) = dy (i) + da * dx (i)
dy (i + 1) = dy (i + 1) + da * dx (i + 1)
dy (i + 2) = dy (i + 2) + da * dx (i + 2)
dy (i + 3) = dy (i + 3) + da * dx (i + 3)
50 продолжить
возвращение
конец
подпрограмма dscal (n, da, dx, incx)
c
c масштабирует вектор константой.c использует развернутые циклы для приращения, равного единице.
c джек донгарра, линпак, 11.03.78.
c
двойной точности da, dx (1)
целое число i, incx, m, mp1, n, nincx
c
если (n.le.0) return
если (incx.eq.1) перейти к 20
c
c код для приращения, не равного 1
c
nincx = n * incx
сделать 10 я = 1, nincx, incx
dx (i) = da * dx (i)
10 продолжить
возвращение
c
c код для приращения, равного 1
c
c
c цикл очистки
c
20 м = mod (n, 5)
если М .
экв. 0) перейти к 40
do 30 i = 1, м
dx (i) = da * dx (i)
30 продолжить
if (n .lt. 5) return
40 мп1 = м + 1
сделать 50 я = mp1, n, 5
dx (i) = da * dx (i)
dx (я + 1) = да * dx (я + 1)
dx (я + 2) = da * dx (я + 2)
dx (i + 3) = da * dx (i + 3)
dx (я + 4) = da * dx (я + 4)
50 продолжить
возвращение
конец
функция двойной точности ddot (n, dx, incx, dy, incy)
c
c образует скалярное произведение двух векторов.c использует развернутые циклы для приращений, равных единице.
c джек донгарра, линпак, 11.03.78.
c
двойная точность dx (1), dy (1), dtemp
целое число i, incx, incy, ix, iy, m, mp1, n
c
ddot = 0,0d0
dtemp = 0,0d0
если (n.le.0) return
if (incx.eq.1.and.incy.eq.1) перейдите к 20
c
c код для неравных или равных приращений
c не равно 1
c
ix = 1
iy = 1
если (incx.lt.0) ix = (-n + 1) * incx + 1
если (incy.lt.0) iy = (-n + 1) * incy + 1
сделать 10 я = 1, п
dtemp = dtemp + dx (ix) * dy (iy)
ix = ix + incx
iy = iy + incy
10 продолжить
ddot = dtemp
возвращение
c
c код для обоих приращений, равных 1
c
c
c цикл очистки
c
20 м = mod (n, 5)
если (m.
экв. 0) перейти к 40
do 30 i = 1, м
dtemp = dtemp + dx (я) * dy (я)
30 продолжить
если (п. 5) перейти к 60
40 мп1 = м + 1
сделать 50 я = mp1, n, 5
dtemp = dtemp + dx (i) * dy (i) + dx (i + 1) * dy (i + 1) +
* dx (i + 2) * dy (i + 2) + dx (i + 3) * dy (i + 3) + dx (i + 4) * dy (i + 4)
50 продолжить
60 ddot = dtemp
возвращение
конец
время подпрограммы (t, i)
реальная * 8 т, тт, ст
вызов gettim (i1, i2, i3, i4)
tt = i1 * 3600 + i2 * 60 + i3 + i4 * 0.01d0
если (ieq.0) st = tt
если (т.е.0) возврат
tt = tt-st
если (tt.lt.0.) tt = tt + 86400
т = тт
возвращение
конец
Корреляция свойств материалов: сравнение между FRAPCON-3.4, FRAPTRAN 1.4 и MATPRO – страница 124 из 143
С помощью программного обеспечения для оптического распознавания символов из изображения на этой странице был автоматически извлечен следующий текст:
cthcon вычисляет теплопроводность оболочки как функцию
c температуры, времени, потока и холодной обработки
c
c ccon = выходная теплопроводность циркалоя-4 (w / mk)
c cdkdt = выходная производная теплопроводности с
c относительно температуры
c
c ctemp = температура точки сетки на входе оболочки (k)
c time = время на входе при температуре и потоке (с)
c flux = входной поток быстрых нейтронов ((нейтронов / м ** 2) / с )
c coldw = входная холодная обработка (безразмерное соотношение площадей)
c
c перейти к соответствующей модели для типа облицовки
перейти к (400, 400, 400, 400, 400, 600, 400, 600, 400, 400, 400 ) icm
c
400 продолжить
c Общее свойство циркалоя из MATPRO
c
c Уравнение, используемое в этом подкоде, основано на данных из
c w.
k. Андерсон, К.Дж. Бек, А. kephart и j.s. theilacker
c astm-stp-314, 1962, стр 62 - 93
c j.n. чиригос, гр. kass, w.s. Кирк и Дж. Дж. salvaggio
c изготовление тепловыделяющих элементов, академическая пресса, 1961, стр 19 - 55
c a.d. feith, gemp-669 (октябрь 1966 г.)
c c.f. удачи и h.w. deem, bmi-1273 (1958) pp 7-9
c a.e. Пауэрс, капл-2146 (1961)
c д.б. scott, wcap-3269-41 (1965) pp 5-9
c
c Эта версия cthcon не использует время, флюс или холодную обработку
c в качестве параметров при расчете теплопроводности циркалоя
c
c Одно стандартное отклонение этой функции = 1.01 w / m-k
c
c Эта версия cthcon была разработана r.l. Миллер, декабрь 1975 г.
c Последнее изменение: e.r. carlson в июне 1978 г.
c Изменено PNNL, январь 1997 г., чтобы очистить кодирование и удалить
c анализ лицензирования и кодирование неопределенности чувствительности
c
ccon = 7,511d0 + ctemp * (2,088d-2 + ctemp * (-1,450d -5 +
& ctemp * 7. 668d-09))
cdkdt = 2.088d-2 + ctemp * (- 2.9d-5 + ctemp * 2.3d-8)
return
c
600 продолжить
c E-110 свойство как закодировано RRC-KI и предоставлено PNNL
c
c Этот кодовый блок вычисляет теплопроводность (w / (m * k)) и производную
c теплопроводности по температуре
c (w / (m * k * л)) из сплава h2.Данные взяты из:
c Волков Б.Ю. et. al. «Библиотека свойств материалов для плакирования из сплава h2
c», препринт IAE-4941/11, 1989.
c
if (ctemp.le.2133.0d0) then
ccon = dexp (0.000461843d0 * ctemp) * 15.0636d0
cdkdt = 0,006957018215d0 * dexp (0,000461843d0 * ctemp)
else
ccon = 36.d0
cdkdt = 0.d0
endif
return
c
end
A.22
k. Андерсон, К.Дж. Бек, А. kephart и j.s. theilacker 
668d-09)) DEXP h2 – Coub – Самая большая платформа видеомема
DEXP h2 – Coub – Самая большая платформа видеомема- Дом
- Горячий
- Случайный
Подробнее.
..Показать меньше
- Мне нравится
- Закладки
- Сообщества
Животные и домашние животные
Ведение блога
Стенд-ап и анекдоты
Мэшап
Аниме
Фильмы и сериалы
Игры
Мультфильмы
Искусство и дизайн
Живые изображения
Музыка
Новости и политика
Спорт
Наука и технологии
Еда и кухня
Знаменитости
Природа и путешествия
Мода и красота
танец
Авто и техника
Мемы
NSFW
Рекомендуемые
Коуб дня
-
Темная тема
..
