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; }