Let $\mathsf{R}$ be a PID and consider a collection of free, finitely generated $\mathsf{R}$-modules $V\_1,\ldots,V\_n$ along with module maps $m\_j:V\_j \to V\_{j+1}$. That is, we have the following sequence 


$$ V\_1 \stackrel{m\_1}{\to} V\_2 \stackrel{m\_2}{\to} \ldots \stackrel{m\_{n-1}}{\to} V\_n $$


In particular, if we fix bases of each $V\_j$ a priori, we may view each $m\_j$ as a matrix with entries in $\mathsf{R}$. Recall that the [Smith Normal Form](http://en.wikipedia.org/wiki/Smith_normal_form) (SNF) of such a matrix is obtained by performing row and column operations - that is, change of bases of $V\_j$ and $V\_{j+1}$- so that the end result is a canonical matrix with zero off-diagonal entries and where each diagonal entry divides the next.


Start at the end of this sequence and note that we can always put the last matrix $m\_{n-1}$ in SNF by performing suitable row and column operations. But now, if we try to do the same with $m\_{n-2}$ for instance, we would generically have to change the basis of $V\_{n-1}$ and in this new basis $m\_{n-1}$ may no longer be in SNF.


Here are the questions:



> 
> Under what conditions can one prove the existence of bases of the $V\_j$s such that each $m\_j$ is in SNF with respect to the chosen basis of $V\_j$ and $V\_{j+1}$?
> 
> 
> 


And when such conditions hold,



> 
> Is there an efficient algorithm for finding these bases?
> 
> 
> 



Simultaneous Smith Normalization of a Composable Matrix Sequence
The Smith Normal Form of a product of matrices is not, in general, the product of the Smith Normal Forms. So most of the time you can't do this. 


There are some results in this direction in the book "Integral Matrices" by M. Newman, which I have found useful in the past. However, none of these results are constructive.


EDIT: 
Let me add one such condition. If $\det(m\_1\cdots m\_{j-1})$ is relatively prime to $\det(m\_j)$ for $j=2,\ldots,n$, then the SNF of $m\_1\cdots m\_n$ is the product of the Smith normal forms of $m\_1,\ldots,m\_n$. 


