CPPLapack Code

--- a/trunk/include/cpplapack.h
+++ b/trunk/include/cpplapack.h
@@ -29,7 +29,7 @@
 #endif//CPPL_EPS
 
 #ifndef CPPL_SS_SECTOR
-#define CPPL_SS_SECTOR 100
+#define CPPL_SS_SECTOR 10000
 #endif//CPPL_SS_SECTOR
 
 //=============================================================================

--- a/trunk/include/dssmatrix-/dssmatrix-io.hpp
+++ b/trunk/include/dssmatrix-/dssmatrix-io.hpp
@@ -45,12 +45,13 @@
   }
 #endif//CPPL_DEBUG
   
+  //// search (i,j) component ////
   for(long c=0; c<VOL; c++){
     if(Indx[c]==i && Jndx[c]==j){ return Array[c]; }
   }
   
   //// (i,j) component not found ////
-  if(VOL==CAP){ 
+  if(VOL==CAP){
     expand(CPPL_SS_SECTOR);
 #ifdef  CPPL_DEBUG
     std::cerr << "[NOTE] dssmatrix::put(long&, long&, double&) "
@@ -68,7 +69,7 @@
 ///////////////////////////////////////////////////////////////////////////////
 
 //=============================================================================
-/*! put value with isListed check and volume cheack */
+/*! put value with volume cheack without isListed check */
 inline void dssmatrix::put(const long& i, const long& j, const double& v)
 {
 #ifdef  CPPL_VERBOSE
@@ -85,11 +86,6 @@
     exit(1);
   }
 #endif//CPPL_DEBUG
-  
-  //// in case (i,j) already exists ////
-  for(long c=0; c<VOL; c++){
-    if(Indx[c]==i && Jndx[c]==j){ Array[c]=v; return; }
-  }
   
   //// in case (i,j) doesnot exist ////
   if(VOL==CAP){