%                                                                              
%FILE  s3dkq4m2.mtx                                                            
%TITLE Cyl shell R/t=1000 unif 150x100 quad mesh DKQ elem with drill rot       
%KEY   s3dkq4m2                                                                
%                                                                              
%                                                                              
%CONTRIBUTOR Reijo Kouhia (reijo.kouhia@hut.fi)                                
%                                                                              
%REFERENCE   M. Benzi, R. Kouhia, M.Tuma: An assesment of some                 
%            preconditioning techniques in shell problems                      
%            Technical Report LA-UR-97-3892, Los Alamos National Laboratory    
%                                                                              
%BEGIN DESCRIPTION                                                             
% Matrix from a static analysis of a cylindrical shell                         
% Radius to thickness ratio R/t = 1000                                         
% Length to radius ratio    R/L = 1                                            
% One octant discretized with uniform 150 x 100 quadrilateral mesh             
% element:                                                                     
% facet-type shell element where the bending part is formulated                
% using the discrete Kirchhoff quadrilateral (Lyons-Crisfield version),        
% the membrane part includes drilling rotations using                          
% the Hughes-Brezzi formulation with (regularizing parameter = G/1000,         
% where G is the shear modulus)                                                
% full 2x2 Gauss-Legendre integration                                          
% --------------------------------------------------------------------------   
% Note:                                                                        
% The sparsity pattern of the matrix is determined from the element            
% connectivity data assuming that the element matrix is full.                  
% Since this case the  material model is linear isotropically elastic          
% and the FE mesh is  uniform there exist some zeros.                          
% Since the removal of those zero elements is trivial                          
% but the reconstruction of the current sparsity                               
% pattern is impossible from the sparsified structure without any further      
% knowledge of the element connectivity, the zeros are retained in this file.  
% ---------------------------------------------------------------------------  
%END DESCRIPTION                                                               
%                                                                              
%