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