diff --git a/src/BOPTest/BOPTest.cxx b/src/BOPTest/BOPTest.cxx
index c69edf1c57..429ddde3d0 100644
--- a/src/BOPTest/BOPTest.cxx
+++ b/src/BOPTest/BOPTest.cxx
@@ -42,6 +42,7 @@ void  BOPTest::AllCommands(Draw_Interpretor& theCommands)
   BOPTest::APICommands       (theCommands);
   BOPTest::OptionCommands    (theCommands);
   BOPTest::HistoryCommands   (theCommands);
+  BOPTest::DebugCommands     (theCommands);
   BOPTest::CellsCommands     (theCommands);
 }
 //=======================================================================
diff --git a/src/BOPTest/BOPTest.hxx b/src/BOPTest/BOPTest.hxx
index 8b966ccccf..1969e079ba 100644
--- a/src/BOPTest/BOPTest.hxx
+++ b/src/BOPTest/BOPTest.hxx
@@ -55,9 +55,9 @@ public:
 
   Standard_EXPORT static void HistoryCommands  (Draw_Interpretor& aDI);
 
+  Standard_EXPORT static void DebugCommands  (Draw_Interpretor& aDI);
   Standard_EXPORT static void CellsCommands  (Draw_Interpretor& aDI);
 
-
 protected:
 
 
diff --git a/src/BOPTest/BOPTest_BOPCommands.cxx b/src/BOPTest/BOPTest_BOPCommands.cxx
index 8807317e08..068c2cb90b 100644
--- a/src/BOPTest/BOPTest_BOPCommands.cxx
+++ b/src/BOPTest/BOPTest_BOPCommands.cxx
@@ -22,7 +22,6 @@
 #include <BOPCol_ListOfShape.hxx>
 #include <BOPDS_DS.hxx>
 #include <BOPTest.hxx>
-#include <BOPTest_DrawableShape.hxx>
 #include <BOPTest_Objects.hxx>
 #include <BRep_Builder.hxx>
 #include <BRepAlgoAPI_BooleanOperation.hxx>
@@ -32,7 +31,6 @@
 #include <BRepAlgoAPI_Section.hxx>
 #include <DBRep.hxx>
 #include <Draw.hxx>
-#include <Draw_Color.hxx>
 #include <DrawTrSurf.hxx>
 #include <Geom2d_Curve.hxx>
 #include <Geom_Curve.hxx>
@@ -77,7 +75,6 @@ static Standard_Integer bfuse     (Draw_Interpretor&, Standard_Integer, const ch
 static Standard_Integer bcommon   (Draw_Interpretor&, Standard_Integer, const char**);
 //
 static Standard_Integer bopcurves (Draw_Interpretor&, Standard_Integer, const char**);
-static Standard_Integer bopnews   (Draw_Interpretor&, Standard_Integer, const char**);
 static Standard_Integer mkvolume   (Draw_Interpretor&, Standard_Integer, const char**);
 
 //=======================================================================
@@ -109,7 +106,6 @@ static Standard_Integer mkvolume   (Draw_Interpretor&, Standard_Integer, const c
   //
   theCommands.Add("bopcurves", "use bopcurves F1 F2 [-2d/-2d1/-2d2]",
                                                       __FILE__, bopcurves, g);
-  theCommands.Add("bopnews"  , "use  bopnews -v[e,f]"      , __FILE__, bopnews, g);
   theCommands.Add("mkvolume", "make solids from set of shapes.\nmkvolume r b1 b2 ... [-c] [-ni]", 
                   __FILE__, mkvolume , g);
 }
@@ -528,77 +524,6 @@ Standard_Integer bsmt (Draw_Interpretor& di,
   return 0;
 }
 //=======================================================================
-//function : bopnews
-//purpose  : 
-//=======================================================================
-Standard_Integer bopnews (Draw_Interpretor& di, 
-                          Standard_Integer n, 
-                          const char** a)
-{
-  if (n!=2) {
-    di << " use bopnews -v[e,f]\n";
-    return 0;
-  }
-  //
-  if (pPF==NULL) {
-    di << " Prepare BOPAlgo_PaveFiller first >bop S1 S2\n";
-    return 0;
-  }
-  //
-  char buf[32];
-  Standard_CString aText;
-  Standard_Integer i, i1, i2, iFound;
-  Draw_Color aTextColor(Draw_cyan);
-  TopAbs_ShapeEnum aT;
-  Handle(BOPTest_DrawableShape) aDShape;
-  //
-  const BOPDS_PDS& pDS=pPF->PDS();
-  //
-  aT=TopAbs_SHAPE;
-  if (!strcmp (a[1], "-f")) {
-    aT=TopAbs_FACE;
-  }
-  else if (!strcmp (a[1], "-e")){
-    aT=TopAbs_EDGE;
-  }
-  else if (!strcmp (a[1], "-v")){
-    aT=TopAbs_VERTEX;
-  }
-  else {
-    di << " use bopnews -v[e,f]\n";
-    return 0;
-  }
-  //
-  iFound=0;
-  i1=pDS->NbSourceShapes();
-  i2=pDS->NbShapes();
-  for (i=i1; i<i2; ++i) {
-    const BOPDS_ShapeInfo& aSI=pDS->ShapeInfo(i);
-    if (aSI.ShapeType()==aT) {
-      const TopoDS_Shape& aS=aSI.Shape();
-      //
-      Sprintf (buf, "z%d", i);
-      aText=buf;
-      aDShape=new BOPTest_DrawableShape (aS, aText, aTextColor);
-      Draw::Set (aText, aDShape);
-      //
-      Sprintf (buf, " z%d", i);
-      di << buf;
-      //
-      iFound=1;
-    }
-  }
-  //
-  if (iFound) {
-    di << "\n";
-  }
-  else {
-    di << " not found\n";
-  }
-  //
-  return 0;
-}
-//=======================================================================
 //function : bopcurves
 //purpose  : 
 //=======================================================================
diff --git a/src/BOPTest/BOPTest_DebugCommands.cxx b/src/BOPTest/BOPTest_DebugCommands.cxx
new file mode 100644
index 0000000000..400c589aed
--- /dev/null
+++ b/src/BOPTest/BOPTest_DebugCommands.cxx
@@ -0,0 +1,1597 @@
+// Created by: Eugeny MALTCHIKOV
+// Copyright (c) 2015 OPEN CASCADE SAS
+//
+// This file is part of Open CASCADE Technology software library.
+//
+// This library is free software; you can redistribute it and/or modify it under
+// the terms of the GNU Lesser General Public License version 2.1 as published
+// by the Free Software Foundation, with special exception defined in the file
+// OCCT_LGPL_EXCEPTION.txt. Consult the file LICENSE_LGPL_21.txt included in OCCT
+// distribution for complete text of the license and disclaimer of any warranty.
+//
+// Alternatively, this file may be used under the terms of Open CASCADE
+// commercial license or contractual agreement.
+
+
+#include <BOPTest.hxx>
+#include <BOPTest_Objects.hxx>
+#include <BOPTest_DrawableShape.hxx>
+
+#include <Draw.hxx>
+#include <Draw_Color.hxx>
+
+#include <DBRep.hxx>
+
+#include <BRep_Builder.hxx>
+
+#include <TopoDS_Compound.hxx>
+
+#include <TopExp_Explorer.hxx>
+
+#include <BOPDS_DS.hxx>
+#include <BOPDS_Iterator.hxx>
+#include <BOPDS_MapOfCommonBlock.hxx>
+
+#include <BOPAlgo_Builder.hxx>
+#include <BOPAlgo_BuilderFace.hxx>
+#include <BOPAlgo_BuilderSolid.hxx>
+
+#include <BOPCol_MapOfShape.hxx>
+#include <BOPCol_DataMapOfShapeShape.hxx>
+#include <BOPCol_DataMapOfShapeListOfShape.hxx>
+
+
+static
+  void GetTypeByName(const char* theName,
+                     TopAbs_ShapeEnum& theType);
+static
+  void GetNameByType(const TopAbs_ShapeEnum& theType,
+                     char* theName);
+
+template <class InterfType> static
+  void DumpInterfs(const BOPCol_NCVector<InterfType>& theVInterf,
+                   Draw_Interpretor& di);
+
+template <class InterfType> static
+  void SearchNewIndex(const char* theCType,
+                      const Standard_Integer theInd,
+                      const BOPCol_NCVector<InterfType>& theVInterf,
+                      Draw_Interpretor& di);
+static 
+  Standard_Integer bopfinfo(Draw_Interpretor& di,
+                            Standard_Integer n,
+                            const char** a,
+                            const Standard_Integer iPriz);
+
+
+// commands
+// 1. filler commands
+// 1.1 DS commands
+static Standard_Integer bopds       (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopiterator (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopinterf   (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopnews     (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopwho      (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopindex    (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopsd       (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopsc       (Draw_Interpretor&, Standard_Integer, const char**);
+
+// 1.2 pave blocks commands
+static Standard_Integer boppb       (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopcb       (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopsp       (Draw_Interpretor&, Standard_Integer, const char**);
+
+// 1.3 face info commands
+static Standard_Integer bopfon      (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopfin      (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopfsc      (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopfav      (Draw_Interpretor&, Standard_Integer, const char**);
+
+// 2. builder commands
+// 2.1 images commands
+static Standard_Integer bopimage    (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer boporigin   (Draw_Interpretor&, Standard_Integer, const char**);
+static Standard_Integer bopfsd      (Draw_Interpretor&, Standard_Integer, const char**);
+
+// 2.2 building faces
+static Standard_Integer bopbface    (Draw_Interpretor&, Standard_Integer, const char**);
+// 2.3 building solids
+static Standard_Integer bopbsolid   (Draw_Interpretor&, Standard_Integer, const char**);
+
+
+//=======================================================================
+//function : DebugCommands
+//purpose  : 
+//=======================================================================
+void BOPTest::DebugCommands(Draw_Interpretor& theCommands)
+{
+  static Standard_Boolean done = Standard_False;
+  if (done) return;
+  done = Standard_True;
+  // Chapter's name
+  const char* g = "BOPTest commands";
+  // Commands
+  theCommands.Add("bopds",
+                  "Shows the shapes from DS. Use: bopds [v/e/w/f/sh/s/cs/c]",
+                  __FILE__, bopds,       g);
+  theCommands.Add("bopiterator", 
+                  "Shows the pairs of interfered shapes. Use: bopiterator [type1 type2]",
+                  __FILE__, bopiterator, g);
+  theCommands.Add("bopinterf", "Shows interferences of given type. Use: bopinterf type1 type2",
+                  __FILE__, bopinterf,   g);
+  theCommands.Add("bopnews", "Shows the newly created shapes. Use: bopnews [v,e,f]",
+                  __FILE__, bopnews,     g);
+  theCommands.Add("bopwho", "Shows where the new shape was created. Use: bopwho #",
+                  __FILE__, bopwho,      g);
+  theCommands.Add("bopindex", "Gets the index of the shape in the DS. Use: bopindex s",
+                  __FILE__, bopindex,    g);
+  theCommands.Add("bopsd", "Gets the Same domain shape. Use: bopsd #",
+                   __FILE__, bopsd,      g);
+  theCommands.Add("bopsc", "Shows the section curves. Use: bopsc [nF1 nF2]",
+                  __FILE__, bopsc,       g);
+  theCommands.Add("boppb", "Shows information about pave blocks. Use: boppb [#e]",
+                  __FILE__, boppb,       g);
+  theCommands.Add("bopcb", "Shows information about common blocks. Use: bopcb [#e]",
+                  __FILE__, bopcb,       g);
+  theCommands.Add("bopsp", "Shows the splits of edges. Use: bopsp [#e]",
+                  __FILE__, bopsp,       g);
+  theCommands.Add("bopfon", "Shows ON information for the face. Use: bopfon #f",
+                  __FILE__, bopfon,      g);
+  theCommands.Add("bopfin", "Shows IN information for the face. Use: bopfin #f",
+                  __FILE__, bopfin,      g);
+  theCommands.Add("bopfsc", "Shows SC information for the face. Use: bopfsc #f",
+                  __FILE__, bopfsc,      g);
+  theCommands.Add("bopfav", "Shows information about alone vertices for the face. Use: bopfav #f",
+                  __FILE__, bopfav,      g);
+  theCommands.Add("bopimage", "Shows split parts of the shape. Use: bopimage s",
+                  __FILE__, bopimage,    g);
+  theCommands.Add("boporigin", "Shows the original shape for the shape. Use: boporigin s",
+                  __FILE__, boporigin,   g);
+  theCommands.Add("bopfsd", "Shows SD faces for the face: Use: bopfsd f",
+                  __FILE__, bopfsd,      g);
+  theCommands.Add("bopbsolid", "Build solids from set of shared faces. Use: bopbsolid sr cx",
+                  __FILE__, bopbsolid,   g);
+  theCommands.Add("bopbface", "Splits the face by set of shared edges. Use: bopbface fr cx",
+                  __FILE__, bopbface,    g);
+}
+
+//=======================================================================
+//function : bopds
+//purpose  : 
+//=======================================================================
+Standard_Integer bopds (Draw_Interpretor& di,
+                        Standard_Integer n,
+                        const char** a)
+{
+  if (n > 2) {
+    di << "Shows the shapes from DS. Use: bopds [v/e/w/f/sh/s/cs/c]\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  char buf[32];
+  Standard_CString aText;
+  Standard_Integer i, aNbS;
+  TopAbs_ShapeEnum aType, aTypeShape;
+  Draw_Color aTextColor(Draw_cyan);
+  //
+  BOPDS_DS& aDS = *pDS;
+  aNbS = aDS.NbSourceShapes();
+  //
+  aType = TopAbs_SHAPE;
+  if (n == 2) {
+    GetTypeByName(a[1], aType);
+  }
+  //
+  for (i = 0; i < aNbS; ++i) {
+    const BOPDS_ShapeInfo& aSI = aDS.ShapeInfo(i);
+    const TopoDS_Shape& aS = aSI.Shape();
+    aTypeShape = aSI.ShapeType();
+    if (n == 1) {
+      if (aTypeShape == TopAbs_COMPOUND ||
+          aTypeShape == TopAbs_COMPSOLID ||
+          aTypeShape == TopAbs_SOLID ||
+          aTypeShape == TopAbs_SHELL || 
+          aTypeShape == TopAbs_WIRE) {
+        continue;
+      }
+    }
+    else {
+      if (aTypeShape != aType) {
+        continue;
+      }
+    }
+    //
+    Sprintf(buf, "z%d", i);
+    aText = buf;
+    Handle(BOPTest_DrawableShape) aDShape = new BOPTest_DrawableShape(aS, aText, aTextColor);
+    Draw::Set(aText, aDShape);
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopiterator
+//purpose  : 
+//=======================================================================
+Standard_Integer bopiterator (Draw_Interpretor& di,
+                              Standard_Integer n,
+                              const char** a)
+{
+  if (n != 1 && n != 3) {
+    di << "Shows the pairs of interfered shapes. Use: bopiterator [type1 type2]\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  Standard_Boolean bFlag;
+  Standard_Integer n1, n2;
+  char buf[64], aST1[10], aST2[10];
+  BOPDS_Iterator aIt;
+  //
+  BOPDS_DS& aDS = *pDS;
+  aIt.SetDS(&aDS);
+  aIt.Prepare();
+  //
+  if (n == 1) {
+    // type has not been defined. show all pairs
+    Standard_Integer i, j;
+    TopAbs_ShapeEnum aT[4] = 
+      {TopAbs_VERTEX, TopAbs_EDGE, TopAbs_FACE, TopAbs_SOLID};
+    for (i = 0; i < 4; ++i) {
+      GetNameByType(aT[i], aST1);
+      //
+      for (j = 0; j < 4; ++j) {
+        GetNameByType(aT[j], aST2);
+        //
+        aIt.Initialize(aT[i], aT[j]);
+        for (; aIt.More(); aIt.Next()) {
+          aIt.Value(n1, n2, bFlag);
+          //
+          Sprintf(buf, "%s/%s: (z%d z%d)\n", aST1, aST2, n1, n2);
+          di << buf;
+        }
+      }
+    }
+  }
+  else if (n == 3) {
+    TopAbs_ShapeEnum aT1, aT2;
+    //
+    GetTypeByName(a[1], aT1);
+    GetTypeByName(a[2], aT2);
+    //
+    GetNameByType(aT1, aST1);
+    GetNameByType(aT2, aST2);
+    //
+    aIt.Initialize(aT1, aT2);
+    for (; aIt.More(); aIt.Next()) {
+      aIt.Value(n1, n2, bFlag);
+      //
+      Sprintf(buf, "%s/%s: (z%d z%d)\n", aST1, aST2, n1, n2);
+      di << buf;
+    }
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopinterf
+//purpose  : 
+//=======================================================================
+Standard_Integer bopinterf (Draw_Interpretor& di,
+                            Standard_Integer n,
+                            const char** a)
+{
+  if (n != 3) {
+    di << "Shows interferences of given type. Use: bopinterf type1 type2\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  TopAbs_ShapeEnum aT1, aT2;
+  GetTypeByName(a[1], aT1);
+  GetTypeByName(a[2], aT2);
+  //
+  if (aT1 == TopAbs_VERTEX && aT2 == TopAbs_VERTEX) {
+    di << "V/V: ";
+    DumpInterfs<BOPDS_InterfVV>(pDS->InterfVV(), di);
+  }
+  else if ((aT1 == TopAbs_VERTEX && aT2 == TopAbs_EDGE) ||
+           (aT2 == TopAbs_VERTEX && aT1 == TopAbs_EDGE)) {
+    di << "V/E: ";
+    DumpInterfs<BOPDS_InterfVE>(pDS->InterfVE(), di);
+  }
+  else if (aT1 == TopAbs_EDGE && aT2 == TopAbs_EDGE) {
+    di << "E/E: ";
+    DumpInterfs<BOPDS_InterfEE>(pDS->InterfEE(), di);
+  }
+  else if ((aT1 == TopAbs_VERTEX && aT2 == TopAbs_FACE) ||
+           (aT2 == TopAbs_VERTEX && aT1 == TopAbs_FACE)) {
+    di << "V/F: ";
+    DumpInterfs<BOPDS_InterfVF>(pDS->InterfVF(), di);
+  }
+  else if ((aT1 == TopAbs_EDGE && aT2 == TopAbs_FACE) ||
+           (aT2 == TopAbs_EDGE && aT1 == TopAbs_FACE)) {
+    di << "E/F: ";
+    DumpInterfs<BOPDS_InterfEF>(pDS->InterfEF(), di);
+  }
+  else if (aT1 == TopAbs_FACE && aT2 == TopAbs_FACE) {
+    di << "F/F: ";
+    DumpInterfs<BOPDS_InterfFF>(pDS->InterfFF(), di);
+  }
+  else if ((aT1 == TopAbs_VERTEX && aT2 == TopAbs_SOLID) ||
+           (aT2 == TopAbs_VERTEX && aT1 == TopAbs_SOLID)) {
+    di << "V/S: ";
+    DumpInterfs<BOPDS_InterfVZ>(pDS->InterfVZ(), di);
+  }
+  else if ((aT1 == TopAbs_EDGE && aT2 == TopAbs_SOLID) ||
+           (aT2 == TopAbs_EDGE && aT1 == TopAbs_SOLID)) {
+    di << "E/S: ";
+    DumpInterfs<BOPDS_InterfEZ>(pDS->InterfEZ(), di);
+  }
+  else if ((aT1 == TopAbs_FACE && aT2 == TopAbs_SOLID) ||
+           (aT2 == TopAbs_FACE && aT1 == TopAbs_SOLID)) {
+    di << "F/S: ";
+    DumpInterfs<BOPDS_InterfFZ>(pDS->InterfFZ(), di);
+  }
+  else if (aT1 == TopAbs_SOLID && aT2 == TopAbs_SOLID) {
+    di << "S/S: ";
+    DumpInterfs<BOPDS_InterfZZ>(pDS->InterfZZ(), di);
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopwho
+//purpose  : 
+//=======================================================================
+Standard_Integer bopwho (Draw_Interpretor& di,
+                         Standard_Integer n,
+                         const char** a)
+{
+  if (n != 2) {
+    di << "Shows where the new shape was created. Use: bopwho #\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS=BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 0;
+  }
+  //
+  Standard_Integer ind = Draw::Atoi(a[1]);
+  if (ind <= 0) {
+    di << " Index must be grater than 0\n";
+    return 1;
+  }
+  //
+  Standard_Integer i1, i2;
+  //
+  i1 = 0;
+  i2 = pDS->NbShapes();
+  if (ind < i1 || ind > i2) {
+    di << " DS does not contain the shape\n";
+    return 1;
+  }
+  //
+  if (!pDS->IsNewShape(ind)) {
+    Standard_Integer iRank = pDS->Rank(ind);
+    di << " Rank: " << iRank << "\n";
+    return 0;
+  }
+  //
+  // the shape is new
+  di << "the shape is new\n";
+  //
+  const BOPDS_ShapeInfo& aSI = pDS->ShapeInfo(ind);
+  if (aSI.ShapeType() != TopAbs_VERTEX) {
+    return 0;
+  }
+  // search among interfs
+  BOPDS_VectorOfInterfVV& aVVs = pDS->InterfVV();
+  BOPDS_VectorOfInterfVE& aVEs = pDS->InterfVE();
+  BOPDS_VectorOfInterfEE& aEEs = pDS->InterfEE();
+  BOPDS_VectorOfInterfVF& aVFs = pDS->InterfVF();
+  BOPDS_VectorOfInterfEF& aEFs = pDS->InterfEF();
+  //
+  SearchNewIndex<BOPDS_InterfVV> ("V/V: ", ind, aVVs, di);
+  SearchNewIndex<BOPDS_InterfVE> ("V/E: ", ind, aVEs, di);
+  SearchNewIndex<BOPDS_InterfEE> ("E/E: ", ind, aEEs, di);
+  SearchNewIndex<BOPDS_InterfVF> ("V/F: ", ind, aVFs, di);
+  SearchNewIndex<BOPDS_InterfEF> ("E/F: ", ind, aEFs, di);
+  //
+  //--------------------------------------FF
+  char buf[64];
+  Standard_Boolean bFound;
+  Standard_Integer i, n1, n2, k, aNb, aNbC, aNbP, nV1, nV2;
+  BOPDS_ListIteratorOfListOfPaveBlock aItLPB;
+  //
+  bFound = Standard_False;
+  BOPDS_VectorOfInterfFF& aFFs = pDS->InterfFF();
+  aNb = aFFs.Extent();
+  for (i = 0; i < aNb; ++i) {
+    const BOPDS_InterfFF& anInt = aFFs(i);
+    anInt.Indices(n1, n2);
+    //
+    const BOPDS_VectorOfCurve& aVNC = anInt.Curves();
+    aNbC = aVNC.Extent();
+    for (k = 0; k < aNbC; ++k) {
+      const BOPDS_Curve& aNC = aVNC(k);
+      const BOPDS_ListOfPaveBlock& aLPB = aNC.PaveBlocks(); 
+      aItLPB.Initialize(aLPB);
+      for(; aItLPB.More(); aItLPB.Next()) {
+        const Handle(BOPDS_PaveBlock)& aPB = aItLPB.Value();
+        aPB->Indices(nV1, nV2);
+        if (ind == nV1 || ind == nV2) {
+          if (!bFound) {
+            di << " FF curves: ";
+            bFound = Standard_True;
+          }
+          Sprintf (buf,"(%d, %d) ", n1, n2);
+          di << buf;
+          break;
+        }
+      }
+    }//for (k=0; k<aNbC; ++k) 
+    if (bFound) {
+      di << "\n";
+    }
+    //
+    bFound = Standard_False;
+    const BOPDS_VectorOfPoint& aVNP = anInt.Points();
+    aNbP = aVNP.Extent();
+    for (k = 0; k < aNbP; ++k) {
+      const BOPDS_Point& aNP = aVNP(k);
+      nV1 = aNP.Index();
+      if (ind == nV1) {
+        if (!bFound) {
+          di << " FF points: ";
+          bFound = Standard_True;
+        }
+        sprintf (buf,"(%d, %d) ", n1, n2);
+        di << buf;
+      }
+    }//for (k=0; k<aNbP; ++k) 
+    if (bFound) {
+      di << "\n";
+    }
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopnews
+//purpose  : 
+//=======================================================================
+Standard_Integer bopnews(Draw_Interpretor& di,
+                         Standard_Integer n,
+                         const char** a)
+{ 
+  if (n != 2) {
+    di << "Shows the newly created shapes. Use: bopnews v/e/f\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  TopAbs_ShapeEnum aType;
+  GetTypeByName(a[1], aType);
+  //
+  if (!(aType == TopAbs_VERTEX ||
+        aType == TopAbs_EDGE ||
+        aType == TopAbs_FACE)) {
+    di << "Use: bopnews v/e/f\n";
+    return 1;
+  }
+  //
+  char buf[32];
+  Standard_CString aText;
+  Standard_Boolean bFound;
+  Standard_Integer i, i1, i2;
+  Draw_Color aTextColor(Draw_cyan);
+  Handle(BOPTest_DrawableShape) aDShape;
+  //
+  bFound = Standard_False;
+  i1 = pDS->NbSourceShapes();
+  i2 = pDS->NbShapes();
+  for (i = i1; i < i2; ++i) {
+    const BOPDS_ShapeInfo& aSI = pDS->ShapeInfo(i);
+    if (aSI.ShapeType() == aType) {
+      const TopoDS_Shape& aS = aSI.Shape();
+      //
+      sprintf (buf, "z%d", i);
+      aText = buf;
+      aDShape = new BOPTest_DrawableShape(aS, aText, aTextColor);
+      Draw::Set(aText, aDShape);
+      //
+      sprintf (buf, "z%d ", i);
+      di << buf;
+      //
+      bFound = Standard_True;
+    }
+  }
+  //
+  if (bFound) {
+    di << "\n";
+  }
+  else {
+    di << " No new shapes found\n";
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopindex
+//purpose  : 
+//=======================================================================
+Standard_Integer bopindex (Draw_Interpretor& di,
+                           Standard_Integer n,
+                           const char** a)
+{
+  if (n != 2) {
+    di << "Gets the index of the shape in the DS. Use: bopindex s\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS=BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  TopoDS_Shape aS = DBRep::Get(a[1]);
+  if (aS.IsNull()) {
+    di << a[1] << " is a null shape\n";
+    return 1;
+  }
+  //
+  Standard_Integer ind = pDS->Index(aS);
+  Standard_Boolean bFound = (ind > 0);
+  if (bFound) {
+    di << " Index: " << ind << "\n";
+  }
+  else {
+    Standard_Integer i1, i2;
+    //
+    i1 = pDS->NbSourceShapes();
+    i2 = pDS->NbShapes();
+    for (ind = i1; ind < i2; ++ind) {
+      const TopoDS_Shape& aSx = pDS->Shape(ind);
+      if (aSx.IsSame(aS)) {
+        di << " Index: " << ind << "\n";
+        bFound = Standard_True;
+        break;
+      }
+    }
+  }
+  //
+  if (!bFound) {
+    di << " DS does not contain the shape\n";
+  }
+  //
+  return 0;
+}
+  
+//=======================================================================
+//function : bopsd
+//purpose  : 
+//=======================================================================
+Standard_Integer bopsd(Draw_Interpretor& di,
+                       Standard_Integer n,
+                       const char** a)
+{ 
+  if (n != 2) {
+    di << "Gets the Same domain shape. Use: bopsd #\n";
+    return 0;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 0;
+  }
+  //
+  char buf[32];
+  Standard_Boolean bHasSD;
+  Standard_Integer ind, i1, i2, iSD;
+  //
+  ind = Draw::Atoi(a[1]);
+  //
+  i1 = 0;
+  i2 = pDS->NbShapes();
+  if (ind < i1 || ind > i2) {
+    di << " DS does not contain the shape\n";
+    return 0;
+  }
+  //
+  bHasSD = pDS->HasShapeSD(ind, iSD);
+  if (bHasSD) {
+    Sprintf(buf, " Shape %d has SD shape %d\n", ind, iSD);
+    di << buf;
+  }
+  else {
+    Sprintf (buf, " Shape: %d has no SD shape\n", ind);
+    di << buf;
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopsc
+//purpose  : 
+//=======================================================================
+Standard_Integer bopsc(Draw_Interpretor& di,
+                       Standard_Integer n,
+                       const char** a)
+{
+  if (n != 1 && n != 3) {
+    di << "Shows the section curves. Use: bopsc [nF1 nF2]\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS=BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 0;
+  }
+  //
+  char buf[32];
+  Standard_CString aText;
+  Draw_Color aTextColor(Draw_cyan);
+  Standard_Integer nSF1, nSF2, nF1, nF2;
+  Standard_Integer aNb, j, iCnt, k, iPriz, aNbC, aNbP, nSp;
+  Standard_Integer iX;
+  Handle(BOPTest_DrawableShape) aDShape;
+  BOPDS_ListIteratorOfListOfPaveBlock aItLPB;
+  //
+  nSF1 = nSF2 = -1;
+  if (n == 3) {
+    nSF1 = Draw::Atoi(a[1]);
+    nSF2 = Draw::Atoi(a[2]);
+  }
+  //
+  BOPDS_VectorOfInterfFF& aFFs = pDS->InterfFF();
+  //
+  iCnt = 0;
+  iPriz = 0;
+  aNb = aFFs.Extent();
+  for (j = 0; j < aNb; ++j) {
+    const BOPDS_InterfFF& aFF = aFFs(j);
+    if (n == 3) {
+      if (!aFF.Contains(nSF1) || !aFF.Contains(nSF2)) {
+        continue;
+      }
+      iPriz = 1;
+    }
+    //
+    aFF.Indices(nF1, nF2);
+    //
+    iX = 0;
+    const BOPDS_VectorOfCurve& aVNC = aFF.Curves();
+    aNbC = aVNC.Extent();
+    for (k = 0; k < aNbC; ++k) {
+      const BOPDS_Curve& aNC = aVNC(k);
+      const BOPDS_ListOfPaveBlock& aLPB = aNC.PaveBlocks();
+      aItLPB.Initialize(aLPB);
+      for(; aItLPB.More(); aItLPB.Next()) {
+        const Handle(BOPDS_PaveBlock)& aPB = aItLPB.Value();
+        if (!aPB->HasEdge(nSp)) {
+          continue;
+        }
+        //
+        if (!iX) {
+          Sprintf (buf, "[%d %d] section edges: ", nF1, nF2);
+          di << buf;
+          iX = 1;
+        }
+        sprintf (buf, "t_%d_%d", k, nSp);
+        di << buf;
+        //
+        const TopoDS_Shape& aSp = pDS->Shape(nSp);
+        aText = buf;
+        aDShape = new BOPTest_DrawableShape(aSp, aText, aTextColor);
+        Draw::Set(aText, aDShape);
+        di << " ";
+        ++iCnt;
+      }
+    }
+    if (iX) {
+      di << "\n";
+    }
+    //
+    iX = 0;
+    const BOPDS_VectorOfPoint& aVNP = aFF.Points();
+    aNbP = aVNP.Extent();
+    for (k = 0; k < aNbP; ++k) {
+      const BOPDS_Point& aNP = aVNP(k);
+      nSp = aNP.Index();
+      if (nSp < 0) {
+        continue;
+      }
+      if (!iX) {
+        sprintf (buf, "[%d %d] section vertices: ", nF1, nF2);
+        di << buf;
+        iX = 1;
+      }
+      sprintf (buf, "p_%d_%d", k, nSp);
+      di << buf;
+      //
+      const TopoDS_Shape& aSp = pDS->Shape(nSp);
+      aText = buf;
+      aDShape = new BOPTest_DrawableShape(aSp, aText, aTextColor);
+      Draw::Set(aText, aDShape);
+      di << " ";
+      ++iCnt;
+    }
+    if (iX) {
+      di << "\n";
+    }
+    //
+    if (iPriz) {
+      break;
+    }
+  }// for (j=0; j<aNb; ++j) {
+  //
+  if (iCnt) {
+    di << "\n";
+  }
+  else {
+    di << " no sections found\n";
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : boppb
+//purpose  : 
+//=======================================================================
+Standard_Integer boppb(Draw_Interpretor& di,
+                       Standard_Integer n,
+                       const char** a)
+{ 
+  if (n > 2) {
+    di << "Shows information about pave blocks. Use: boppb [#e]\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  Standard_Boolean bHasPaveBlocks;
+  Standard_Integer ind, i1, i2;
+  TopAbs_ShapeEnum aType;
+  BOPDS_ListIteratorOfListOfPaveBlock aItPB;
+  //
+  i1 = 0;
+  i2 = pDS->NbSourceShapes();
+  if (n == 2) {
+    ind = Draw::Atoi(a[1]);
+    i1 = ind;
+    i2 = ind + 1;
+  }
+  //
+  for (ind = i1; ind < i2; ++ind) {
+    const BOPDS_ShapeInfo& aSI = pDS->ShapeInfo(ind);
+    aType = aSI.ShapeType();
+    if (aType != TopAbs_EDGE) {
+      continue;
+    }
+    //
+    bHasPaveBlocks = pDS->HasPaveBlocks(ind);
+    if (!bHasPaveBlocks) {
+      continue;
+    }
+    //
+    const BOPDS_ListOfPaveBlock& aLPB = pDS->PaveBlocks(ind);
+    aItPB.Initialize(aLPB);
+    for (; aItPB.More(); aItPB.Next()) {
+      const Handle(BOPDS_PaveBlock)& aPB = aItPB.Value();
+      aPB->Dump();
+      printf("\n");
+    }
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopcb
+//purpose  : 
+//=======================================================================
+Standard_Integer bopcb(Draw_Interpretor& di,
+                       Standard_Integer n,
+                       const char** a)
+{ 
+  if (n > 2) {
+    di << "Shows information about common blocks. Use: bopcb [#e]\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  Standard_Boolean bHasPaveBlocks;
+  Standard_Integer ind, i1, i2;
+  TopAbs_ShapeEnum aType;
+  BOPDS_ListIteratorOfListOfPaveBlock aItPB;
+  BOPDS_MapOfCommonBlock aMCB;
+  //
+  i1 = 0;
+  i2 = pDS->NbSourceShapes();
+  if (n == 2) {
+    ind = Draw::Atoi(a[1]);
+    i1 = ind;
+    i2 = ind + 1;
+  }
+  //
+  for (ind = i1; ind < i2; ++ind) {
+    const BOPDS_ShapeInfo& aSI = pDS->ShapeInfo(ind);
+    aType = aSI.ShapeType();
+    if (aType != TopAbs_EDGE) {
+      continue;
+    }
+    //
+    bHasPaveBlocks = pDS->HasPaveBlocks(ind);
+    if (!bHasPaveBlocks) {
+      continue;
+    }
+    //
+    const BOPDS_ListOfPaveBlock& aLPB = pDS->PaveBlocks(ind);
+    aItPB.Initialize(aLPB);
+    for (; aItPB.More(); aItPB.Next()) {
+      const Handle(BOPDS_PaveBlock)& aPB = aItPB.Value();
+      if (pDS->IsCommonBlock(aPB)) {
+        const Handle(BOPDS_CommonBlock)& aCB = pDS->CommonBlock(aPB);
+        if(aMCB.Add(aCB)) {
+          aCB->Dump();
+          printf("\n");
+        }
+      }
+    }
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopsp
+//purpose  : 
+//=======================================================================
+Standard_Integer bopsp(Draw_Interpretor& di,
+                       Standard_Integer n,
+                       const char** a)
+{
+  if (n > 2) {
+    di << "Shows the splits of edges. Use: bopsp [#e]\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  char buf[32];
+  Standard_Boolean bHasPaveBlocks;
+  Standard_Integer ind, i1, i2, nSp;
+  TopAbs_ShapeEnum aType;
+  BOPDS_ListIteratorOfListOfPaveBlock aItPB;
+  Standard_CString aText;
+  Draw_Color aTextColor(Draw_cyan);
+  Handle(BOPTest_DrawableShape) aDShape;
+  //
+  i1 = 0;
+  i2 = pDS->NbSourceShapes();
+  if (n == 2) {
+    ind = Draw::Atoi(a[1]);
+    i1 = ind;
+    i2 = ind + 1;
+  }
+  //
+  for (ind = i1; ind < i2; ++ind) {
+    const BOPDS_ShapeInfo& aSI = pDS->ShapeInfo(ind);
+    aType = aSI.ShapeType();
+    if (aType != TopAbs_EDGE) {
+      continue;
+    }
+    //
+    bHasPaveBlocks = pDS->HasPaveBlocks(ind);
+    if (!bHasPaveBlocks) {
+      continue;
+    }
+    //
+    di << "Edge " << ind << ": ";
+    //
+    const BOPDS_ListOfPaveBlock& aLPB = pDS->PaveBlocks(ind);
+    aItPB.Initialize(aLPB);
+    for (; aItPB.More(); aItPB.Next()) {
+      const Handle(BOPDS_PaveBlock)& aPB = aItPB.Value();
+      nSp = aPB->Edge();
+      const TopoDS_Shape& aSp = pDS->Shape(nSp);
+      //
+      Sprintf(buf, "z%d_%d", ind, nSp);
+      aText = buf;
+      aDShape = new BOPTest_DrawableShape(aSp, aText, aTextColor);
+      Draw::Set(aText, aDShape);
+      di << buf << " ";
+    }
+    di << "\n";
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopfon
+//purpose  : 
+//=======================================================================
+Standard_Integer bopfon(Draw_Interpretor& di,
+                        Standard_Integer n,
+                        const char** a)
+{
+  return bopfinfo(di, n, a, 0);
+}
+
+//=======================================================================
+//function : bopfin
+//purpose  : 
+//=======================================================================
+Standard_Integer bopfin(Draw_Interpretor& di,
+                        Standard_Integer n,
+                        const char** a)
+{
+  return bopfinfo(di, n, a, 1);
+}
+
+//=======================================================================
+//function : bopfspsc
+//purpose  : 
+//=======================================================================
+Standard_Integer bopfsc(Draw_Interpretor& di,
+                        Standard_Integer n,
+                        const char** a)
+{
+  return bopfinfo(di, n, a, 2);
+}
+
+//=======================================================================
+//function : bopfinfo
+//purpose  : 
+//=======================================================================
+Standard_Integer bopfinfo(Draw_Interpretor& di,
+                          Standard_Integer n,
+                          const char** a,
+                          const Standard_Integer iPriz)
+{
+  if (n != 2) {
+    di << "Shows " << ((iPriz == 0) ? "ON" : ((iPriz == 1) ? "IN" : "SC")) <<
+      " information for the face. Use: bopf* #f\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  char aText[32];
+  Standard_Integer nF, i1, i2, nV, i, aNb;
+  Handle(BOPDS_PaveBlock) aPB;
+  //
+  nF = Draw::Atoi(a[1]);
+  i1 = 0;
+  i2 = pDS->NbSourceShapes();
+  if (nF < i1 || nF > i2) {
+    di << " DS does not contain the shape\n";
+    return 1;
+  }
+  //
+  if (pDS->ShapeInfo(nF).ShapeType() != TopAbs_FACE) {
+    di << " The shape is not a face\n";
+    return 1;
+  }
+  //
+  if (!pDS->HasFaceInfo(nF)) {
+    di << " The face has no face information\n";
+    return 0;
+  }
+  //
+  BOPDS_FaceInfo& aFI = pDS->ChangeFaceInfo(nF);
+  //
+  BOPDS_IndexedMapOfPaveBlock aMPB;
+  BOPCol_MapOfInteger aMI;
+  if (iPriz == 0) {
+    strcpy(aText, "On");
+    aMPB = aFI.ChangePaveBlocksOn();
+    aMI = aFI.ChangeVerticesOn();
+  }
+  else if (iPriz == 1) {
+    strcpy(aText, "In");
+    aMPB = aFI.ChangePaveBlocksIn();
+    aMI = aFI.ChangeVerticesIn();
+  }
+  else if (iPriz == 2) {
+    strcpy(aText, "Sc");
+    aMPB = aFI.ChangePaveBlocksSc();
+    aMI = aFI.ChangeVerticesSc();
+  }
+  //
+  if (aMPB.Extent()) {
+    printf(" pave blocks %s:\n", aText);
+    aNb = aMPB.Extent();
+    for (i = 1; i <= aNb; ++i) {
+      aPB = aMPB(i);
+      aPB->Dump();
+      printf(" \n" );
+    }
+  }
+  else {
+    printf(" no pave blocks %s found\n", aText);
+  }
+  //
+  if (aMI.Extent()) {
+    printf(" vertices %s:\n", aText);
+    BOPCol_MapIteratorOfMapOfInteger aItMI(aMI);
+    for (; aItMI.More(); aItMI.Next()) {
+      nV = aItMI.Value();
+      printf(" %d", nV);
+    }
+    printf(" \n" );
+  }
+  else {
+    printf(" no verts %s found\n", aText);
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopfav
+//purpose  : alone vertices on face
+//=======================================================================
+Standard_Integer bopfav(Draw_Interpretor& di,
+                        Standard_Integer n,
+                        const char** a)
+{
+  if (n != 2) {
+    di << "Shows information about alone vertices for the face. Use: bopfav #f\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 0;
+  }
+  //
+  Standard_Integer i1, i2, nF, nV;
+  //
+  nF = Draw::Atoi(a[1]);
+  i1 = 0;
+  i2 = pDS->NbSourceShapes();
+  if (nF < i1 || nF > i2) {
+    di << "DS does not contain the shape\n";
+    return 1;
+  }
+  //
+  if (pDS->ShapeInfo(nF).ShapeType() != TopAbs_FACE) {
+    di << " The shape is not a face\n";
+    return 1;
+  }
+  //
+  if (!pDS->HasFaceInfo(nF)) {
+    di << " The face has no face information\n";
+    return 0;
+  }
+  //
+  BOPCol_ListOfInteger aLI;
+  pDS->AloneVertices(nF, aLI);
+  if (!aLI.Extent()) {
+    di << " no alone vertices found\n";
+    return 0;
+  }
+  //
+  di << " alone vertices: \n";
+  BOPCol_ListIteratorOfListOfInteger aItLI(aLI);
+  for (; aItLI.More(); aItLI.Next()) {
+    nV = aItLI.Value();
+    di << nV << " ";
+  }
+  di <<"\n";
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopimage
+//purpose  : 
+//=======================================================================
+Standard_Integer bopimage(Draw_Interpretor& di,
+                          Standard_Integer n,
+                          const char** a)
+{ 
+  if (n != 2) {
+    di << "Shows split parts of the shape. Use: bopimage s\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  TopoDS_Shape aS = DBRep::Get(a[1]);
+  if (aS.IsNull()) {
+    di << a[1] << " is a null shape\n";
+    return 1;
+  }
+  //
+  BOPAlgo_Builder& aBuilder = BOPTest_Objects::Builder();
+  const BOPCol_DataMapOfShapeListOfShape& anImages = aBuilder.Images();
+  if (!anImages.IsBound(aS)) {
+    di << " no images found\n"; 
+    return 0;
+  }
+  //
+  char buf[32];
+  Standard_Integer i;
+  BRep_Builder aBB;
+  TopoDS_Compound aC;
+  //
+  aBB.MakeCompound(aC);
+  //
+  const BOPCol_ListOfShape& aLSIm = anImages.Find(aS);
+  BOPCol_ListIteratorOfListOfShape aIt(aLSIm);
+  for (i = 0; aIt.More(); aIt.Next(), ++i) {
+    const TopoDS_Shape& aSIm = aIt.Value();
+    aBB.Add(aC, aSIm);
+  }
+  //
+  di << i << " images found\n";
+  sprintf(buf, "%s_im", a[1]);
+  DBRep::Set(buf, aC);
+  di << buf << "\n";
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : boporigin
+//purpose  : 
+//=======================================================================
+Standard_Integer boporigin(Draw_Interpretor& di,
+                           Standard_Integer n,
+                           const char** a)
+{ 
+  if (n != 2) {
+    di << "Shows the original shape for the shape. Use: boporigin s\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  TopoDS_Shape aS = DBRep::Get(a[1]);
+  if (aS.IsNull()) {
+    di << a[1] << " is a null shape\n";
+    return 0;
+  }
+  //
+  char buf[32];
+  //
+  BOPAlgo_Builder& aBuilder = BOPTest_Objects::Builder();
+  const BOPCol_DataMapOfShapeShape& aDMI = aBuilder.Origins();
+  if (!aDMI.IsBound(aS)) {
+    di << " no origins found\n"; 
+    return 0;
+  }
+  //
+  const TopoDS_Shape& aSx = aDMI.Find(aS);
+  //
+  sprintf(buf, "%s_or", a[1]);
+  DBRep::Set(buf, aSx);
+  //
+  di << buf << "\n";
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopfsd
+//purpose  : 
+//=======================================================================
+Standard_Integer bopfsd(Draw_Interpretor& di,
+                        Standard_Integer n,
+                        const char** a)
+{ 
+  if (n != 2) {
+    di << "Shows SD faces for the face: Use: bopfsd f\n";
+    return 1;
+  }
+  //
+  BOPDS_PDS pDS = BOPTest_Objects::PDS();
+  if (!pDS) {
+    di << " prepare PaveFiller first\n";
+    return 1;
+  }
+  //
+  TopoDS_Shape aS = DBRep::Get(a[1]);
+  if (aS.IsNull()) {
+    di << a[1] << " is a null shape\n";
+    return 1;
+  }
+  //
+  BOPAlgo_Builder& aBuilder = BOPTest_Objects::Builder();
+  const BOPCol_DataMapOfShapeShape& aDMSD = aBuilder.ShapesSD();
+  if (!aDMSD.IsBound(aS)) {
+    di << " shape has no sd shape\n"; 
+    return 0;
+  }
+  //
+  char buf[32];
+  Standard_Integer i;
+  BRep_Builder aBB;
+  TopoDS_Compound aC;
+  //
+  aBB.MakeCompound(aC);
+  //
+  BOPCol_DataMapIteratorOfDataMapOfShapeShape aItSD;
+  aItSD.Initialize(aDMSD);
+  for (i = 0; aItSD.More(); aItSD.Next()) {
+    const TopoDS_Shape& aSK = aItSD.Key();
+    const TopoDS_Shape& aSV = aItSD.Value();
+    if (aSK.IsEqual(aS)) {
+      if (!aSV.IsEqual(aS)) {
+        aBB.Add(aC, aS);
+        ++i;
+      }
+    }
+    //
+    else if (aSV.IsEqual(aS)) {
+      if (!aSK.IsEqual(aS)) {
+        aBB.Add(aC, aS);
+        ++i;
+      }
+    }
+  }
+  //
+  di << i << " SD shapes found\n";
+  //
+  sprintf(buf, "%s_sd", a[1]);
+  DBRep::Set(buf, aC);
+  //
+  di << buf << "\n";
+  return 0;
+}
+
+//=======================================================================
+//function : bopbface
+//purpose  : 
+//=======================================================================
+Standard_Integer bopbface (Draw_Interpretor& di,
+                           Standard_Integer n,
+                           const char** a)
+{
+  if ( n!= 3) {
+    di << "Build faces from surface and set of shared edges. Use: bopbface fr cx\n";
+    return 1;
+  }
+  //
+  TopoDS_Shape aS = DBRep::Get(a[2]);
+  if (aS.IsNull()) {
+    di << a[1] << " is a null shape\n";
+    return 1;
+  }
+  //
+  TopoDS_Face aF;
+  BOPCol_ListOfShape aLE;
+  Standard_Integer i;
+  //
+  TopoDS_Iterator aItS(aS);
+  for (i=0; aItS.More(); aItS.Next(), ++i) {
+    const TopoDS_Shape& aSx = aItS.Value();
+    if (!i) {
+      if (aSx.ShapeType() != TopAbs_FACE) {
+        di << " shape " << i << " is not a face\n";
+        return 1;
+      }
+      aF = *(TopoDS_Face*)&aSx;
+    }
+    else {
+      if (aSx.ShapeType() != TopAbs_EDGE) {
+        di << " shape " << i << " is not an edge\n";
+        return 1;
+      }
+      aLE.Append(aSx);
+    }
+  }
+  //
+  BOPAlgo_BuilderFace aBF;
+  aBF.SetFace(aF);
+  aBF.SetShapes(aLE);
+  aBF.Perform();
+  //
+  Standard_Integer iErr = aBF.ErrorStatus();
+  if (iErr != 0) {
+    di << " Error: " << iErr << "\n";
+    return 1;
+  }
+  //
+  char buf[128];
+  const BOPCol_ListOfShape& aLFR = aBF.Areas();
+  BOPCol_ListIteratorOfListOfShape aIt(aLFR);
+  for (i = 1; aIt.More(); aIt.Next(), ++i) {
+    const TopoDS_Shape& aFR = aIt.Value();
+    sprintf(buf, "%s_%d", a[1], i);
+    DBRep::Set(buf, aFR);
+    di << " " << buf;
+  }
+  //
+  i = aLFR.Extent();
+  if (i) {
+    di << "\n " << i << " faces were built\n";
+  }
+  else {
+    di << " No faces were built\n";
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : bopbsolid
+//purpose  : 
+//=======================================================================
+Standard_Integer bopbsolid (Draw_Interpretor& di,
+                            Standard_Integer n,
+                            const char** a)
+{
+  if ( n!= 3) {
+    di << "Build solids from set of shared faces. Use: bopbsolid sr cx\n";
+    return 1;
+  }
+  //
+  TopoDS_Shape aS = DBRep::Get(a[2]);
+  if (aS.IsNull()) {
+    di << a[1] << " is a null shape\n";
+    return 1;
+  }
+  //
+  BOPCol_ListOfShape aLF;
+  TopExp_Explorer aExp(aS, TopAbs_FACE);
+  for (; aExp.More(); aExp.Next()) {
+    const TopoDS_Shape& aF = aExp.Current();
+    aLF.Append(aF);
+  }
+  //
+  if (aLF.IsEmpty()) {
+    di << " No faces to build solids\n";
+    return 1;
+  }
+  //
+  BOPAlgo_BuilderSolid aBS;
+  aBS.SetShapes(aLF);
+  aBS.Perform();
+  //
+  Standard_Integer iErr = aBS.ErrorStatus();
+  if (iErr != 0) {
+    di << " Error: " << iErr << "\n";
+    return 1;
+  }
+  //
+  Standard_Integer i;
+  TopoDS_Compound aSolids;
+  BRep_Builder aBB;
+  //
+  aBB.MakeCompound(aSolids);
+  //
+  char buf[128];
+  const BOPCol_ListOfShape& aLSR = aBS.Areas();
+  BOPCol_ListIteratorOfListOfShape aIt(aLSR);
+  for (i = 1; aIt.More(); aIt.Next(), ++i) {
+    const TopoDS_Shape& aSR = aIt.Value();
+    sprintf(buf, "%s_%d", a[1], i);
+    DBRep::Set(buf, aSR);
+    di << " " << buf;
+  }
+  //
+  i = aLSR.Extent();
+  if (i) {
+    di << "\n " << i << " solids were built\n";
+  }
+  else {
+    di << " No solids were built\n";
+  }
+  //
+  return 0;
+}
+
+//=======================================================================
+//function : GetTypeByName
+//purpose  : 
+//=======================================================================
+void GetTypeByName(const char* theName,
+                   TopAbs_ShapeEnum& theType)
+{
+  if (!strcmp (theName, "v") ||
+      !strcmp (theName, "V")) {
+    theType = TopAbs_VERTEX;
+  }
+  else if (!strcmp (theName, "e") ||
+           !strcmp (theName, "E")) {
+    theType = TopAbs_EDGE;
+  }
+  else if (!strcmp (theName, "w") ||
+           !strcmp (theName, "W")) {
+    theType = TopAbs_WIRE;
+  }
+  else if (!strcmp (theName, "f") ||
+           !strcmp (theName, "F")) {
+    theType = TopAbs_FACE;
+  }
+  else if (!strcmp (theName, "sh") ||
+           !strcmp (theName, "Sh") ||
+           !strcmp (theName, "SH")) {
+    theType = TopAbs_SHELL;
+  }
+  else if (!strcmp (theName, "s") ||
+           !strcmp (theName, "S")) {
+    theType = TopAbs_SOLID;
+  }
+  else if (!strcmp (theName, "cs") ||
+           !strcmp (theName, "Cs") ||
+           !strcmp (theName, "CS")) {
+    theType = TopAbs_COMPSOLID;
+  }
+  else if (!strcmp (theName, "c") ||
+           !strcmp (theName, "C")) {
+    theType = TopAbs_COMPOUND;
+  }
+  else {
+    theType = TopAbs_SHAPE;
+  }
+}
+
+//=======================================================================
+//function : GetNameByType
+//purpose  : 
+//=======================================================================
+void GetNameByType(const TopAbs_ShapeEnum& theType,
+                   char* theName)
+{
+  switch (theType) {
+  case TopAbs_VERTEX:
+    strcpy (theName, "V");
+    break;
+  case TopAbs_EDGE:
+    strcpy (theName, "E");
+    break; 
+  case TopAbs_WIRE:
+    strcpy (theName, "w");
+    break;
+  case TopAbs_FACE:
+    strcpy (theName, "F");
+    break;
+  case TopAbs_SHELL:
+    strcpy (theName, "Sh");
+    break; 
+  case TopAbs_SOLID:
+    strcpy (theName, "S");
+    break; 
+  case TopAbs_COMPSOLID:
+    strcpy (theName, "Cs");
+    break; 
+  case TopAbs_COMPOUND:
+    strcpy (theName, "c");
+    break; 
+  default:
+    strcpy (theName, "Shape");
+    break; 
+  }
+}
+
+//=======================================================================
+//function : DumpInterfs
+//purpose  : 
+//=======================================================================
+template <class InterfType> void DumpInterfs
+  (const BOPCol_NCVector<InterfType>& theVInterf,
+   Draw_Interpretor& di)
+{
+  Standard_Integer i, aNb, n1, n2, nNew;
+  char buf[64];
+  //
+  aNb = theVInterf.Extent();
+  if (aNb == 0) {
+    di << "Not found\n";
+    return;
+  }
+  //
+  di << aNb << " interference(s) found\n";
+  for (i = 0; i < aNb; ++i) {
+    const InterfType& anInt = theVInterf(i);
+    anInt.Indices(n1, n2);
+    if (anInt.HasIndexNew()) {
+      nNew = anInt.IndexNew();
+      Sprintf(buf, " (%d, %d, %d)\n", n1, n2, nNew);
+    }
+    else {
+      Sprintf(buf, " (%d, %d)\n", n1, n2);
+    }
+    di << buf;
+  }
+}
+
+//=======================================================================
+//function : SearchNewIndex
+//purpose  : 
+//=======================================================================
+template <class InterfType> void SearchNewIndex
+  (const char* theCType,
+   const Standard_Integer theInd,
+   const BOPCol_NCVector<InterfType>& theVInterf,
+   Draw_Interpretor& di)
+{
+  char buf[64];
+  Standard_Boolean bFound;
+  Standard_Integer i, aNb, n1, n2, nNew;
+  //
+  bFound = Standard_False;
+  aNb = theVInterf.Extent();
+  for (i = 0 ; i < aNb; ++i) {
+    const InterfType& anInt = theVInterf(i);
+    nNew = anInt.IndexNew();
+    if (theInd == nNew) {
+      if (!bFound) {
+        di << theCType;
+        bFound = Standard_True;
+      }
+      //
+      anInt.Indices(n1, n2);
+      sprintf(buf,"(%d, %d) ", n1, n2);
+      di << buf;
+    }
+  }
+  if (bFound) {
+    di << "\n";
+  }
+}
diff --git a/src/BOPTest/FILES b/src/BOPTest/FILES
index d3aeb3e4ef..3c62a01a28 100755
--- a/src/BOPTest/FILES
+++ b/src/BOPTest/FILES
@@ -13,4 +13,4 @@ BOPTest_OptionCommands.cxx
 BOPTest_PartitionCommands.cxx
 BOPTest_TolerCommands.cxx
 BOPTest_HistoryCommands.cxx
-BOPTest_CellsCommands.cxx
\ No newline at end of file
+BOPTest_DebugCommands.cxxBOPTest_CellsCommands.cxx
\ No newline at end of file
diff --git a/tests/bugs/modalg_2/bug23008 b/tests/bugs/modalg_2/bug23008
index 8b3c7251c7..c9ff6b2de0 100755
--- a/tests/bugs/modalg_2/bug23008
+++ b/tests/bugs/modalg_2/bug23008
@@ -10,8 +10,10 @@ set BugNumber OCC23008
 
 restore [locate_data_file bug23008_ez15] b1 
 restore [locate_data_file bug23008_fz37] b2 
-bop b1 b2
-set info [ bopnews -v ]
+bclearobjects; bcleartools
+baddobjects b1; baddtools b2
+bfillds
+set info [ bopnews v ]
 for {set i 0} {$i <= 100} {incr i} {
     set word [lsearch $info z$i]
     if {$word == -1} {