xref: /petsc/src/dm/interface/ftn-custom/zdmf.c (revision 7d5fd1e4d9337468ad3f05b65b7facdcd2dfd2a4)
1 #include <petsc/private/fortranimpl.h>
2 #include <petscdm.h>
3 #include <petscviewer.h>
4 
5 #if defined(PETSC_HAVE_FORTRAN_CAPS)
6 #define dmcreateinterpolation_       DMCREATEINTERPOLATION
7 #define dmview_                      DMVIEW
8 #define dmsetoptionsprefix_          DMSETOPTIONSPREFIX
9 #define dmsettype_                   DMSETTYPE
10 #define dmgettype_                   DMGETTYPE
11 #define dmsetmattype_                DMSETMATTYPE
12 #define dmsetvectype_                DMSETVECTYPE
13 #define dmgetmattype_                DMGETMATTYPE
14 #define dmgetvectype_                DMGETVECTYPE
15 #define dmlabelview_                 DMLABELVIEW
16 #define dmcreatelabel_               DMCREATELABEL
17 #define dmhaslabel_                  DMHASLABEL
18 #define dmgetlabelvalue_             DMGETLABELVALUE
19 #define dmsetlabelvalue_             DMSETLABELVALUE
20 #define dmgetlabelsize_              DMGETLABELSIZE
21 #define dmgetlabelidis_              DMGETLABELIDIS
22 #define dmgetlabelname_              DMGETLABELNAME
23 #define dmgetlabel_                  DMGETLABEL
24 #define dmgetstratumsize_            DMGETSTRATUMSIZE
25 #define dmgetstratumis_              DMGETSTRATUMIS
26 #define dmsetstratumis_              DMSETSTRATUMIS
27 #define dmremovelabel_               DMREMOVELABEL
28 #define dmviewfromoptions_           DMVIEWFROMOPTIONS
29 #define dmcreatesuperdm_             DMCREATESUPERDM
30 #define dmdestroy_                   DMDESTROY
31 #elif !defined(PETSC_HAVE_FORTRAN_UNDERSCORE)
32 #define dmcreateinterpolation_       dmcreateinterpolation
33 #define dmview_                      dmview
34 #define dmsetoptionsprefix_          dmsetoptionsprefix
35 #define dmsettype_                   dmsettype
36 #define dmgettype_                   dmgettype
37 #define dmsetmattype_                dmsetmattype
38 #define dmsetvectype_                dmsetvectype
39 #define dmgetmattype_                dmgetmattype
40 #define dmgetvectype_                dmgetvectype
41 #define dmlabelview_                 dmlabelview
42 #define dmcreatelabel_               dmcreatelabel
43 #define dmhaslabel_                  dmhaslabel
44 #define dmgetlabelvalue_             dmgetlabelvalue
45 #define dmsetlabelvalue_             dmsetlabelvalue
46 #define dmgetlabelsize_              dmlabelsize
47 #define dmgetlabelidis_              dmlabelidis
48 #define dmgetlabelname_              dmgetlabelname
49 #define dmgetlabel_                  dmgetlabel
50 #define dmgetstratumsize_            dmgetstratumsize
51 #define dmgetstratumis_              dmgetstratumis
52 #define dmsetstratumis_              dmsetstratumis
53 #define dmremovelabel_               dmremovelabel
54 #define dmviewfromoptions_           dmviewfromoptions
55 #define dmcreatesuperdm_             dmreatesuperdm
56 #define dmdestroy_                   dmdestroy
57 #endif
58 
59 PETSC_EXTERN void dmgetmattype_(DM *mm,char* name,PetscErrorCode *ierr,PETSC_FORTRAN_CHARLEN_T len)
60 {
61   const char *tname;
62 
63   *ierr = DMGetMatType(*mm,&tname);if (*ierr) return;
64   if (name != PETSC_NULL_CHARACTER_Fortran) {
65     *ierr = PetscStrncpy(name,tname,len);if (*ierr) return;
66   }
67   FIXRETURNCHAR(PETSC_TRUE,name,len);
68 }
69 
70 PETSC_EXTERN void dmgetvectype_(DM *mm,char* name,PetscErrorCode *ierr,PETSC_FORTRAN_CHARLEN_T len)
71 {
72   const char *tname;
73 
74   *ierr = DMGetVecType(*mm,&tname);if (*ierr) return;
75   if (name != PETSC_NULL_CHARACTER_Fortran) {
76     *ierr = PetscStrncpy(name,tname,len);if (*ierr) return;
77   }
78   FIXRETURNCHAR(PETSC_TRUE,name,len);
79 }
80 
81 PETSC_EXTERN void dmview_(DM *da,PetscViewer *vin,PetscErrorCode *ierr)
82 {
83   PetscViewer v;
84   PetscPatchDefaultViewers_Fortran(vin,v);
85   *ierr = DMView(*da,v);
86 }
87 
88 PETSC_EXTERN void dmsetoptionsprefix_(DM *dm,char* prefix, PetscErrorCode *ierr,PETSC_FORTRAN_CHARLEN_T len)
89 {
90   char *t;
91 
92   FIXCHAR(prefix,len,t);
93   *ierr = DMSetOptionsPrefix(*dm,t);if (*ierr) return;
94   FREECHAR(prefix,t);
95 }
96 
97 PETSC_EXTERN void dmsettype_(DM *x,char* type_name,PetscErrorCode *ierr,PETSC_FORTRAN_CHARLEN_T len)
98 {
99   char *t;
100 
101   FIXCHAR(type_name,len,t);
102   *ierr = DMSetType(*x,t);if (*ierr) return;
103   FREECHAR(type_name,t);
104 }
105 
106 PETSC_EXTERN void dmgettype_(DM *mm,char* name,PetscErrorCode *ierr,PETSC_FORTRAN_CHARLEN_T len)
107 {
108   const char *tname;
109 
110   *ierr = DMGetType(*mm,&tname);if (*ierr) return;
111   if (name != PETSC_NULL_CHARACTER_Fortran) {
112     *ierr = PetscStrncpy(name,tname,len);if (*ierr) return;
113   }
114   FIXRETURNCHAR(PETSC_TRUE,name,len);
115 
116 }
117 
118 PETSC_EXTERN void dmsetmattype_(DM *dm,char* prefix, PetscErrorCode *ierr,PETSC_FORTRAN_CHARLEN_T len)
119 {
120   char *t;
121 
122   FIXCHAR(prefix,len,t);
123   *ierr = DMSetMatType(*dm,t);if (*ierr) return;
124   FREECHAR(prefix,t);
125 }
126 
127 PETSC_EXTERN void dmsetvectype_(DM *dm,char* prefix, PetscErrorCode *ierr,PETSC_FORTRAN_CHARLEN_T len)
128 {
129   char *t;
130 
131   FIXCHAR(prefix,len,t);
132   *ierr = DMSetVecType(*dm,t);if (*ierr) return;
133   FREECHAR(prefix,t);
134 }
135 
136 PETSC_EXTERN void dmcreatelabel_(DM *dm, char* name, int *ierr,PETSC_FORTRAN_CHARLEN_T lenN)
137 {
138   char *lname;
139 
140   FIXCHAR(name, lenN, lname);
141   *ierr = DMCreateLabel(*dm, lname);if (*ierr) return;
142   FREECHAR(name, lname);
143 }
144 
145 PETSC_EXTERN void dmhaslabel_(DM *dm, char* name, PetscBool *hasLabel, int *ierr,PETSC_FORTRAN_CHARLEN_T lenN)
146 {
147   char *lname;
148 
149   FIXCHAR(name, lenN, lname);
150   *ierr = DMHasLabel(*dm, lname, hasLabel);if (*ierr) return;
151   FREECHAR(name, lname);
152 }
153 
154 PETSC_EXTERN void dmgetlabelvalue_(DM *dm, char* name, PetscInt *point, PetscInt *value, int *ierr,PETSC_FORTRAN_CHARLEN_T lenN)
155 {
156   char *lname;
157 
158   FIXCHAR(name, lenN, lname);
159   *ierr = DMGetLabelValue(*dm, lname, *point, value);if (*ierr) return;
160   FREECHAR(name, lname);
161 }
162 
163 PETSC_EXTERN void dmsetlabelvalue_(DM *dm, char* name, PetscInt *point, PetscInt *value, int *ierr,PETSC_FORTRAN_CHARLEN_T lenN)
164 {
165   char *lname;
166 
167   FIXCHAR(name, lenN, lname);
168   *ierr = DMSetLabelValue(*dm, lname, *point, *value);if (*ierr) return;
169   FREECHAR(name, lname);
170 }
171 
172 PETSC_EXTERN void dmgetlabelsize_(DM *dm, char* name, PetscInt *size, int *ierr,PETSC_FORTRAN_CHARLEN_T lenN)
173 {
174   char *lname;
175 
176   FIXCHAR(name, lenN, lname);
177   *ierr = DMGetLabelSize(*dm, lname, size);if (*ierr) return;
178   FREECHAR(name, lname);
179 }
180 
181 PETSC_EXTERN void dmgetlabelidis_(DM *dm, char* name, IS *ids, int *ierr,PETSC_FORTRAN_CHARLEN_T lenN)
182 {
183   char *lname;
184 
185   FIXCHAR(name, lenN, lname);
186   *ierr = DMGetLabelIdIS(*dm, lname, ids);if (*ierr) return;
187   FREECHAR(name, lname);
188 }
189 
190 PETSC_EXTERN void dmgetlabelname_(DM *dm,PetscInt *n,char* name,PetscErrorCode *ierr,PETSC_FORTRAN_CHARLEN_T len)
191 {
192   const char *tmp;
193   *ierr = DMGetLabelName(*dm,*n,&tmp);
194   *ierr = PetscStrncpy(name,tmp,len);if (*ierr) return;
195   FIXRETURNCHAR(PETSC_TRUE,name,len);
196 }
197 
198 PETSC_EXTERN void dmgetlabel_(DM *dm, char* name, DMLabel *label, int *ierr,PETSC_FORTRAN_CHARLEN_T lenN)
199 {
200   char *lname;
201 
202   FIXCHAR(name, lenN, lname);
203   *ierr = DMGetLabel(*dm, lname, label);if (*ierr) return;
204   FREECHAR(name, lname);
205 }
206 
207 PETSC_EXTERN void dmgetstratumsize_(DM *dm, char* name, PetscInt *value, PetscInt *size, int *ierr,PETSC_FORTRAN_CHARLEN_T lenN)
208 {
209   char *lname;
210 
211   FIXCHAR(name, lenN, lname);
212   *ierr = DMGetStratumSize(*dm, lname, *value, size);if (*ierr) return;
213   FREECHAR(name, lname);
214 }
215 
216 PETSC_EXTERN void dmgetstratumis_(DM *dm, char* name, PetscInt *value, IS *is, int *ierr,PETSC_FORTRAN_CHARLEN_T lenN)
217 {
218   char *lname;
219 
220   FIXCHAR(name, lenN, lname);
221   *ierr = DMGetStratumIS(*dm, lname, *value, is);if (*ierr) return;
222   if (is && !*is) *is = (IS)0;
223   FREECHAR(name, lname);
224 }
225 
226 PETSC_EXTERN void dmsetstratumis_(DM *dm, char* name, PetscInt *value, IS *is, int *ierr,PETSC_FORTRAN_CHARLEN_T lenN)
227 {
228   char *lname;
229 
230   FIXCHAR(name, lenN, lname);
231   *ierr = DMSetStratumIS(*dm, lname, *value, *is);if (*ierr) return;
232   FREECHAR(name, lname);
233 }
234 
235 PETSC_EXTERN void dmremovelabel_(DM *dm, char* name, DMLabel *label, int *ierr,PETSC_FORTRAN_CHARLEN_T lenN)
236 {
237   char *lname;
238 
239   FIXCHAR(name, lenN, lname);
240   *ierr = DMRemoveLabel(*dm, lname, label);if (*ierr) return;
241   FREECHAR(name, lname);
242 }
243 
244 PETSC_EXTERN void dmviewfromoptions_(DM *dm,PetscObject obj,char* type,PetscErrorCode *ierr,PETSC_FORTRAN_CHARLEN_T len)
245 {
246   char *t;
247 
248   FIXCHAR(type,len,t);
249   CHKFORTRANNULLOBJECT(obj);
250   *ierr = DMViewFromOptions(*dm,obj,t);if (*ierr) return;
251   FREECHAR(type,t);
252 }
253 
254 PETSC_EXTERN void dmcreateinterpolation_(DM *dmc,DM *dmf,Mat *mat,Vec *vec, int *ierr)
255 {
256   CHKFORTRANNULLOBJECT(vec);
257   *ierr = DMCreateInterpolation(*dmc,*dmf,mat,vec);
258 }
259 
260 PETSC_EXTERN void dmcreatesuperdm_(DM dms[], PetscInt *len, IS ***is, DM *superdm, int *ierr)
261 {
262   *ierr = DMCreateSuperDM(dms, *len, *is, superdm);
263 }
264 
265 PETSC_EXTERN void dmdestroy_(DM *x,int *ierr)
266 {
267   PETSC_FORTRAN_OBJECT_F_DESTROYED_TO_C_NULL(x);
268   *ierr = DMDestroy(x); if (*ierr) return;
269   PETSC_FORTRAN_OBJECT_C_NULL_TO_F_DESTROYED(x);
270 }
271