#yspWGvfmhvI.m
function spvF=yspWGvfmhvI(baasi0,l,GIt,EIw)
# sparse fransfer matrix of the of the first order analysis 
# for a beam warping (Sign Convention 2)
# Zp=F*Zv, Zp'=[Wp Fip Qp Mp]
# 
#  l - length of the element
#  GIt - torsional stiffness
#  EIw - warping stiffness
#  baasi0 - scaling multiplier for the displacements (i = GIt)
#
#
##=========================================================================
## This Program is written by Andres Lahe,   2015-03-30
##                    e-mail: andres.lahe@ttu.ee
## LAST MODIFIED: Andres Lahe,    2015-03-30
## Copyright (c)  2015  by Tallinn University of Technology 
##                Department of Mechanics
##                http://www.ttu.ee/
##
## This program is free software; you can redistribute it and/or
## modify it under the terms of the GNU General Public License
## as published by the Free Software Foundation; either version 2
## of the License, or (at your option) any later version.

## This program is distributed in the hope that it will be useful,
## but WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
## GNU General Public License for more details.

## You should have received a copy of the GNU General Public License
## along with this program; if not, write to the Free Software
## Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
## http://www.gnu.org/licenses/old-licenses/gpl-2.0.html
##=========================================================================
#
#%%%%%%%%%%%%%%%%%
#if nargin != 4
if ~(nargin==4)
error(' function yspWGvfmhvI have wrong number of input arguments!') 
end
#%%%%%%%%%%%%%%%%%
#
#
x=l;
spvF=sparse(4,8);
spvI=speye(4);
#spvI1=spvI.*(-1);
#spvF=spInsertBtoA(spvF,1,1,spvI);
spvU=yspWGfhlin(baasi0,l,x,GIt,EIw);
#spvU=spvU1.*(-1);
#spvF=spInsertBtoA(spvF,1,5,spvU);
#
spvF=spInsertBtoA(spvF,1,1,spvU);
spvI1=spvI.*(-1);
spvF=spInsertBtoA(spvF,1,5,spvI1);
#
endfunction

