https://github.com/basp1/spmat выложил тут мелкую библиотеку для работы с разреженными матрицами. C++. простая арифметика, LLt, LDLt, rank1 update, solve, если вы понимаете, о чем я :)

Предположим, у нас есть такой класс:

struct A
{
string *ps;
A(const string &s = string()): ps(new string(s)) { }
~A() { delete ps; }
};
( Конструкторы копирования )