LinkedList tailadd(LinkedList l1, LinkedList l2){
if(l1==NULL)
return l2;
if(l2==NULL)
return l1;
int re=0,tsum=0;
LinkedList thead(0);
LinkedList pre=thead;
while(l1!=NULL && l2!=NULL){
tsum=re+l1.val+l2.val;
re=tsum/10;
tsum=tsum%10;
LinkedList temp(tsum);
pre.next=temp;
pre=temp;
l1=l1.next;
l2=l2.next;
}
LinkedList rest=NULL;
if (l1!=NULL)
rest=l1;
if(l2!=NULL )
rest=l2;
while(res!=NULL){
tsum=re+res.val;
re=tsum/10;
tsum=tsum%10;
LinkedList temp(tsum);
pre.next=temp;
pre=temp;
res=res.next;
}
if(re!=0){
LinkedList temp(re);
pre.next=temp;
}
return thead.next;
}
int Delta::execute()
{
Options sourceOptions;
{
sourceOptions.add<std::string>("filename", m_sourceFile);
sourceOptions.add<bool>("debug", isDebug());
sourceOptions.add<boost::uint32_t>("verbose", getVerboseLevel());
}
Stage* source = AppSupport::makeReader(sourceOptions);
source->initialize();
boost::uint32_t totalPointCount(source->getNumPoints());
PointBuffer source_data(source->getSchema(), totalPointCount);
StageSequentialIterator* source_iter = source->createSequentialIterator(source_data);
boost::uint32_t numRead = source_iter->read(source_data);
assert(numRead == source_data.getNumPoints());
delete source_iter;
delete source;
Options candidateOptions;
{
candidateOptions.add<std::string>("filename", m_candidateFile);
candidateOptions.add<bool>("debug", isDebug());
candidateOptions.add<boost::uint32_t>("verbose", getVerboseLevel());
}
Stage* candidate = AppSupport::makeReader(candidateOptions);
candidate->initialize();
IndexedPointBuffer candidate_data(candidate->getSchema(), totalPointCount);
StageSequentialIterator* candidate_iter = candidate->createSequentialIterator(candidate_data);
numRead = candidate_iter->read(candidate_data);
assert(numRead == candidate_data.getNumPoints());
delete candidate_iter;
if (source_data.getNumPoints() != candidate_data.getNumPoints())
{
std::cerr << "Source and candidate files do not have the same point count, testing each source point only!" << std::endl;
}
// m_summary_x(xd);
// m_summary_y(yd);
// m_summary_z(zd);
if (m_outputFileName.size())
{
m_outputStream = FileUtils::createFile(m_outputFileName);
}
candidate_data.build(m_3d);
boost::uint32_t count(std::min(source_data.getNumPoints(), candidate_data.getNumPoints()));
boost::scoped_ptr<std::map<Point, Point> > points(cumulatePoints(source_data, candidate_data));
if (m_OutputDetail)
{
outputDetail(source_data, candidate_data, points.get());
return 0;
}
std::map<Point, Point>::const_iterator i;
for(i = points->begin(); i != points->end(); ++i)
{
Point const& s = i->first;
Point const& c = i->second;
double xd = s.x - c.x;
double yd = s.y - c.y;
double zd = s.z - c.z;
m_summary_x(xd);
m_summary_y(yd);
m_summary_z(zd);
}
std::string headline("------------------------------------------------------------------------------------------");
std::cout << headline << std::endl;
std::cout << " Delta summary for source '" << m_sourceFile << "' and candidate '" << m_candidateFile <<"'" << std::endl;
std::cout << headline << std::endl;
std::cout << std::endl;
std::string thead("----------- --------------- --------------- --------------");
std::cout << thead << std::endl;
std::cout << " Dimension X Y Z " << std::endl;
std::cout << thead << std::endl;
boost::format fmt("%.4f");
double sminx = (boost::accumulators::min)(m_summary_x);
double sminy = (boost::accumulators::min)(m_summary_y);
double sminz = (boost::accumulators::min)(m_summary_z);
double smaxx = (boost::accumulators::max)(m_summary_x);
double smaxy = (boost::accumulators::max)(m_summary_y);
double smaxz = (boost::accumulators::max)(m_summary_z);
double smeanx = (boost::accumulators::mean)(m_summary_x);
double smeany = (boost::accumulators::mean)(m_summary_y);
double smeanz = (boost::accumulators::mean)(m_summary_z);
std::cout << " Min " << fmt % sminx << " " << fmt % sminy << " " << fmt % sminz<<std::endl;
std::cout << " Min " << fmt % smaxx << " " << fmt % smaxy << " " << fmt % smaxz<<std::endl;
std::cout << " Mean " << fmt % smeanx << " " << fmt % smeany << " " << fmt % smeanz<<std::endl;
std::cout << thead << std::endl;
return 0;
}
