Имеется фортрановская библиотека LAPACK - пакет работы с линейной алгеброй. Конкретно в данный момент интересует её функция DGTSV (метод прогонки для трёхдиагональных матриц). Её сигнатура в фортране такая:
subroutine dgtsv (
integer N,
integer NRHS,
double precision, dimension( * ) DL,
double precision, dimension( * ) D,
double precision, dimension( * ) DU,
double precision, dimension( ldb, * ) B,
integer LDB,
integer INFO
)
В компонентном паскале (он только 32-bit, поэтому версия LAPACK - тоже 32-bit) я подключал LAPACK таким образом (звездочка имеет смысл метки для видимости извне, аналог public в C#):
MODULE LinalgLAPACK ["liblapack.dll"];
PROCEDURE DGTSV* ["dgtsv_"] (VAR n : LAPACK_Integer; VAR nrhs : LAPACK_Integer;
VAR dl : Double; VAR d : Double; VAR du : Double; VAR b : Double;
VAR ldb : LAPACK_Integer; VAR info : LAPACK_Integer);
Тут есть нюанс - liblapack.dll требует libblas.dll, а еще libgfortran-3.dll, libgcc_s_dw2-1.dll и libquadmath-0.dll. Все они - 32-битные. Зависимости исследовались Dependency Walker'ом.
Поскольку в фортране ВСЕ аргументы передаются ТОЛЬКО по ссылке, поэтому аргументы-указатели на массивы заменены на ссылки на начальные элементы массивов. В основном модуле вызов выглядел так:
MODULE LinalgTest;
IMPORT LAPACK := LinalgLAPACK;
...
PROCEDURE Test*;
VAR
size, nrhs, info : INTEGER;
c0, c1, c2, v : POINTER TO ARRAY OF REAL;
BEGIN
size := 3;
nrhs := 1;
NEW (c0, size);
NEW (c1, size);
NEW (c2, size);
NEW (v, size);
...
LAPACK.DGTSV (size, nrhs, c0[1], c1[0], c2[0], v[0], size, info);
END Test;
Главная диагональ имеет длину size, над- и под- диагонали - соответственно size-1. Во избежание смещения индексов используется маленькая хитрость: по ссылке (а фортрановская DGTSV интерпретирует сcылку как указатель) передается элемент с индексом 1 из массива длиной size, а не элемент с индексом 0 из массива длиной size - 1. И все прекрасно работает...
Теперь пытаюсь проделать то же в C#. Используя скомпилированные в x64 .dll библиотеки - это x64 версия liblapack.dll, которая явно требует x64 версии libblas.dll, а неявно еще и libgfortran_64-3.dll, libgcc_s_seh_64-1.dll, libquadmath-0.dll, решил разобраться с простым консольным x64 приложением на C#.
using System;
using System.Runtime.InteropServices;
namespace Linalg
{
class Program
{
[DllImport("liblapack.dll", CharSet = CharSet.Ansi, SetLastError = true)]
public static extern void dgtsv_(in int size1, in int nrhs, in double c0, in double c1, in double c2, in double v, in int size2, out int info);
static void Main(string[] args)
{
int size = 3;
int nrhs = 1;
double[] c0 = new double[] { 0.0, 1.0, 0.0 };
double[] c1 = new double[] { 1.0, -2.0, 1.0 };
double[] c2 = new double[] { 0.0, 1.0, 0.0 };
double[] v = new double[] { 1.0, 2.0, 3.0 };
int info = 0;
dgtsv_(in size, in nrhs, in c0[1], in c1[0], in c2[0], in v[0], in size, out info);
Console.WriteLine("info = {0}", info);
}
}
}
При запуске - System.BadImageFormatException: "Была сделана попытка загрузить программу, имеющую неверный формат. (0x8007000B)"
Люди, добрые - что не так? )))