Actual source code: petscmat.h

  1: !
  2: !  $Id: petscmat.h,v 1.73 2001/08/06 21:19:26 bsmith Exp $;
  3: !
  4: !  Include file for Fortran use of the Mat package in PETSc
  5: !
  6: #if !defined (__PETSCMAT_H)

  9: #define Mat PetscFortranAddr
 10: #define MatFDColoring PetscFortranAddr
 11: #define MatPartitioning PetscFortranAddr
 12: #define MatAIJIndices PetscFortranAddr
 13: #define MatType character*(80)
 14: #define MatOption integer 
 15: #define MatAssemblyType integer
 16: #define MatOrderingType character*(80)
 17: #define MatSORType integer
 18: #define MatInfoType integer
 19: #define MatReuse integer
 20: #define MatOperation integer
 21: #define MatColoringType character*(80)
 22: #define MatInfo double precision
 23: #define MatILUInfo double precision
 24: #define MatLUInfo double precision
 25: #define MatCholeskyInfo double precision
 26: #define MatDuplicateOption integer      
 27: #define MatStructure integer
 28: #define MatPartitioningType character*(80)
 29: #define MatNullSpace PetscFortranAddr

 31: #define MatStencil integer
 32: #define MatStencil_k 1
 33: #define MatStencil_j 2
 34: #define MatStencil_i 3
 35: #define MatStencil_c 4

 37: #define MATPARTITIONING_CURRENT 'current'
 38: #define MATPARTITIONING_PARMETIS 'parmetis'

 40: #define MATCOLORING_NATURAL 'natural'
 41: #define MATCOLORING_SL 'sl'
 42: #define MATCOLORING_LF 'lf'
 43: #define MATCOLORING_ID 'id'

 45: #define MATORDERING_NATURAL 'natural'
 46: #define MATORDERING_ND 'nd'
 47: #define MATORDERING_1WD '1wd'
 48: #define MATORDERING_RCM 'rcm'
 49: #define MATORDERING_QMD 'qmd'
 50: #define MATORDERING_ROWLENGTH 'rowlength'
 51: #define MATORDERING_DSC_ND 'dsc_nd'
 52: #define MATORDERING_DSC_MMD 'dsc_mmd'
 53: #define MATORDERING_DSC_MDF 'dsc_mdf'

 55: #endif

 57: !
 58: !  Matrix types
 59: !
 60: #define MATSAME     'same'
 61: #define MATSEQMAIJ  'seqmaij'
 62: #define MATMPIMAIJ  'mpimaij'
 63: #define MATIS       'is'
 64: #define MATMPIROWBS 'mpirowbs'
 65: #define MATSEQDENSE 'seqdense'
 66: #define MATSEQAIJ   'seqaij'
 67: #define MATMPIAIJ   'mpiaij'
 68: #define MATSHELL    'shell'
 69: #define MATSEQBDIAG 'seqbdiag'
 70: #define MATMPIBDIAG 'mpibdiag'
 71: #define MATMPIDENSE 'mpidense'
 72: #define MATSEQBAIJ  'seqbaij'
 73: #define MATMPIBAIJ  'mpibaij'
 74: #define MATMPIADJ   'mpiadj'
 75: #define MATSEQSBAIJ 'seqsbaij'
 76: #define MATMPISBAIJ 'mpisbaij'


 79: #if !defined (PETSC_AVOID_DECLARATIONS)

 81: !
 82: !  Flag for matrix assembly
 83: !
 84:       integer MAT_FLUSH_ASSEMBLY,MAT_FINAL_ASSEMBLY

 86:       parameter(MAT_FLUSH_ASSEMBLY=1,MAT_FINAL_ASSEMBLY=0)
 87: !
 88: !  Matrix options; must match those in include/petscmat.h
 89: !
 90:       integer MAT_ROW_ORIENTED,MAT_COLUMN_ORIENTED,MAT_ROWS_SORTED
 91:       integer MAT_COLUMNS_SORTED,MAT_NO_NEW_NONZERO_LOCATIONS
 92:       integer MAT_YES_NEW_NONZERO_LOCATIONS,MAT_SYMMETRIC
 93:       integer MAT_STRUCTURALLY_SYMMETRIC,MAT_NO_NEW_DIAGONALS
 94:       integer MAT_YES_NEW_DIAGONALS,MAT_INODE_LIMIT_1
 95:       integer MAT_INODE_LIMIT_2,MAT_INODE_LIMIT_3,MAT_INODE_LIMIT_4
 96:       integer MAT_INODE_LIMIT_5,MAT_IGNORE_OFF_PROC_ENTRIES
 97:       integer MAT_ROWS_UNSORTED,MAT_COLUMNS_UNSORTED
 98:       integer MAT_NEW_NONZERO_LOCATION_ERR
 99:       integer MAT_NEW_NONZERO_ALLOCATION_ERR,MAT_USE_HASH_TABLE
100:       integer MAT_KEEP_ZEROED_ROWS,MAT_IGNORE_ZERO_ENTRIES
101:       integer MAT_USE_INODES,MAT_DO_NOT_USE_INODES

103:       parameter (MAT_ROW_ORIENTED=1,MAT_COLUMN_ORIENTED=2)
104:       parameter (MAT_ROWS_SORTED=4,MAT_COLUMNS_SORTED=8)
105:       parameter (MAT_NO_NEW_NONZERO_LOCATIONS=16)
106:       parameter (MAT_YES_NEW_NONZERO_LOCATIONS=32)
107:       parameter (MAT_SYMMETRIC=64,MAT_STRUCTURALLY_SYMMETRIC=65)
108:       parameter (MAT_NO_NEW_DIAGONALS=66,MAT_YES_NEW_DIAGONALS=67)
109:       parameter (MAT_INODE_LIMIT_1=68,MAT_INODE_LIMIT_2=69)
110:       parameter (MAT_INODE_LIMIT_3=70,MAT_INODE_LIMIT_4=71)
111:       parameter (MAT_INODE_LIMIT_5=72,MAT_IGNORE_OFF_PROC_ENTRIES=73)
112:       parameter (MAT_ROWS_UNSORTED=74,MAT_COLUMNS_UNSORTED=75)
113:       parameter (MAT_NEW_NONZERO_LOCATION_ERR=76)
114:       parameter (MAT_NEW_NONZERO_ALLOCATION_ERR=77)
115:       parameter (MAT_USE_HASH_TABLE=78)
116:       parameter (MAT_KEEP_ZEROED_ROWS=79)
117:       parameter (MAT_IGNORE_ZERO_ENTRIES=80)
118:       parameter (MAT_USE_INODES=81,MAT_DO_NOT_USE_INODES=82)
119: !
120: !  MatDuplicateOption
121: !
122:       integer MAT_DO_NOT_COPY_VALUES,MAT_COPY_VALUES
123:       parameter (MAT_DO_NOT_COPY_VALUES=0,MAT_COPY_VALUES=1)
124: !
125: !  Flags for PCSetOperators()
126: !
127:       integer SAME_NONZERO_PATTERN,DIFFERENT_NONZERO_PATTERN
128:       integer SAME_PRECONDITIONER

130:       parameter (SAME_NONZERO_PATTERN = 0,DIFFERENT_NONZERO_PATTERN = 1)
131:       parameter (SAME_PRECONDITIONER = 2)

133: !
134: !  Note: MAT_INFO_SIZE must equal # elements in MatInfo structure
135: !  (See petsc/include/petscmat.h)
136: !
137:       integer   MAT_INFO_SIZE

139:       parameter (MAT_INFO_SIZE=14)

141:       integer MAT_INFO_ROWS_GLOBAL,MAT_INFO_COLUMNS_GLOBAL
142:       integer MAT_INFO_ROWS_LOCAL,MAT_INFO_COLUMNS_LOCAL
143:       integer MAT_INFO_BLOCK_SIZE,MAT_INFO_NZ_ALLOCATED
144:       integer MAT_INFO_NZ_USED,MAT_INFO_NZ_UNNEEDED
145:       integer MAT_INFO_MEMORY,MAT_INFO_ASSEMBLIES
146:       integer MAT_INFO_MALLOCS,MAT_INFO_FILL_RATIO_GIVEN
147:       integer MAT_INFO_FILL_RATIO_NEEDED,MAT_INFO_FACTOR_MALLOCS

149:       parameter (MAT_INFO_ROWS_GLOBAL=1,MAT_INFO_COLUMNS_GLOBAL=2)
150:       parameter (MAT_INFO_ROWS_LOCAL=3,MAT_INFO_COLUMNS_LOCAL=4)
151:       parameter (MAT_INFO_BLOCK_SIZE=5,MAT_INFO_NZ_ALLOCATED=6)
152:       parameter (MAT_INFO_NZ_USED=7,MAT_INFO_NZ_UNNEEDED=8)
153:       parameter (MAT_INFO_MEMORY=9,MAT_INFO_ASSEMBLIES=10)
154:       parameter (MAT_INFO_MALLOCS=11,MAT_INFO_FILL_RATIO_GIVEN=12)
155:       parameter (MAT_INFO_FILL_RATIO_NEEDED=13)
156:       parameter (MAT_INFO_FACTOR_MALLOCS=14)
157: !
158: !  MatReuse
159: !
160:       integer MAT_INITIAL_MATRIX,MAT_REUSE_MATRIX

162:       parameter (MAT_INITIAL_MATRIX=0,MAT_REUSE_MATRIX=1)

164: !
165: !  MatInfoType
166: !
167:       integer MAT_LOCAL,MAT_GLOBAL_MAX,MAT_GLOBAL_SUM

169:       parameter (MAT_LOCAL=1,MAT_GLOBAL_MAX=2,MAT_GLOBAL_SUM=3)

171: !
172: !  Note: MAT_ILUINFO_SIZE must equal # elements in MatILUInfo structure
173: !  (See petsc/include/petscmat.h)
174: !
175:       integer   MAT_ILUINFO_SIZE

177:       parameter (MAT_ILUINFO_SIZE=109)

179:       integer MAT_ILUINFO_LEVELS
180:       integer MAT_ILUINFO_FILL
181:       integer MAT_ILUINFO_DIAGONAL_FILL
182:       integer MAT_ILUINFO_DT
183:       integer MAT_ILUINFO_DTCOL
184:       integer MAT_ILUINFO_DTCOUNT
185:       integer MAT_ILUINFO_DAMPING
186:       integer MAT_ILUINFO_DAMP
187:       integer MAT_ILUINFO_PIVOT_ZERO
188:       integer MAT_ILUINFO_PIVOT_IN_BLOCKS

190:       parameter (MAT_ILUINFO_LEVELS = 1)
191:       parameter (MAT_ILUINFO_FILL = 2)
192:       parameter (MAT_ILUINFO_DIAGONAL_FILL = 3)
193:       parameter (MAT_ILUINFO_DT = 4)
194:       parameter (MAT_ILUINFO_DTCOL = 5)
195:       parameter (MAT_ILUINFO_DTCOUNT = 6)
196:       parameter (MAT_ILUINFO_DAMPING = 7)
197:       parameter (MAT_ILUINFO_DAMP = 8)
198:       parameter (MAT_ILUINFO_PIVOT_ZERO = 9)
199:       parameter (MAT_ILUINFO_PIVOT_IN_BLOCKS = 10)

201: !
202: !  Note: MAT_LUINFO_SIZE must equal # elements in MatLUInfo structure
203: !  (See petsc/include/petscmat.h)
204: !
205:       integer   MAT_LUINFO_SIZE

207:       parameter (MAT_LUINFO_SIZE = 6)

209:       integer MAT_LUINFO_FILL
210:       integer MAT_LUINFO_DTCOL
211:       integer MAT_LUINFO_DAMPING
212:       integer MAT_LUINFO_DAMP
213:       integer MAT_LUINFO_PIVOT_ZERO
214:       integer MAT_LUINFO_PIVOT_IN_BLOCKS

216:       parameter (MAT_LUINFO_FILL  = 1)
217:       parameter (MAT_LUINFO_DTCOL = 2)
218:       parameter (MAT_LUINFO_DAMPING = 3)
219:       parameter (MAT_LUINFO_DAMP = 4)
220:       parameter (MAT_LUINFO_PIVOT_ZERO = 5)
221:       parameter (MAT_LUINFO_PIVOT_IN_BLOCKS = 6)


224: !
225: !  Note: MAT_CHOLESKYINFO_SIZE must equal # elements in MatCholeskyInfo structure
226: !  (See petsc/include/petscmat.h)
227: !
228:       integer   MAT_CHOLESKYINFO_SIZE

230:       parameter (MAT_CHOLESKYINFO_SIZE = 4)

232:       integer MAT_CHOLESKYINFO_FILL
233:       integer MAT_CHOLESKYINFO_DAMPING
234:       integer MAT_CHOLESKYINFO_DAMP
235:       integer MAT_CHOLESKYINFO_PIVOTINBLOCKS

237:       parameter (MAT_CHOLESKYINFO_FILL  = 1)
238:       parameter (MAT_CHOLESKYINFO_DAMPING = 2)
239:       parameter (MAT_CHOLESKYINFO_DAMP = 3)
240:       parameter (MAT_CHOLESKYINFO_PIVOTINBLOCKS = 4)

242: !
243: !  Options for SOR and SSOR
244: !  MatSorType may be bitwise ORd together, so do not change the numbers
245: !
246:       integer SOR_FORWARD_SWEEP,SOR_BACKWARD_SWEEP,SOR_SYMMETRIC_SWEEP
247:       integer SOR_LOCAL_FORWARD_SWEEP,SOR_LOCAL_BACKWARD_SWEEP
248:       integer SOR_LOCAL_SYMMETRIC_SWEEP,SOR_ZERO_INITIAL_GUESS
249:       integer SOR_EISENSTAT,SOR_APPLY_UPPER,SOR_APPLY_LOWER

251:       parameter (SOR_FORWARD_SWEEP=1,SOR_BACKWARD_SWEEP=2)
252:       parameter (SOR_SYMMETRIC_SWEEP=3,SOR_LOCAL_FORWARD_SWEEP=4)
253:       parameter (SOR_LOCAL_BACKWARD_SWEEP=8)
254:       parameter (SOR_LOCAL_SYMMETRIC_SWEEP=12)
255:       parameter (SOR_ZERO_INITIAL_GUESS=16,SOR_EISENSTAT=32)
256:       parameter (SOR_APPLY_UPPER=64,SOR_APPLY_LOWER=128)
257: !
258: !  MatOperation
259: !
260:       integer MATOP_SET_VALUES
261:       integer MATOP_GET_ROW
262:       integer MATOP_RESTORE_ROW
263:       integer MATOP_MULT
264:       integer MATOP_MULT_ADD
265:       integer MATOP_MULT_TRANSPOSE
266:       integer MATOP_MULT_TRANSPOSE_ADD
267:       integer MATOP_SOLVE
268:       integer MATOP_SOLVE_ADD
269:       integer MATOP_SOLVE_TRANSPOSE
270:       integer MATOP_SOLVE_TRANSPOSE_ADD
271:       integer MATOP_LUFACTOR
272:       integer MATOP_CHOLESKYFACTOR
273:       integer MATOP_RELAX
274:       integer MATOP_TRANSPOSE
275:       integer MATOP_GETINFO
276:       integer MATOP_EQUAL
277:       integer MATOP_GET_DIAGONAL
278:       integer MATOP_DIAGONAL_SCALE
279:       integer MATOP_NORM
280:       integer MATOP_ASSEMBLY_BEGIN
281:       integer MATOP_ASSEMBLY_END
282:       integer MATOP_COMPRESS
283:       integer MATOP_SET_OPTION
284:       integer MATOP_ZERO_ENTRIES
285:       integer MATOP_ZERO_ROWS
286:       integer MATOP_LUFACTOR_SYMBOLIC
287:       integer MATOP_LUFACTOR_NUMERIC
288:       integer MATOP_CHOLESKY_FACTOR_SYMBOLIC
289:       integer MATOP_CHOLESKY_FACTOR_NUMERIC
290:       integer MATOP_GET_SIZE
291:       integer MATOP_GET_LOCAL_SIZE
292:       integer MATOP_GET_OWNERSHIP_RANGE
293:       integer MATOP_ILUFACTOR_SYMBOLIC
294:       integer MATOP_ICCFACTOR_SYMBOLIC
295:       integer MATOP_GET_ARRAY
296:       integer MATOP_RESTORE_ARRAY

298:       integer MATOP_CONVERT_SAME_TYPE
299:       integer MATOP_FORWARD_SOLVE
300:       integer MATOP_BACKWARD_SOLVE
301:       integer MATOP_ILUFACTOR
302:       integer MATOP_ICCFACTOR
303:       integer MATOP_AXPY
304:       integer MATOP_GET_SUBMATRICES
305:       integer MATOP_INCREASE_OVERLAP
306:       integer MATOP_GET_VALUES
307:       integer MATOP_COPY
308:       integer MATOP_PRINT_HELP
309:       integer MATOP_SCALE
310:       integer MATOP_SHIFT
311:       integer MATOP_DIAGONAL_SHIFT
312:       integer MATOP_ILUDT_FACTOR
313:       integer MATOP_GET_BLOCK_SIZE

315:       integer MATOP_GET_ROW_IJ
316:       integer MATOP_RESTORE_ROW_IJ
317:       integer MATOP_GET_COLUMN_IJ
318:       integer MATOP_RESTORE_COLUMN_IJ
319:       integer MATOP_FDCOLORING_CREATE
320:       integer MATOP_COLORING_PATCH
321:       integer MATOP_SET_UNFACTORED
322:       integer MATOP_PERMUTE
323:       integer MATOP_SET_VALUES_BLOCKED

325:       integer MATOP_DESTROY
326:       integer MATOP_VIEW

328:       parameter(MATOP_SET_VALUES=0)
329:       parameter(MATOP_GET_ROW=1)
330:       parameter(MATOP_RESTORE_ROW=2)
331:       parameter(MATOP_MULT=3)
332:       parameter(MATOP_MULT_ADD=4)
333:       parameter(MATOP_MULT_TRANSPOSE=5)
334:       parameter(MATOP_MULT_TRANSPOSE_ADD=6)
335:       parameter(MATOP_SOLVE=7)
336:       parameter(MATOP_SOLVE_ADD=8)
337:       parameter(MATOP_SOLVE_TRANSPOSE=9)
338:       parameter(MATOP_SOLVE_TRANSPOSE_ADD=10)
339:       parameter(MATOP_LUFACTOR=11)
340:       parameter(MATOP_CHOLESKYFACTOR=12)
341:       parameter(MATOP_RELAX=13)
342:       parameter(MATOP_TRANSPOSE=14)
343:       parameter(MATOP_GETINFO=15)
344:       parameter(MATOP_EQUAL=16)
345:       parameter(MATOP_GET_DIAGONAL=17)
346:       parameter(MATOP_DIAGONAL_SCALE=18)
347:       parameter(MATOP_NORM=19)
348:       parameter(MATOP_ASSEMBLY_BEGIN=20)
349:       parameter(MATOP_ASSEMBLY_END=21)
350:       parameter(MATOP_COMPRESS=22)
351:       parameter(MATOP_SET_OPTION=23)
352:       parameter(MATOP_ZERO_ENTRIES=24)
353:       parameter(MATOP_ZERO_ROWS=25)
354:       parameter(MATOP_LUFACTOR_SYMBOLIC=26)
355:       parameter(MATOP_LUFACTOR_NUMERIC=27)
356:       parameter(MATOP_CHOLESKY_FACTOR_SYMBOLIC=28)
357:       parameter(MATOP_CHOLESKY_FACTOR_NUMERIC=29)
358:       parameter(MATOP_GET_SIZE=30)
359:       parameter(MATOP_GET_LOCAL_SIZE=31)
360:       parameter(MATOP_GET_OWNERSHIP_RANGE=32)
361:       parameter(MATOP_ILUFACTOR_SYMBOLIC=33)
362:       parameter(MATOP_ICCFACTOR_SYMBOLIC=34)
363:       parameter(MATOP_GET_ARRAY=35)
364:       parameter(MATOP_RESTORE_ARRAY=36)

366:       parameter(MATOP_CONVERT_SAME_TYPE=37)
367:       parameter(MATOP_FORWARD_SOLVE=38)
368:       parameter(MATOP_BACKWARD_SOLVE=39)
369:       parameter(MATOP_ILUFACTOR=40)
370:       parameter(MATOP_ICCFACTOR=41)
371:       parameter(MATOP_AXPY=42)
372:       parameter(MATOP_GET_SUBMATRICES=43)
373:       parameter(MATOP_INCREASE_OVERLAP=44)
374:       parameter(MATOP_GET_VALUES=45)
375:       parameter(MATOP_COPY=46)
376:       parameter(MATOP_PRINT_HELP=47)
377:       parameter(MATOP_SCALE=48)
378:       parameter(MATOP_SHIFT=49)
379:       parameter(MATOP_DIAGONAL_SHIFT=50)
380:       parameter(MATOP_ILUDT_FACTOR=51)
381:       parameter(MATOP_GET_BLOCK_SIZE=52)

383:       parameter(MATOP_GET_ROW_IJ=53)
384:       parameter(MATOP_RESTORE_ROW_IJ=54)
385:       parameter(MATOP_GET_COLUMN_IJ=55)
386:       parameter(MATOP_RESTORE_COLUMN_IJ=56)
387:       parameter(MATOP_FDCOLORING_CREATE=57)
388:       parameter(MATOP_COLORING_PATCH=58)
389:       parameter(MATOP_SET_UNFACTORED=59)
390:       parameter(MATOP_PERMUTE=60)
391:       parameter(MATOP_SET_VALUES_BLOCKED=61)


394:       parameter(MATOP_DESTROY=250)
395:       parameter(MATOP_VIEW=251)
396: !
397: !
398: !
399:       integer MATRIX_BINARY_FORMAT_DENSE
400:       parameter (MATRIX_BINARY_FORMAT_DENSE=-1)
401: !
402: !  End of Fortran include file for the Mat package in PETSc
403: !
404: #endif