A key for reusable parallel linear algebra software

摘要

We propose an object-oriented design which enables very good code reuse for both se- quential and parallel linear algebra applications. A linear algebra class library called LAKe is implemented using our design method. We introduce a new reuse mechanism called matrix shape which enables us to derive the implementation of both the sequential and the parallel version of the iterative methods of Linear Algebra Kernels (LAKe). We show that poly- morphism is insufficient to achieve our goal and that both genericity and polymorphism are needed. We propose a new design pattern as a part of the solution. Some numerical experi- ments validate our approach and show that efficiency is not sacrificed.