У меня есть программа, состоящая в основном из Fortran 77 с оболочкой C++, которая читает и записывает в базу данных. Две части приложения совместно используют данные, используя функцию, которая, если вы используете глобальную структуру C/C++, названную как общий блок Fortran. Я почти уверен, что этот подход к интеграции C++/Fortran с псевдонимами 1) работает для многих наборов компиляторов, 2) не является стандартным. Я стараюсь поддерживать свой код так, чтобы он максимально использовал стандартные компоненты. Кроме того, эта интеграция оказалась хрупкой.
В утд.ч:
/*************************************************************
* This struct must follow the common block points.
* See the Fortran include file points.i
*************************************************************/
typedef struct ALIGN points
{
double point[3][MAX_PTS];
double edge[3][MAX_PTS];
double edgenorm[3][MAX_PTS];
double edgerho[MAX_PTS];
int nfacets[MAX_PTS];
double facet1[3][MAX_PTS];
double facet2[3][MAX_PTS];
double gaussk[MAX_PTS];
int sbstin[MAX_PTS];
int coatin[MAX_PTS];
int sbstout[MAX_PTS];
int coatout[MAX_PTS];
int ncorners[MAX_PTS];
double cnrpoint[3][MAX_CNRS][MAX_PTS];
int ncnredgs[MAX_CNRS][MAX_PTS];
double cnredge[3][MAX_CNREDS][MAX_CNRS][MAX_PTS];
int cnrsbst[MAX_CNREDS][MAX_CNRS][MAX_PTS];
int cnrcoat[MAX_CNREDS][MAX_CNRS][MAX_PTS];
int npoints;
} POINTS;
extern POINTS points_;
В утд.cpp:
POINTS points_;
В пунктах.i:
! maxpnt - maximum number of points in a path.
integer maxpnt
parameter ( maxpnt = 1000 )
integer npoints, nfacets(maxpnt)
integer ncorners(maxpnt), ncnredgs(maxpnt,maxcorners)
integer sbstin(maxpnt), coatin(maxpnt)
integer sbstout(maxpnt), coatout(maxpnt)
double precision point(maxpnt,3)
double precision edge(maxpnt,3), edgenorm(maxpnt,3)
double precision edgerho(maxpnt)
double precision facet1(maxpnt,3), facet2(maxpnt,3)
double precision gaussk(maxpnt), cnrpoint(maxpnt,maxcorners,3)
double precision cnredge(maxpnt,maxcorners,maxcnredges,3)
integer cnrsbst(maxpnt,maxcorners,maxcnredges)
integer cnrcoat(maxpnt,maxcorners,maxcnredges)
common /points/ point, edge, edgenorm, edgerho,
* nfacets, facet1, facet2, gaussk,
* sbstin, coatin, sbstout, coatout,
* ncorners, cnrpoint, ncnredgs,
* cnredge, cnrsbst, cnrcoat,
* npoints
Есть ли способ лучше? Я мог бы преобразовать общие блоки в модули. Тогда я уверен, что этот бизнес с псевдонимом общего блока с глобальной структурой невозможен. Как вы создаете модуль Fortran из C++? Как вы читаете данные из такого модуля?
Что вы посоветуете по поводу интеграции C++/Fortran?