← Back to team overview

dolfin team mailing list archive

Re: Re: Linear Algebra

 



One example is the module elasticity-updated in DOLFIN for which
assembly time dominates solution time (using the old form evaluation
system in DOLFIN). This would be a typical example of a system which
is more complicated than Poisson and where assembly time may be
significant.

This is a good point. There are these complex-model linear problems
where you just form the matrix and solve. Is that what you were
thinking? Since you can vary the elasticity tensor these days with
special materials, one could even imagine optimizing over varying
elasticity tensors, in an extreme case.

                                Ridg

Yes. In this case, you form the matrix, solve, and modify (move) the
mesh in each time step. The domininating cost is the assembly.

/Anders

From my experience, solving time-dependent Navier-Stokes in 3d (linear
elements on tetrahedrons), using PetSc both for inserting elements (as blocks) in assembling the "momentum matrix" and solving the resulting linear system with preconditioned GMRES, the amount of time spent on one assemble of the matrix can be roughly about 3 times the time for solving the system.

/Johan





Follow ups

References