// =============================================================================
// Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
// Copyright (C) ????-2008 - INRIA
//
//  This file is distributed under the same license as the Scilab package.
// =============================================================================
// <-- CLI SHELL MODE -->
mode(5)
lines(0)
ij=[1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6];
v=[1;2;3;4;5;6;7;8;9];
ij1=[1 6;1 2;6 5];
v1=[-1;-1;-1];
vc=v+%i*(21:29)';
v1c=v1+%i*[0;0.3;-1.2];
zer=sparse([],[],[6,6]);
//
a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
b=sparse([1 6;1 2;6 5],[-1;-1;-1],[6 6]);
ma=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],-(1:9),[6 6]);
mb=sparse([1 6;1 2;6 5],-[-1;-1;-1],[6 6]);
ac=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],(1:9)+%i*(21:29),[6 6]);
bc=sparse([1 6;1 2;6 5],[-1;-1;-1]+%i*[0;0.3;-1.2],[6 6]);
mac=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],-(1:9)-%i*(21:29),[6 6]);
mbc=sparse([1 6;1 2;6 5],[1;1;1]+%i*[0;-0.3;1.2],[6 6]);
//
if %t then
//-------------------------------------------------------------
//test des primitives sparse, full
//--------------------------------------------------------------
//  -- for real matrix
if or(full(sparse(0.3))<>0.3) then bugmes();quit;end
v=0*ones(1,3);v(3)=5;v(1)=1;
if or(full(sparse([1 1;1 3],[1 5]))<>v) then bugmes();quit;end
v=0*ones(1,300);v(300)=5;v(1)=1;
if or(full(sparse([1 1;1 300],[1 5]))<>v) then bugmes();quit;end
v=0*ones(3,1);v(3)=5;v(1)=1;
if or(full(sparse([1 1;3 1],[1 5]))<>v) then bugmes();quit;end
v=0*ones(300,1);v(300)=5;v(1)=1;
if or(full(sparse([1 1;300 1],[1 5]))<>v) then bugmes();quit;end
sp=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
a=[0,0,3,0,2,1;;5,0,0,4,0,0;0,0,0,0,0,0;0,0,7,6,0,0;8,0,0,0,0,0;0,0,0,0,0,9];
if or(full(sp)<>a) then bugmes();quit;end
sp=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[8 6]);a(8,6)=0;
if or(full(sp)<>a) then bugmes();quit;end
if or(full(sparse([],[],[4 10]))<>0*ones(4,10)) then bugmes();quit;end
v=sparse([2 1;3 1;4 1;6 1],[10;11;12;13],[6,1]);
a=[0;10;11;12;0;13];
if or(full(v)<>a) then bugmes();quit;end
v=sparse([1 2;1 3;1 4;1 6],[10;11;12;13],[1,6]);
a=[0,10,11,12,0,13];
if or(full(v)<>a) then bugmes();quit;end
//  -- for complex matrix
sp=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],(1:9)+%i*(21:29),[6 6]);
a=[0,0,3+%i*23,0,2+%i*22,1+%i*21
   5+%i*25,0,0,4+%i*24,0,0
   0,0,0,0,0,0
   0,0,7+%i*27,6+%i*26,0,0
   8+%i*28,0,0,0,0,0
   0,0,0,0,0,9+%i*29];
if or(full(sp)<>a) then bugmes();quit;end
sp=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],(1:9)+%i*(21:29),[8 6]);
a(8,6)=0;
if or(full(sp)<>a) then bugmes();quit;end
v=sparse([2 1;3 1;4 1;6 1],[10-3*%i;11;12+5*%i;13+0.5*%i],[6,1]);
a=[0;10-%i*3;11;12+%i*5;0;13+%i*0.5];
if or(full(v)<>a) then bugmes();quit;end
v=sparse([1 2;1 3;1 4;1 6],[10-3*%i;11;12+5*%i;13+0.5*%i],[1,6]);
a=[0,10-%i*3,11,12+%i*5,0,13+%i*0.5];
if or(full(v)<>a) then bugmes();quit;end
//-------------------------------------------------------------
//test de la primitives spget
//--------------------------------------------------------------
sp=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
[ij,v]=spget(sp);
if or(ij<>[1,3;1,5;1,6;2,1;2,4;4,3;4,4;5,1;6,6])|or(v<>[3;2;1;5;4;7;6;8;9]) then bugmes();quit;end
zer=sparse([],[],[6,6]);[ij,v]=spget(zer);
if ij<>[]|v<>[] then bugmes();quit;end
//-----------------------------------------------
// addition and subtraction tests
//-----------------------------------------------
// -- real real
a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
b=sparse([1 6;1 2;6 5],[-1;-1;-1],[6 6]);
v=sparse([2 1;3 1;4 1;6 1],[10;11;12;13],[6,1]);
vt=sparse([1 2;1 3;1 4;1 6],[10;11;12;13],[1,6]);
ma=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],-(1:9),[6 6]);
mb=sparse([1 6;1 2;6 5],-[-1;-1;-1],[6 6]);
zer=sparse([],[],[6,6]);
ac=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],(1:9)+%i*(21:29),[6 6]);
bc=sparse([1 6;1 2;6 5],[-1;-1;-1]+%i*[0;0.3;-1.2],[6 6]);
vc=sparse([2 1;3 1;4 1;6 1],[10-3*%i;11;12+5*%i;13+0.5*%i],[6,1]);
vct=sparse([1 2;1 3;1 4;1 6],[10-3*%i;11;12+5*%i;13+0.5*%i],[1,6]);
mac=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],-(1:9)-%i*(21:29),[6 6]);
mbc=sparse([1 6;1 2;6 5],[1;1;1]+%i*[0;-0.3;1.2],[6 6]);
if or(full(v+sparse([],[],[6,1]))<>full(v)) then bugmes();quit;end
if or(full(sparse([],[],[6,1])+v)<>full(v)) then bugmes();quit;end
if or(full(v+v)<>full(v)+full(v)) then bugmes();quit;end
if or(full(v-v)<>full(v)-full(v)) then bugmes();quit;end
if or(full(vt+sparse([],[],[1,6]))<>full(vt)) then bugmes();quit;end
if or(full(vt+vt)<>full(vt)+full(vt)) then bugmes();quit;end
if or(full(vt-vt)<>full(vt)-full(vt)) then bugmes();quit;end
if or(full(zer+zer)<>0*ones(6,6)) then bugmes();quit;end
if or(full(a+a)<>full(a)+full(a)) then bugmes();quit;end
if or(full(b+b)<>full(b)+full(b)) then bugmes();quit;end
if or(full(a+zer)<>full(a)) then bugmes();quit;end
if or(full(zer+a)<>full(a)) then bugmes();quit;end
if or(full(b+a)<>full(b)+full(a)) then bugmes();quit;end
if or(full(a+b)<>full(b)+full(a)) then bugmes();quit;end
if or(full(a+ma)<>full(a)+full(ma)) then bugmes();quit;end
if or(full(a-a)<>full(a)-full(a)) then bugmes();quit;end
if or(full(a-ma)<>full(a)-full(ma)) then bugmes();quit;end
if or(full(b-mb)<>full(b)-full(mb)) then bugmes();quit;end
if or(full(a-zer)<>full(a)) then bugmes();quit;end
if or(full(zer-a)<>-full(a)) then bugmes();quit;end
if or(full(a-mb)<>full(a)-full(mb)) then bugmes();quit;end
//  -- real complex
if or(full(sparse([],[],[6,1])+vc)<>full(vc)) then bugmes();quit;end
if or(full(v+vc)<>full(v)+full(vc)) then bugmes();quit;end
if or(full(v-vc)<>full(v)-full(vc)) then bugmes();quit;end
if or(full(vt+vct)<>full(vt)+full(vct)) then bugmes();quit;end
if or(full(vt-vct)<>full(vt)-full(vct)) then bugmes();quit;end
if or(full(a+ac)<>full(a)+full(ac)) then bugmes();quit;end
if or(full(b+bc)<>full(b)+full(bc)) then bugmes();quit;end
if or(full(a+bc)<>full(a)+full(bc)) then bugmes();quit;end
if or(full(zer+ac)<>full(zer)+full(ac)) then bugmes();quit;end
if or(full(b+ac)<>full(b)+full(ac)) then bugmes();quit;end
if or(full(a-ac)<>full(a)-full(ac)) then bugmes();quit;end
if or(full(b-bc)<>full(b)-full(bc)) then bugmes();quit;end
if or(full(a-bc)<>full(a)-full(bc)) then bugmes();quit;end
if or(full(zer-ac)<>full(zer)-full(ac)) then bugmes();quit;end
if or(full(b-ac)<>full(b)-full(ac)) then bugmes();quit;end
// -- complex real
if or(full(vc+v)<>full(vc)+full(v)) then bugmes();quit;end
if or(full(vc-v)<>full(vc)-full(v)) then bugmes();quit;end
if or(full(vct+vt)<>full(vct)+full(vt)) then bugmes();quit;end
if or(full(vct-vt)<>full(vct)-full(vt)) then bugmes();quit;end
if or(full(ac+a)<>full(ac)+full(a)) then bugmes();quit;end
if or(full(bc+b)<>full(bc)+full(b)) then bugmes();quit;end
if or(full(ac+b)<>full(ac)+full(b)) then bugmes();quit;end
if or(full(ac+zer)<>full(zer)+full(ac)) then bugmes();quit;end
if or(full(bc+a)<>full(bc)+full(a)) then bugmes();quit;end
if or(full(ac-a)<>full(ac)-full(a)) then bugmes();quit;end
if or(full(bc-b)<>full(bc)-full(b)) then bugmes();quit;end
if or(full(ac-b)<>full(ac)-full(b)) then bugmes();quit;end
if or(full(ac-zer)<>-full(zer)+full(ac)) then bugmes();quit;end
if or(full(bc-a)<>full(bc)-full(a)) then bugmes();quit;end
// -- complex complex
if or(full(vc+vc)<>full(vc)+full(vc)) then bugmes();quit;end
if or(full(vc-vc)<>full(vc)-full(vc)) then bugmes();quit;end
if or(full(vct+vct)<>full(vct)+full(vct)) then bugmes();quit;end
if or(full(vct-vct)<>full(vct)-full(vct)) then bugmes();quit;end
if or(full(ac+ac)<>full(ac)+full(ac)) then bugmes();quit;end
if or(full(bc+bc)<>full(bc)+full(bc)) then bugmes();quit;end
if or(full(ac+bc)<>full(ac)+full(bc)) then bugmes();quit;end
if or(full(bc+ac)<>full(bc)+full(ac)) then bugmes();quit;end
if or(full(ac-ac)<>full(zer)) then bugmes();quit;end
if or(full(bc-bc)<>full(zer)) then bugmes();quit;end
if or(full(ac-bc)<>full(ac)-full(bc)) then bugmes();quit;end
if or(full(bc-ac)<>full(bc)-full(ac)) then bugmes();quit;end
// sparse full and full sparse operation (soft coded apoerations)
if or(full(a+1)<>full(a)+1) then bugmes();quit;end
if or(full(1+a)<>full(a)+1) then bugmes();quit;end
if or(full(a+ones(6,6))<>full(a)+ones(6,6)) then bugmes();quit;end
if or(full(ones(6,6)+a)<>full(a)+ones(6,6)) then bugmes();quit;end
if or(full(a+2*eye())<>full(a)+2*eye()) then bugmes();quit;end
if or(full(2*eye()+a)<>full(a)+2*eye()) then bugmes();quit;end
if or(full(a-1)<>full(a)-1) then bugmes();quit;end
if or(full(1-a)<>1-full(a)) then bugmes();quit;end
if or(full(a-ones(6,6))<>full(a)-ones(6,6)) then bugmes();quit;end
if or(full(ones(6,6)-a)<>ones(6,6)-full(a)) then bugmes();quit;end
if or(full(a-2*eye())<>full(a)-2*eye()) then bugmes();quit;end
if or(full(2*eye()-a)<>2*eye()-full(a)) then bugmes();quit;end
if or(full(ac+1)<>full(ac)+1) then bugmes();quit;end
if or(full(1+ac)<>full(ac)+1) then bugmes();quit;end
if or(full(ac+ones(6,6))<>full(ac)+ones(6,6)) then bugmes();quit;end
if or(full(ones(6,6)+ac)<>full(ac)+ones(6,6)) then bugmes();quit;end
if or(full(ac+2*eye())<>full(ac)+2*eye()) then bugmes();quit;end
if or(full(2*eye()+ac)<>full(ac)+2*eye()) then bugmes();quit;end
if or(full(ac-1)<>full(ac)-1) then bugmes();quit;end
if or(full(1-ac)<>1-full(ac)) then bugmes();quit;end
if or(full(ac-ones(6,6))<>full(ac)-ones(6,6)) then bugmes();quit;end
if or(full(ones(6,6)-ac)<>ones(6,6)-full(ac)) then bugmes();quit;end
if or(full(ac-2*eye())<>full(ac)-2*eye()) then bugmes();quit;end
if or(full(2*eye()-ac)<>2*eye()-full(ac)) then bugmes();quit;end
if or(full(ac+full(bc))<>full(ac)+full(bc)) then bugmes();quit;end
if or(full(ac-full(bc))<>full(ac)-full(bc)) then bugmes();quit;end
if or(full(full(ac)+full(bc))<>full(ac)+full(bc)) then bugmes();quit;end
end
if %t then
//-----------------------------------------------
// multiplication tests
//-----------------------------------------------
// real x real
// sparse scalar , saclar sparse
if or(full(a*2)<>full(a)*2) then bugmes();quit;end
if or(full(2*a)<>full(a)*2) then bugmes();quit;end
if  a*[]<>[] then bugmes();quit;end
if  []*a<>[] then bugmes();quit;end
c=rand(6,6);
if norm(a*c-full(a)*c) >100*%eps then bugmes();quit;end
if norm(c*a-c*full(a)) >100*%eps then bugmes();quit;end
// sparse sparse
if or(full(zer*zer)<>full(zer)) then bugmes();quit;end
if or(full(a*zer)<>full(zer)) then bugmes();quit;end
if or(full(zer*a)<>full(zer)) then bugmes();quit;end
if  norm(full(a*a)-full(a)*full(a),1)>100*%eps then bugmes();quit;end
if  norm(full(b*b)-full(b)*full(b),1)>100*%eps then bugmes();quit;end
if  norm(full(a*b)-full(a)*full(b),1)>100*%eps then bugmes();quit;end
if  norm(full(b*a)-full(b)*full(a),1)>100*%eps then bugmes();quit;end
// complex x real real x complex
// sparse scalar , scalar sparse
if or(full(ac*2)<>full(ac)*2) then bugmes();quit;end
if or(full(2*ac)<>full(ac)*2) then bugmes();quit;end
if norm(full((2+%i)*a)-(2+%i)*full(a),1) >100*%eps then bugmes();quit;end
if norm(full(a*(2+%i))-(2+%i)*full(a),1) >100*%eps then bugmes();quit;end
if  ac*[]<>[] then bugmes();quit;end
if  []*ac<>[] then bugmes();quit;end
c=rand(6,6);
cc=c+rand(6,6)*%i;
if norm(ac*c-full(ac)*c) >100*%eps then bugmes();quit;end
if norm(cc*a-cc*full(a)) >100*%eps then bugmes();quit;end
// sparse sparse
if or(full(ac*zer)<>full(zer)) then bugmes();quit;end
if or(full(zer*ac)<>full(zer)) then bugmes();quit;end
if  norm(full(ac*a)-full(ac)*full(a),1)>100*%eps then bugmes();quit;end
if  norm(full(bc*b)-full(bc)*full(b),1)>100*%eps then bugmes();quit;end
if  norm(full(a*bc)-full(a)*full(bc),1)>100*%eps then bugmes();quit;end
if  norm(full(b*ac)-full(b)*full(ac),1)>100*%eps then bugmes();quit;end
// // complex x complex
if norm(ac*cc-full(ac)*cc) >100*%eps then bugmes();quit;end
if norm(cc*ac-cc*full(ac)) >100*%eps then bugmes();quit;end
// sparse sparse
if  norm(full(ac*ac)-full(ac)*full(ac),1)>100*%eps then bugmes();quit;end
if  norm(full(bc*bc)-full(bc)*full(bc),1)>100*%eps then bugmes();quit;end
if  norm(full(bc*ac)-full(bc)*full(ac),1)>100*%eps then bugmes();quit;end
//----------------------------------------------------------
// element wise multiplication tests
//----------------------------------------------------------
if or(full(ac.*2)<>full(ac)*2) then bugmes();quit;end
if or(full((2).*ac)<>full(ac)*2) then bugmes();quit;end
if  a.*[]<>[] then bugmes();quit;end
if  [].*a<>[] then bugmes();quit;end
c=rand(6,6);
//if norm(a*c-full(a)*c) >100*%eps then bugmes();quit;end
//if norm(c*a-c*full(a)) >100*%eps then bugmes();quit;end
// sparse sparse
if or(full(zer.*zer)<>full(zer)) then bugmes();quit;end
if or(full(a.*zer)<>full(zer)) then bugmes();quit;end
if or(full(zer.*a)<>full(zer)) then bugmes();quit;end
if  norm(full(a.*a)-full(a).*full(a),1)>100*%eps then bugmes();quit;end
if  norm(full(b.*b)-full(b).*full(b),1)>100*%eps then bugmes();quit;end
if  norm(full(a.*b)-full(a).*full(b),1)>100*%eps then bugmes();quit;end
if  norm(full(b.*a)-full(b).*full(a),1)>100*%eps then bugmes();quit;end
// complex x real real x complex
// sparse scalar , scalar sparse
if norm(full((2+%i).*a)-(2+%i).*full(a),1) >100*%eps then bugmes();quit;end
if norm(full(a.*(2+%i))-(2+%i).*full(a),1) >100*%eps then bugmes();quit;end
if  ac.*[]<>[] then bugmes();quit;end
if  [].*ac<>[] then bugmes();quit;end
c=rand(6,6);
cc=c+rand(6,6)*%i;
if norm(full(ac.*c)-full(ac).*c) >100*%eps then bugmes();quit;end
if norm(full(cc.*a)-cc.*full(a)) >100*%eps then bugmes();quit;end
// sparse sparse
if or(full(ac.*zer)<>full(zer)) then bugmes();quit;end
if or(full(zer.*ac)<>full(zer)) then bugmes();quit;end
if  norm(full(ac.*a)-full(ac).*full(a),1)>100*%eps then bugmes();quit;end
if  norm(full(bc.*b)-full(bc).*full(b),1)>100*%eps then bugmes();quit;end
if  norm(full(a.*bc)-full(a).*full(bc),1)>100*%eps then bugmes();quit;end
if  norm(full(b.*ac)-full(b).*full(ac),1)>100*%eps then bugmes();quit;end
// // complex x complex
if norm(ac.*cc-full(ac).*cc) >100*%eps then bugmes();quit;end
if norm(cc.*ac-cc.*full(ac)) >100*%eps then bugmes();quit;end
// sparse sparse
if  norm(full(ac.*ac)-full(ac).*full(ac),1)>100*%eps then bugmes();quit;end
if  norm(full(bc.*bc)-full(bc).*full(bc),1)>100*%eps then bugmes();quit;end
if  norm(full(bc.*ac)-full(bc).*full(ac),1)>100*%eps then bugmes();quit;end
// ----------------------------------------------------------
// test de la transposition
//-----------------------------------------------------------
if or(full(a')<>full(a)') then bugmes();quit;end
if or(full(ac')<>full(ac)') then bugmes();quit;end
if or(full(zer')<>full(zer)' ) then bugmes();quit;end
v=sparse([2 1;3 1;4 1;6 1],[10;11;12;13],[6,1]);
if or(full(v')<>full(v)') then bugmes();quit;end
if or(full((v')')<>full(v)) then bugmes();quit;end
vc=sparse([2 1;3 1;4 1;6 1],[10-3*%i;11;12+5*%i;13+0.5*%i],[6,1]);
if or(full(vc')<>full(vc)') then bugmes();quit;end
if or(full((vc')')<>full(vc)) then bugmes();quit;end
// ----------------------------------------------------------
// test des concatenation
//-----------------------------------------------------------
if or(full([a])<>full(a)) then bugmes();quit;end
if or(full([a b])<>[full(a) full(b)]) then bugmes();quit;end
if or(full([a;b])<>[full(a);full(b)]) then bugmes();quit;end
if or(full([a []])<>full(a)) then bugmes();quit;end
if or(full([a;[]])<>full(a)) then bugmes();quit;end
if or(full([a zer])<>[full(a) full(zer)]) then bugmes();quit;end
if or(full([zer;b])<>[full(zer);full(b)]) then bugmes();quit;end
if or(full([ac])<>full(ac)) then bugmes();quit;end
if or(full([ac b])<>[full(ac) full(b)]) then bugmes();quit;end
if or(full([ac;b])<>[full(ac);full(b)]) then bugmes();quit;end
if or(full([ac []])<>full(ac)) then bugmes();quit;end
if or(full([ac;[]])<>full(ac)) then bugmes();quit;end
if or(full([a bc])<>[full(a) full(bc)]) then bugmes();quit;end
if or(full([a;bc])<>[full(a);full(bc)]) then bugmes();quit;end
if or(full([ac bc])<>[full(ac) full(bc)]) then bugmes();quit;end
if or(full([ac;bc])<>[full(ac);full(bc)]) then bugmes();quit;end
// ----------------------------------------------------------
// test des extractions
//-----------------------------------------------------------
af=full(a);
if or(full(a(1,3))<>af(1,3)) then bugmes();quit;end
if or(full(a(1,4))<>af(1,4)) then bugmes();quit;end
if or(full(a(1,:))<>af(1,:)) then bugmes();quit;end
v=sparse([2 1;3 1;4 1;6 1],[10;11;12;13],[6,1]);vf=full(v);
if or(full(v(:))<>vf(:)) then bugmes();quit;end
if or(full(v(3:4))<>vf(3:4)) then bugmes();quit;end
if or(full(v([1 5]))<>vf([1 5])) then bugmes();quit;end
if or(full(v([4 3]))<>vf([4 3])) then bugmes();quit;end
if or(full(v([4 4]))<>vf([4 4])) then bugmes();quit;end
if or(full(v([1 1]))<>vf([1 1])) then bugmes();quit;end
v=v';vf=vf';
if or(full(v(:))<>vf(:)) then bugmes();quit;end
if or(full(v(3:4))<>vf(3:4)) then bugmes();quit;end
if or(full(v([1 5]))<>vf([1 5])) then bugmes();quit;end
if or(full(v([4 3]))<>vf([4 3])) then bugmes();quit;end
if or(full(v([4 4]))<>vf([4 4])) then bugmes();quit;end
if or(full(v([1 1]))<>vf([1 1])) then bugmes();quit;end
acff=full(ac);
if or(full(ac(1,3))<>acff(1,3)) then bugmes();quit;end
if or(full(ac(1,4))<>acff(1,4)) then bugmes();quit;end
if or(full(ac(1,:))<>acff(1,:)) then bugmes();quit;end
vc=sparse([2 1;3 1;4 1;6 1],[10-3*%i;11;12+5*%i;13+0.5*%i],[6,1]);vcf=full(vc);
if or(full(vc(:))<>vcf(:)) then bugmes();quit;end
if or(full(vc(3:4))<>vcf(3:4)) then bugmes();quit;end
if or(full(vc([1 5]))<>vcf([1 5])) then bugmes();quit;end
if or(full(vc([4 3]))<>vcf([4 3])) then bugmes();quit;end
if or(full(vc([4 4]))<>vcf([4 4])) then bugmes();quit;end
if or(full(vc([1 1]))<>vcf([1 1])) then bugmes();quit;end
vc=vc';vcf=vcf';
if or(full(vc(:))<>vcf(:)) then bugmes();quit;end
if or(full(vc(3:4))<>vcf(3:4)) then bugmes();quit;end
if or(full(vc([1 5]))<>vcf([1 5])) then bugmes();quit;end
if or(full(vc([4 3]))<>vcf([4 3])) then bugmes();quit;end
if or(full(vc([4 4]))<>vcf([4 4])) then bugmes();quit;end
if or(full(vc([1 1]))<>vcf([1 1])) then bugmes();quit;end
// ----------------------------------------------------------
// test des insertions
//-----------------------------------------------------------
end
a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
vt=sparse([1 2;1 3;1 4;1 6],[10;11;12;13],[1,6]);
// full line insertion
//----------------------
a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(7,:)=vt;A=full(a);Vt=full(vt);A(7,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(8,:)=vt;A=full(a);Vt=full(vt);A(8,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 3],:)=[vt;2*vt];A=full(a);Vt=full(vt);A([1 3],:)=[Vt;2*Vt];
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([3 1],:)=[vt;2*vt];A=full(a);Vt=full(vt);A([3 1],:)=[Vt;2*Vt];
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 3 1],:)=[vt;2*vt;3*vt];A=full(a);Vt=full(vt);A([1 3 1],:)=[Vt;2*Vt;3*Vt];
if or(full(a1)<>A) then bugmes();quit;end
//  insert zero vector
vt=sparse([],[],[1,6]);
a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(7,:)=vt;A=full(a);Vt=full(vt);A(7,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(8,:)=vt;A=full(a);Vt=full(vt);A(8,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 3],:)=[vt;2*vt];A=full(a);Vt=full(vt);A([1 3],:)=[Vt;2*Vt];
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([3 1],:)=[vt;2*vt];A=full(a);Vt=full(vt);A([3 1],:)=[Vt;2*Vt];
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 3 1],:)=[vt;2*vt;3*vt];A=full(a);Vt=full(vt);A([1 3 1],:)=[Vt;2*Vt;3*Vt];
if or(full(a1)<>A) then bugmes();quit;end
a=sparse([],[],[6,6]);
a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(7,:)=vt;A=full(a);Vt=full(vt);A(7,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(8,:)=vt;A=full(a);Vt=full(vt);A(8,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([1 1;1 3;1 6;2 1;2 2;2 4],10:15,[2,6]);
a1=a;a1([1 3],:)=b;A=full(a);B=full(b);A([1 3],:)=B;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([3 1],:)=b;A=full(a);B=full(b);A([3 1],:)=B;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([1 1;1 3;1 6;2 1;2 2;2 4;3 3;3 5],10:17,[3,6]);
a1=a;a1([1 3 1],:)=b;A=full(a);B=full(b);A([1 3 1],:)=B;
if or(full(a1)<>A) then bugmes();quit;end
//  insert zero vector
vt=sparse([],[],[1,6]);
a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1,:)=vt;A=full(a);Vt=full(vt);A(1,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);Vt=full(vt);A(3,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(7,:)=vt;A=full(a);Vt=full(vt);A(7,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(8,:)=vt;A=full(a);Vt=full(vt);A(8,:)=Vt;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([1 1;1 3;1 6;2 1;2 2;2 4],10:15,[2,6]);
a1=a;a1([1 3],:)=b;A=full(a);B=full(b);A([1 3],:)=B;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([3 1],:)=b;A=full(a);B=full(b);A([3 1],:)=B;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([1 1;1 3;1 6;2 1;2 2;2 4;3 3;3 5],10:17,[3,6]);
a1=a;a1([1 3 1],:)=b;A=full(a);B=full(b);A([1 3 1],:)=B;
if or(full(a1)<>A) then bugmes();quit;end
// full column insertion
//----------------------
a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
v=sparse([2 1;3 1;4 1;6 1],[10;11;12;13],[6,1]);
a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,3)=v;A=full(a);V=full(v);A(:,3)=V;
if or(full(a1)<>A) then bugmes();quit;end
//
a1=a;a1(:,7)=v;A=full(a);V=full(v);A(:,7)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,8)=v;A=full(a);V=full(v);A(:,8)=V;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([1 2;2 1;3 1;3 2;4 1;6 2],10:15,[6,2]);
a1=a;a1(:,[1 3])=b;A=full(a);B=full(b);A(:,[1 3])=B;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,[3 1])=b;A=full(a);B=full(b);A(:,[3 1])=B;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([1 2;2 1;2 3;3 1;3 2;4 1;5 3;6 2],10:17,[6,3]);
a1=a;a1(:,[1 3 1])=b;A=full(a);B=full(b);A(:,[1 3 1])=B;
if or(full(a1)<>A) then bugmes();quit;end
v=sparse([],[],[6,1]);
a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,3)=v;A=full(a);V=full(v);A(:,3)=V;
if or(full(a1)<>A) then bugmes();quit;end
//
a1=a;a1(:,7)=v;A=full(a);V=full(v);A(:,7)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,8)=v;A=full(a);V=full(v);A(:,8)=V;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([],[],[6,2]);
a1=a;a1(:,[1 3])=b;A=full(a);B=full(b);A(:,[1 3])=B;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,[3 1])=b;A=full(a);B=full(b);A(:,[3 1])=B;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([],[],[6,3]);
a1=a;a1(:,[1 3 1])=b;A=full(a);B=full(b);A(:,[1 3 1])=B;
if or(full(a1)<>A) then bugmes();quit;end
a=sparse([],[],[6,6]);
a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,3)=v;A=full(a);V=full(v);A(:,3)=V;
if or(full(a1)<>A) then bugmes();quit;end
//
a1=a;a1(:,7)=v;A=full(a);V=full(v);A(:,7)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,8)=v;A=full(a);V=full(v);A(:,8)=V;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([1 2;2 1;3 1;3 2;4 1;6 2],10:15,[6,2]);
a1=a;a1(:,[1 3])=b;A=full(a);B=full(b);A(:,[1 3])=B;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,[3 1])=b;A=full(a);B=full(b);A(:,[3 1])=B;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([1 2;2 1;2 3;3 1;3 2;4 1;5 3;6 2],10:17,[6,3]);
a1=a;a1(:,[1 3 1])=b;A=full(a);B=full(b);A(:,[1 3 1])=B;
if or(full(a1)<>A) then bugmes();quit;end
v=sparse([],[],[6,1]);
a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,3)=v;A=full(a);V=full(v);A(:,3)=V;
if or(full(a1)<>A) then bugmes();quit;end
//
a1=a;a1(:,7)=v;A=full(a);V=full(v);A(:,7)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,8)=v;A=full(a);V=full(v);A(:,8)=V;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([1 2;2 1;3 1;3 2;4 1;6 2],10:15,[6,2]);
a1=a;a1(:,[1 3])=b;A=full(a);B=full(b);A(:,[1 3])=B;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,[3 1])=b;A=full(a);B=full(b);A(:,[3 1])=B;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([1 2;2 1;2 3;3 1;3 2;4 1;5 3;6 2],10:17,[6,3]);
a1=a;a1(:,[1 3 1])=b;A=full(a);B=full(b);A(:,[1 3 1])=B;
if or(full(a1)<>A) then bugmes();quit;end
v=sparse([],[],[6,1]);
a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,1)=v;A=full(a);V=full(v);A(:,1)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);V=full(v);A(:,2)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,3)=v;A=full(a);V=full(v);A(:,3)=V;
if or(full(a1)<>A) then bugmes();quit;end
//
a1=a;a1(:,7)=v;A=full(a);V=full(v);A(:,7)=V;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,8)=v;A=full(a);V=full(v);A(:,8)=V;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([],[],[6,2]);
a1=a;a1(:,[1 3])=b;A=full(a);B=full(b);A(:,[1 3])=B;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,[3 1])=b;A=full(a);B=full(b);A(:,[3 1])=B;
if or(full(a1)<>A) then bugmes();quit;end
b=sparse([],[],[6,3]);
a1=a;a1(:,[1 3 1])=b;A=full(a);B=full(b);A(:,[1 3 1])=B;
if or(full(a1)<>A) then bugmes();quit;end
// row column insertion
a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
b=sparse([1 2;1 3;3 3],-(1:3),[3,3]);
a1=a;a1(1,1)=sparse(30);A=full(a);A(1,1)=30;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1,6)=sparse(30);A=full(a);A(1,6)=30;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1,8)=sparse(30);A=full(a);A(1,8)=30;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1:3,1:3)=b;A=full(a);A(1:3,1:3)=full(b);
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1:3,6:8)=b;A=full(a);A(1:3,6:8)=full(b);
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(6:8,1:3)=b;A=full(a);A(6:8,1:3)=full(b);
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([3 2 1],1:3)=b;A=full(a);A([3 2 1],1:3)=full(b);
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 2 1],1:3)=b;A=full(a);A([1 2 1],1:3)=full(b);
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([3 2 1],[3 2 1])=b;A=full(a);A([3 2 1],[3 2 1])=full(b);
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 2 1],[3 2 1])=b;A=full(a);A([1 2 1],[3 2 1])=full(b);
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 2 1],[1 2 1])=b;A=full(a);A([1 2 1],[1 2 1])=full(b);
if or(full(a1)<>A) then bugmes();quit;end
//sparse full
a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
vt=11:16;
// full line insertion
//----------------------
a1=a;a1(1,:)=vt;A=full(a);A(1,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);A(3,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1,:)=vt;A=full(a);A(1,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);A(3,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(7,:)=vt;A=full(a);A(7,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(8,:)=vt;A=full(a);A(8,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 3],:)=[vt;2*vt];A=full(a);A([1 3],:)=[vt;2*vt];
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([3 1],:)=[vt;2*vt];A=full(a);A([3 1],:)=[vt;2*vt];
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 3 1],:)=[vt;2*vt;3*vt];A=full(a);A([1 3 1],:)=[vt;2*vt;3*vt];
if or(full(a1)<>A) then bugmes();quit;end
a=sparse([],[],[6,6]);
a1=a;a1(1,:)=vt;A=full(a);A(1,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);A(3,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1,:)=vt;A=full(a);A(1,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(3,:)=vt;A=full(a);A(3,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(7,:)=vt;A=full(a);A(7,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(8,:)=vt;A=full(a);A(8,:)=vt;
if or(full(a1)<>A) then bugmes();quit;end
b=[1:6;11:16];
a1=a;a1([1 3],:)=b;A=full(a);A([1 3],:)=b;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([3 1],:)=b;A=full(a);A([3 1],:)=b;
if or(full(a1)<>A) then bugmes();quit;end
b=[1:6;11:16;21:26];
a1=a;a1([1 3 1],:)=b;A=full(a);A([1 3 1],:)=b;
if or(full(a1)<>A) then bugmes();quit;end
// full column insertion
//----------------------
a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
v=(1:6)';
a1=a;a1(:,1)=v;A=full(a);A(:,1)=v;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);A(:,2)=v;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,1)=v;A=full(a);A(:,1)=v;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,2)=v;A=full(a);A(:,2)=v;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,3)=v;A=full(a);A(:,3)=v;
if or(full(a1)<>A) then bugmes();quit;end
//
a1=a;a1(:,7)=v;A=full(a);A(:,7)=v;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,8)=v;A=full(a);A(:,8)=v;
if or(full(a1)<>A) then bugmes();quit;end
b=[(1:6)' (11:16)'];
a1=a;a1(:,[1 3])=b;A=full(a);A(:,[1 3])=b;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(:,[3 1])=b;A=full(a);A(:,[3 1])=b;
if or(full(a1)<>A) then bugmes();quit;end
b=[(1:6)' (11:16)' (21:26)'];
a1=a;a1(:,[1 3 1])=b;A=full(a);A(:,[1 3 1])=b;
if or(full(a1)<>A) then bugmes();quit;end
//********
// row column insertion
a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
b=-[1 2 3;4 5 6;7 8 9];
a1=a;a1(1,1)=30;A=full(a);A(1,1)=30;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1,6)=30;A=full(a);A(1,6)=30;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1,8)=30;A=full(a);A(1,8)=30;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1:3,1:3)=b;A=full(a);A(1:3,1:3)=b;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(1:3,6:8)=b;A=full(a);A(1:3,6:8)=b;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1(6:8,1:3)=b;A=full(a);A(6:8,1:3)=b;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([3 2 1],1:3)=b;A=full(a);A([3 2 1],1:3)=b;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 2 1],1:3)=b;A=full(a);A([1 2 1],1:3)=b;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([3 2 1],[3 2 1])=b;A=full(a);A([3 2 1],[3 2 1])=b;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 2 1],[3 2 1])=b;A=full(a);A([1 2 1],[3 2 1])=b;
if or(full(a1)<>A) then bugmes();quit;end
a1=a;a1([1 2 1],[1 2 1])=b;A=full(a);A([1 2 1],[1 2 1])=b;
if or(full(a1)<>A) then bugmes();quit;end
// vector insertion
v=sparse([1 1;3 1;6 1],[10 11 12],[6 1]);
v1=v;v1(1)=33;V=full(v);V(1)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(2)=33;V=full(v);V(2)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(8)=33;V=full(v);V(8)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1([1 2 8])=[5;10;33];V=full(v);V([1 2 8])=[5;10;33];
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(:)=[];
if or(full(v1)<>[]) then bugmes();quit;end
v1=v;v1(1)=sparse(33);V=full(v);V(1)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(2)=sparse(33);V=full(v);V(2)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(8)=sparse(33);V=full(v);V(8)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1([1 2 8])=sparse([5;10;33]);V=full(v);V([1 2 8])=[5;10;33];
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1([1 2 1])=sparse([5;10;33]);V=full(v);V([1 2 1])=[5;10;33];
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(:)=[];
if or(full(v1)<>[]) then bugmes();quit;end
v1=v;v1(:)=sparse([2 1],44,[6 1]);V=full(v);V(:)=[0;44;0;0;0;0];
if or(full(v1)<>V) then bugmes();quit;end
v=v';
v1=v;v1(1)=33;V=full(v);V(1)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(2)=33;V=full(v);V(2)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(8)=33;V=full(v);V(8)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1([1 2 8])=[5 10 33];V=full(v);V([1 2 8])=[5 10 33];
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(1)=sparse(33);V=full(v);V(1)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(2)=sparse(33);V=full(v);V(2)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(8)=sparse(33);V=full(v);V(8)=33;
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1([1 2 8])=sparse([5 10 33]);V=full(v);V([1 2 8])=[5 10 33];
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1([1 2 1])=sparse([5 10 33]);V=full(v);V([1 2 1])=[5 10 33];
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(:)=sparse([1 2],44,[1,6]);V=full(v);V(:)=[0 44 0 0 0 0];
if or(full(v1)<>V) then bugmes();quit;end
v1=v;v1(1)=[];V=full(v);V(1)=[];
if or(full(v1)<>V) then bugmes();quit;end
//test des comparaisons
a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
b=sparse([1 6;1 2;6 5],[10;-1;-1],[6 6]);
if full(a==a)<>full(a)==full(a) then bugmes();quit;end
if full(a<>a)<>(full(a)<>full(a)) then bugmes();quit;end
if full(a>sparse(5))<>(full(a)>5) then bugmes();quit;end
if full(sparse(5)>a)<>(5>full(a)) then bugmes();quit;end
if full(b>a)<>(full(b)>full(a))  then bugmes();quit;end
if full(a==full(a))<>full(a)==full(a) then bugmes();quit;end
if full(a<>full(a))<>(full(a)<>full(a)) then bugmes();quit;end
if full(a>5)<>(full(a)>5) then bugmes();quit;end
if full(5>a)<>(5>full(a)) then bugmes();quit;end
if full(b>full(a))<>(full(b)>full(a))  then bugmes();quit;end
if full(full(a)==a)<>full(a)==full(a) then bugmes();quit;end
if full(full(a)<>a)<>(full(a)<>full(a)) then bugmes();quit;end
if full(full(a)>sparse(5))<>(full(a)>5) then bugmes();quit;end
if full(full(b)>a)<>(full(b)>full(a))  then bugmes();quit;end
a=sparse([1 1;3 1;6 1],[10 11 12],[6 1]);
if full(a==a)<>full(a)==full(a) then bugmes();quit;end
if full(a<>a)<>(full(a)<>full(a)) then bugmes();quit;end
if full(a>sparse(5))<>(full(a)>5) then bugmes();quit;end
if full(sparse(5)>a)<>(5>full(a)) then bugmes();quit;end
if full(a==full(a))<>full(a)==full(a) then bugmes();quit;end
if full(a<>full(a))<>(full(a)<>full(a)) then bugmes();quit;end
if full(a>5)<>(full(a)>5) then bugmes();quit;end
if full(5>a)<>(5>full(a)) then bugmes();quit;end
if full(full(a)==a)<>full(a)==full(a) then bugmes();quit;end
if full(full(a)<>a)<>(full(a)<>full(a)) then bugmes();quit;end
if full(full(a)>sparse(5))<>(full(a)>5) then bugmes();quit;end
a=a';
if full(a==a)<>full(a)==full(a) then bugmes();quit;end
if full(a<>a)<>(full(a)<>full(a)) then bugmes();quit;end
if full(a>sparse(5))<>(full(a)>5) then bugmes();quit;end
if full(sparse(5)>a)<>(5>full(a)) then bugmes();quit;end
if full(a==full(a))<>full(a)==full(a) then bugmes();quit;end
if full(a<>full(a))<>(full(a)<>full(a)) then bugmes();quit;end
if full(a>5)<>(full(a)>5) then bugmes();quit;end
if full(5>a)<>(5>full(a)) then bugmes();quit;end
if full(full(a)==a)<>full(a)==full(a) then bugmes();quit;end
if full(full(a)<>a)<>(full(a)<>full(a)) then bugmes();quit;end
if full(full(a)>sparse(5))<>(full(a)>5) then bugmes();quit;end
a=sparse([1 6;1 5;1 3;2 4;2 1;4 4;4 3;5 1;6 6],1:9,[6 6]);
b=sparse([1 6;1 2;5 5],[10;-1;-1],[5 6]);
if a==b<>%f then bugmes();quit;end
if a<>b<>%t then bugmes();quit;end
a=sparse([%f %t %f %t]);
b=sparse([%t %f %t %f]);
assert_checkequal(~a, b);