best way to model matrices whose dimensions increase dynamically in CUDA 4.0

I need to solve problems that manipulate matrices whose dimensions have to increase dynamically,
as the algorithm advances.

What is the best way to model such matrices in CUDA 4.0?

And where can one find eaxmples of such code for CUDA 4.0?