#include "stdafx.h"

// for bspline surface
#include <Geom_BSplineSurface.hxx>
#include <TColgp_Array2OfPnt.hxx>
#include <TColStd_Array1OfReal.hxx>;
#include <TColStd_Array1OfInteger.hxx>;

//end of bspline surface



void CImportExportDoc::OnArchHw()
{
// TODO: Add your command handler code here
//return 0;
TColgp_Array2OfPnt poles(0,7,0,5);//ucol,vrow
TColStd_Array1OfReal UKnots(0,12);
TColStd_Array1OfReal VKnots(0,10);

TColStd_Array1OfReal UKnots1(0,5);
TColStd_Array1OfReal VKnots1(0,3);
UKnots1.SetValue(0,0);UKnots1.SetValue(1,1);UKnots1.SetValue(2,2);UKnots1.SetValue(3,3);UKnots1.SetValue(4,4);UKnots1.SetValue(5,5);
VKnots1.SetValue(0,0);VKnots1.SetValue(1,1);VKnots1.SetValue(2,2);VKnots1.SetValue(3,3);
TColStd_Array1OfInteger UMults(0,5);
UMults.SetValue(0,4);UMults.SetValue(1,1);UMults.SetValue(2,1);UMults.SetValue(3,1);UMults.SetValue(4,1);
UMults.SetValue(5,4);
TColStd_Array1OfInteger VMults(0,3);
VMults.SetValue(0,4);VMults.SetValue(1,1);VMults.SetValue(2,1);
VMults.SetValue(3,4);
int n=7;//(ctr pnts-1)
int k=4;//degree+1
int n1=5;int k1=4;int t=0;int a=0;int b=0;

for (int j = 0;j<=n+k;j++)
{
    if (j<k)
{
        t = 0;a++; UKnots.SetValue(a,t);
}
  
    if (k<=j && j<=n)
{
        t = j-k+1;a++; UKnots.SetValue(a,t);
}
  
    if (j>n)
{
        t = n-k+2;a++; UKnots.SetValue(a,t);
}
}
 UKnots.SetValue(a,t);

for (int s = 0;s<=n1+k1;s++)
{
    if (s<k1)
{
        t = 0; b++; VKnots.SetValue(b,t);
}
  
    if (k1<=s && s<=n1)
{
        t = s-k1+1;b++;VKnots.SetValue(b,t);
}
  
    if (s>n1)
{
        t = n1-k1+2;b++; VKnots.SetValue(b,t);
}
}
 
cout<<"\n"<<a;
cout<<"\n"<<b;
for(int h1=1;h1<=UKnots.Length()-1;h1++)
cout<<"\n UKnots value"<<"\t"<<UKnots.Value(h1);

for(int g1=1;g1<=VKnots.Length()-1;g1++)
cout<<"\n VKnots value"<<"\t"<<VKnots.Value(g1);
 /* 
knots = [knots ; t];
end

knots_matrix = knots;
*/
//gp_Pnt P;
//poles.SetValue(1,1,P);
//poles.SetValue(1,1,P);
poles.SetValue(0,0,gp_Pnt(0,0,0));
poles.SetValue(1,0,gp_Pnt(0.4,0,0));
poles.SetValue(2,0,gp_Pnt(0.8,0,0));
poles.SetValue(3,0,gp_Pnt(1.2,0,0));
poles.SetValue(4,0,gp_Pnt(1.6,0,0));
poles.SetValue(5,0,gp_Pnt(2.0,0,0));
poles.SetValue(6,0,gp_Pnt(2.5,0,0));
poles.SetValue(7,0,gp_Pnt(3,0,0));

poles.SetValue(0,1,gp_Pnt(0,0.4,0));
poles.SetValue(1,1,gp_Pnt(0.4,0.4,0.2));
poles.SetValue(2,1,gp_Pnt(0.8,0.4,0.3));
poles.SetValue(3,1,gp_Pnt(1.2,0.4,0.1));
poles.SetValue(4,1,gp_Pnt(1.6,0.4,0));
poles.SetValue(5,1,gp_Pnt(2.0,0.4,0.2));
poles.SetValue(6,1,gp_Pnt(2.5,0.4,0.1));
poles.SetValue(7,1,gp_Pnt(3,0.4,0));

poles.SetValue(0,2,gp_Pnt(0,0.8,0));
poles.SetValue(1,2,gp_Pnt(0.4,0.8,0.5));
poles.SetValue(2,2,gp_Pnt(0.8,0.8,0.6));
poles.SetValue(3,2,gp_Pnt(1.2,0.8,0.2));
poles.SetValue(4,2,gp_Pnt(1.6,0.8,0.0));
poles.SetValue(5,2,gp_Pnt(2.0,0.8,0));
poles.SetValue(6,2,gp_Pnt(2.5,0.8,0.3));
poles.SetValue(7,2,gp_Pnt(3,0.8,0));

poles.SetValue(0,3,gp_Pnt(0,1.2,0));
poles.SetValue(1,3,gp_Pnt(0.4,1.2,0.4));
poles.SetValue(2,3,gp_Pnt(0.8,1.2,0.4));
poles.SetValue(3,3,gp_Pnt(1.2,1.2,0.2));
poles.SetValue(4,3,gp_Pnt(1.6,1.2,0.0));
poles.SetValue(5,3,gp_Pnt(2.0,1.2,0.0));
poles.SetValue(6,3,gp_Pnt(2.5,1.2,0.5));
poles.SetValue(7,3,gp_Pnt(3,1.2,0));

poles.SetValue(0,4,gp_Pnt(0,1.6,0));
poles.SetValue(1,4,gp_Pnt(0.4,1.6,0.3));
poles.SetValue(2,4,gp_Pnt(0.8,1.6,0.2));
poles.SetValue(3,4,gp_Pnt(1.2,1.6,0.1));
poles.SetValue(4,4,gp_Pnt(1.6,1.6,0.3));
poles.SetValue(5,4,gp_Pnt(2.0,1.6,0.2));
poles.SetValue(6,4,gp_Pnt(2.5,1.6,0.1));
poles.SetValue(7,4,gp_Pnt(3,1.6,0));

poles.SetValue(0,5,gp_Pnt(0,2.0,0));
poles.SetValue(1,5,gp_Pnt(0.4,2.0,0));
poles.SetValue(2,5,gp_Pnt(0.8,2.0,0));
poles.SetValue(3,5,gp_Pnt(1.2,2.0,0));
poles.SetValue(4,5,gp_Pnt(1.6,2.0,0));
poles.SetValue(5,5,gp_Pnt(2.0,2.0,0));
poles.SetValue(6,5,gp_Pnt(2.5,2.0,0));
poles.SetValue(7,5,gp_Pnt(3,2.0,0));

Handle (Geom_BSplineSurface) Surf= new Geom_BSplineSurface(poles,UKnots1,VKnots1,UMults,VMults,3,3,Standard_False,Standard_False);// = new Geom_BSplineSurface();

TopoDS_Face aFace = BRepBuilderAPI_MakeFace(Surf);
Handle (AIS_Shape) aShape = new AIS_Shape(aFace);
myAISContext->SetColor(aShape,Quantity_NOC_YELLOW,Standard_False);//projected wire
myAISContext->Display(aShape,Standard_False);

// TopoDS_Shape sh=Surf;
// Handle(Geom_Surf) plz= new Geom_Surf(Surf);
// Handle (AIS_Shape) ais= new AIS_Shape(sh);
// Handle(Geom_Surface) Sur = new Geom_Surface(Surf);
// Surf.Bounds(0,2,0,1);
// myAISContext->DisplayAll(1,1);//>Display(Surf);
//Standard_Integer d=Surf->NbVPoles();Standard_Integer e=Surf.NbUPoles();
//cout<<"\n\n No of poles"<<"\t"<<d<<"\t"<<e;
}