Next: Exact solutions for Up: Actual program operation Previous: Eigenstrain programs
The finite difference electrical conductivity programs can be operated quite similarly to the finite element electrical conductivity programs, with some small differences given below.
(1) Since the finite difference programs only handle diagonal conductivity tensors, the local phase conductivity variable, sigma, is not a full tensor but only has space for the diagonal elements.
(2) The subroutine DEMBX is not exited until the gtest criterion is met. The total number of conjugate gradient cycles possible within DEMBX is given by ncgsteps, a parameter that is set by the user. Intermediate results are printed out every ncheck cycles, from within DEMBX.
(3) The main variables, like the voltage vector u, are dimensioned ns2 = (nx + 2) x (ny + 2) x (nz + 2), rather than ns = nx x ny x nz like in the finite element programs. Also, these dimensions occur explicitly only in the main program, so only have to be changed there.
There are other differences, but these are not in sections that must be changed by the user. Sections 2 and 3 of this manual and the comments in the actual programs should be consulted if the user wishes to know more about these differences.