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