Our FEA code has its own local serial solver which we use at the final step to solve the subspace etc. The serial solver is both direct and iterative (one chooses). One could just as easily use a PETSc solver -- serial or parallel.
We did it this way since we had a legacy code that already had the algorithms programed up. To convert to parallel we only had to program the operator and inverse operator applications which we did with PETSc. On 3/17/13 12:40 PM, Sonya Blade wrote: > >> For our FEA code we have used PETSc to solve mechanical eigen problems using >> subspace as well as parallel-ARPACK. It works well and is not that hard >> to program. > Which package exactly solved the subspace or you used them in conjunction? > Can you give pseudo-code for that, for example do I have to call a set of > functions or only the one function is associated with subspace iteration in > Petsc or Arpack? > > How about the other methods Jacobi, simultaneous iteration etc.. is it > possible > to use them in Petsc? > >> but it is minor for your average mechanical problem. > If you call 5000x5000 matrix small I really wonder what you will call huge ? > > Regards, > > > -- ----------------------------------------------- Sanjay Govindjee, PhD, PE Professor of Civil Engineering Vice Chair for Academic Affairs 779 Davis Hall Structural Engineering, Mechanics and Materials Department of Civil Engineering University of California Berkeley, CA 94720-1710 Voice: +1 510 642 6060 FAX: +1 510 643 5264 s_g at berkeley.edu http://www.ce.berkeley.edu/~sanjay ----------------------------------------------- New Books: Engineering Mechanics of Deformable Solids: A Presentation with Exercises http://www.oup.com/us/catalog/general/subject/Physics/MaterialsScience/?view=usa&ci=9780199651641 http://ukcatalogue.oup.com/product/9780199651641.do http://amzn.com/0199651647 Engineering Mechanics 3 (Dynamics) http://www.springer.com/materials/mechanics/book/978-3-642-14018-1 http://amzn.com/3642140181 -----------------------------------------------
