|
Текст подпрограммы и версий ( Фортран ) ias2r.zip |
Тексты тестовых примеров ( Фортран ) tias2r.zip |
|
Текст подпрограммы и версий ( Си ) ias2r_c.zip |
Тексты тестовых примеров ( Си ) tias2r_c.zip |
|
Текст подпрограммы и версий ( Паскаль ) ias2r_p.zip |
Тексты тестовых примеров ( Паскаль ) tias2r_p.zip |
Вычисление значения бикубического аппроксимирующего сплайна или его смешанных производных до третьего порядка при заданных значениях аргументов на заданных (возможно неравномерных) сетках при заданных коэффициентах.
B силу свойств бикубических сплайнов решение этой задачи приводится к решению пяти одномерных задач, решение которых изложено в описании подпрограммы - функции IAS1R.
Hа первых четырех этапах фиксируется индекс по Y, а на пятом этапе используются четыре полученных значения.
REAL FUNCTION IAS2R (TX, NX, TY, NY, CXY, NCX, NCY, IDX,
IDY, XX, YY)
Параметры
| TX - | вещественный вектоp длины NX + 4, содержащий значения узлов сетки сплайна по X; |
| NX - | заданное число узлов исходной сетки по X, NX ≥ 4 (тип: целый); |
| TY - | вещественный вектоp длины NY + 4, содержащий значения узлов сетки сплайна по Y; |
| NY - | заданное число узлов исходной сетки по Y, NY ≥ 4 (тип: целый); |
| CXY - | вещественный двумерный массив, содержащий значения коэффициентов сплайна, размерностью NCX * NCY; |
| NCX - | заданное число стpок массива CXY, NCX ≥ NX + 2 (тип: целый); |
| NCY - | заданное число столбцов массива CXY, NCY ≥ NY + 2 (тип: целый); |
| IDX - | заданный порядок вычисляемой производной по X, 0 ≤ IDX ≤ 3 (тип: целый); |
| IDY - | заданный порядок вычисляемой производной по Y, 0 ≤ IDY ≤ 3 (тип: целый); |
| XX - | заданное значение аргумента по X, при котоpом производятся вычисления, TX (3) ≤ XX ≤ TX (NX + 2) (тип: вещественный); |
| YY - | заданное значение аргумента по Y, при котоpом производятся вычисления, TY (3) ≤ YY ≤ TY (NY + 2) (тип: вещественный). |
Версии: нет
Вызываемые подпрограммы: нет
Замечания по использованию
|
Подпрограмма - функция использует рабочие подпрограммы-функции IAS1R1 и IAS1R2. Для вычисления сеток и коэффициентов сплайна рекомендуется воспользоваться подпрограммой IAC2R. |
REAL IAS2R
DIMENSION X(4), Y(4), FXY(4, 4), TX(8), TY(8), CY(8), CXY(8, 8),
* FXA(4), FXB(4)
DATA X /2., 3., 4., 5./, Y /2., 3., 5., 6./, FXA /24., 36., 60., 72./
DATA FXB /60., 90., 150., 180./
DO 10 I = 1, 4
DO 10 J = 1, 4
FXY(I, J) = Y(J) * (X(I)**3)
10 CONTINUE
CALL IAC2R (X, 4, Y, 4, FXY, 4, 4, 1, 1, 1, FXA, 2, FXB,
* 0, X, 0, X, TX, TY, CXY, 8, 8, CY)
S = IAS2R (TX, 4, TY, 4, CXY, 8, 8, 2, 0, 2.4, 3.6)
Результат: S = 51.840