 details about this method will be discussed in Section \[sec:SolvePDEs\] and will be verified in our numerical experiments. Furthermore, this approach enables us to conduct geometric understanding of data without global coordinates reconstruction. Examples include global pattern extraction, comparison and classification as many existing methods conducted using results of differential equations [@Shi:08a; @Lai:2010CVPR; @lai2017multi]. As a byproduct of global information from PDEs, we also propose a new method for reconstructing manifolds by stitching its local patches. This new method is much more efficient than the way of direct global reconstruction using matrix completion method. Moreover, it can also overcome possible reconstruction failure using global matrix completion due to the coherent missing information.

The rest of this paper is organized as follows. In section \[sec:CoordinateRecon\], we propose a low-rank matrix completion model to reconstruct local coordinates and design an algorithm to solve the proposed convex optimization problem based on operator splitting and the alternating direction method. After that, section \[sec:SolvePDEs\] is devoted to discuss point-wisely approximating differential operators based on the intrinsic methods proposed in [@lai2013local; @Liang:CVPR2012; @liang2013solving]. We also discuss our two model PDEs, the Laplace-Beltrami (LB) eigenvalue problem and the Eiknoal equation for solve geodesic distance on manifolds. As an application of global information using solutions of LB eigenvalue problem, we propose a new method for reconstructing manifolds by stitching local patches in the LB frequency space in section \[sec:MfdStitching\]. All numerical experiments about the effectiveness of our local coordinates reconstruction model, the accuracy and robustness of the PDEs solvers, and demonstration of manifold reconstruction are discussed in section \[sec:experiments\]. Finally, we conclude our work in section \[sec:conclusion\].

Local coordinates reconstruction via matrix completion {#sec:CoordinateRecon}
======================================================

In this section, we first review the concept of the classical multidimensional scaling (MDS). Inspired by the classical MDS, we propose a method of reconstructing local coordinates from the given incomplete inter-point distance matrix using the low-rank matrix completion [@CandesRe2008]. A numerical algorithm based on the augmented Lagrangian is also designed to solve the proposed convex optimization problem. As we described before, our strategy is to conduct local coordinates reconstruction for the KNN $
details about this method will be discussed in Section   \[sec: SolvePDEs\ ] and will be verified in our numeric experiment. Furthermore, this approach enables us to conduct geometric sympathy of data without global coordinate reconstruction. Examples include global form extraction, comparison and classification as many exist methods conducted use result of differential equations   [ @Shi:08a; @Lai:2010CVPR; @lai2017multi ]. As a byproduct of global information from PDEs, we also project a new method for reconstructing manifold by stitching its local patches. This new method acting is much more effective than the way of direct global reconstruction using matrix completion method acting. Moreover, it can also overcome possible reconstruction failure using global matrix completion due to the coherent missing information. 

 The rest of this paper is organized as follows. In section   \[sec: CoordinateRecon\ ], we propose a low - rank matrix completion exemplar to reconstruct local coordinates and plan an algorithm to resolve the aim convex optimization problem based on operator splitting and the alternating direction method. After that, section   \[sec: SolvePDEs\ ] is devoted to discourse point - wisely approximating differential operators based on the intrinsic methods proposed in   [ @lai2013local; @Liang: CVPR2012; @liang2013solving ]. We besides discuss our two model PDEs, the Laplace - Beltrami (LB) eigenvalue problem and the Eiknoal equation for solve geodesic distance on manifolds. As an lotion of ball-shaped data using solutions of LB eigenvalue problem, we propose a new method for reconstructing manifolds by sew local patches in the LB frequency distance in section   \[sec: MfdStitching\ ]. All numerical experiment about the effectiveness of our local coordinates reconstruction exemplar, the accuracy and robustness of the PDEs solvers, and presentation of manifold reconstruction are discussed in section   \[sec: experiments\ ]. Finally, we conclude our oeuvre in section   \[sec: conclusion\ ]. 

 Local coordinates reconstruction via matrix completion { # sec: CoordinateRecon } 
 = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = = 

 In this section, we first review the concept of the classical multidimensional scaling (MDS). inspire by the classical MDS, we aim a method of reconstructing local coordinates from the given incomplete inter - point distance matrix use the low - rank matrix completion   [ @CandesRe2008 ]. A numeral algorithm based on the augmented Lagrangian is also design to solve the proposed convex optimization problem. As we identify before, our strategy is to conduct local coordinates reconstruction for the KNN $
 dehails about this method dill be discussgd in Sertion \[sed:SolvePDDs\] and will be verified in onr nymerixal experiments. Furthefmore, thid approaxh eiables us to convhct geometric mndervvanding of data without gnobal coordinader xeconstruction. Examples include globwl pattrrj extraction, cjmpagifon znd classification as many existinf methovs conducted usong results of differentiap eqkations [@Shi:08a; @Lai:2010CVPG; @lai2017multi]. Qs a vyproduct of global information frkm PDEs, we also propose a new mdthod for reconwteuchhng manifolvs by ftitching its local [atches. This new methpd ms mych more efficient then the way of direct global racknstruction using mqtrix com[letkin oetgov. Mkreoveg, iv can also kvercome powsible reconstructipn dailure using globaj iatrix completion due to the coherent mpssihg information.

The rest if this paper is orgajized as sollows. In section \[sec:CoordinateRecon\], we propose a low-reny mcbvix zimoletion model to reconstruct local coordinatef amd design an algjrithm to splge jhe proposed cunvex optjmization problem hased og opeeator splyttimg and the alternating dirextion method. Qfter that, section \[dec:SolvePDEr\] is devoyed to discuss point-wiszly aplroximating differenfkal operators bared ot the inufinsic methods prjposed in [@oai2013lpcal; @Lixng:CFPR2012; @liwng2013solving]. We also discuss our two modep [DEs, the Lwplace-Beltrami (LB) eigenvalue problem and the Giktoan equatijn fov solve geodesis distance on kanifolbs. As xn applicanion of gnobal inforiation using vllutions of NB eigendaluw priblem, wd propose a new method for reconsteucting manifolds ny stkfching local pacehts in the LB freauevcy siacx in faction \[sec:MfdVtitzhivb\]. All numerical cxpdrimrnts about the effecdivehess of our local vojrdinatew reconseruction modek, the accuracy and robuvtnxss of thg PDEs solvers, and demonstratioh of manivolq reconstruceion are discusszd in section \[sec:experiments\]. Finally, we cinclude our work ib section \[sec:conclusnok\].

Local coordmnates reconstrgction via matrix conpletion {#sec:CoordlnateRecon}
======================================================

In this sectjon, we firsh review the concept of the classical multidimensional scaling (MDS). Inspirwd by vhq classical MDS, we pgopjxe a mqthov uf reconstructinb local coordinates from the giten incompnece inter-point distance matrix uring the low-rxnk matrix completion [@CansesRe2008]. A mumerical algorithm based on thr augmented Lagrangian ix also dwsignsd to solvx the proposed convec optikizatuon provlem. As we descriced bcfore, our dtrauegy is to conduct local coordinates reconwtguction for the NNN $
details about this method will be discussed \[sec:SolvePDEs\] will be in our numerical us conduct geometric understanding data without global reconstruction. Examples include global pattern extraction, and classification as many existing methods conducted using results of differential equations [@Shi:08a; @lai2017multi]. As a byproduct of global information from PDEs, we also propose a method reconstructing by its local patches. This new method is much more efficient than the way of direct global using matrix completion method. Moreover, it can also possible reconstruction failure using matrix completion due to the missing The rest this is as follows. In \[sec:CoordinateRecon\], we propose a low-rank matrix completion model to reconstruct local coordinates and design an algorithm to the proposed problem based operator and alternating direction method. section \[sec:SolvePDEs\] is devoted to discuss operators based on the intrinsic methods proposed in @Liang:CVPR2012; @liang2013solving]. also discuss our two model PDEs, Laplace-Beltrami (LB) eigenvalue problem and the Eiknoal equation solve geodesic distance on manifolds. As an application of global information using solutions of LB we propose a new for reconstructing manifolds stitching patches the frequency space section \[sec:MfdStitching\]. All numerical experiments about the effectiveness of our local reconstruction model, the accuracy and robustness of the PDEs solvers, of reconstruction are discussed section \[sec:experiments\]. Finally, we our in section \[sec:conclusion\]. Local via completion this we review the concept of classical multidimensional scaling (MDS). Inspired the classical MDS, we local coordinates from the given incomplete inter-point distance using the low-rank matrix completion [@CandesRe2008]. A algorithm based on the augmented Lagrangian is also designed to solve the convex optimization we described before, our strategy is to conduct coordinates reconstruction for the $
 details about this method wilL be discussEd in SEctIon \[SeC:SolVePDes\] and will be verIFied In our numerical experimeNts. FuRtHErmoRE, tHis apProach eNAbLES us To CoNduCt GEoMetriC unDerstanDing of data WitHoUt global coorDInAtes reconsTruCtion. ExampleS inClude gLoBal PAtterN exTractIon, comPArison And classiFiCAtion aS Many exiSTInG metHods conducted usinG ReSUlts of differenTial eqUaTIoNS [@shi:08A; @LaI:2010CVPR; @lai2017muLtI]. As a bYProduct OF gLOBAl iNFormation from pDEs, we also pROpoSe a new MeThoD For recOnstrUcTIng Manifolds by StitChing its lOcal paTChes. ThiS New methOd is muCh mOre EffiCIeNt ThaN tHE waY Of DirECt gLobal recOnStRuctiOn usING MAtriX coMpleTion mEthod. Moreover, It cAn alSO ovErcomE possIble ReConstRuctioN failUrE using global matRix cOmpletion Due To The CoHerenT MissinG inForMation.

THe rest oF ThiS pAPER iS organized as followS. IN SEcTion \[sec:COordinATereCOn\], we propOsE a lOw-raNK MatriX comPLeTion modeL to recONsTrUct locaL cOordinAtEs aNd dEsign AN algOrithm To solve tHe proPOsed convex optiMIzation probleM BaSED oN OperAtoR splitting aNd thE AlteRnatINg DirECtion MethoD. AFTeR That, section \[sec:SolvepDes\] is deVoted To discuss poinT-wisely appROXImating dIffeREnTIal operators baSed on The intrinsIC methods PropoSed in [@lai2013Local; @LianG:cvPR2012; @liang2013SolVinG]. We AlsO DIsCuss our two modEL pDEs, ThE LaplacE-BeLtrami (Lb) eiGenValUe pRoBlem and thE Eiknoal EqUaTiOn For Solve GEodesic dIsTanCe On mAnifoLDs. As an AppliCatiOn Of GLobAl inforMAtION usiNg SoLutiOns Of lB eigEnvaLUe pRoblem, wE propose a New MEthoD fOr ReconstRucting manifoLdS by stitchiNg LocAl patcHES in the LB Frequency space in section \[SEc:MfdStItcHing\]. ALl nuMerical exPerIments AboUT the efFectivEness Of Our LOCal coORDiNatEs ReconstrucTIOn mOdel, tHe AccuRacy and Robustness of the PDES SolVers, and demonsTraTion OF MaNifOLd REcoNsTRucTIOn are discussed iN section \[seC:eXPeRiments\]. FinALly, We ConcludE our worK in seCTion \[sec:ConclusioN\].

Local cooRdInatES RecOnstructioN via matrIx completIOn {#sec:cOoRdinaTeREcon}
======================================================

In ThIs sEctioN, we firST reView tHe concEpT of the ClassIcAl multidImensional scaling (MDS). InsPired bY the cLasSical MDS, wE prOPosE a method oF recOnstructinG loCal CoordInaTEs froM the GIvEn iNComplEte iNTer-point dIStAncE MAtRix using the LOW-RanK matrIx cOMpletiOn [@CaNdesRe2008]. A numerical aLGorithm based on The aUGMenTed lAgraNgIan is also desigNed To SOLve the prOpOsed convex oPtimizatIoN ProblEm. As we DescriBed befoRE, OuR StrateGy is To cOnduct locAl cOoRDinates ReCoNStructIon fOr The KNN $
 details about this method  will be d iscus sed  in  S ecti on \ [sec:SolvePDEs \ ] an d will be verified inour n um e rica l  e xperi ments.F ur t h erm or e,  th is ap proac h e nablesus to cond uct  g eometric und e rs tanding of  da ta without g lob al coo rd ina t es re con struc tion.E xample s include  g l obal p a ttern e x t ra ctio n, comparison and cl a ssification as  manyex i st i n g m eth ods conduc te d usi n g resul t so f dif f erential equa tions [@Shi : 08a ; @Lai :2 010 C VPR; @ lai20 17 m ult i]. As a by prod uct of gl obal i n formati o n fromPDEs,weals o pr o po se  ane w  me t ho d f o r r econstru ct in g man ifol d s b y st itc hing  itslocal patches . T hisn ewmetho d ismuch  m ore e fficie nt th an  the way of dir ectglobal re con st ruc ti on us i ng mat rix  co mpletio n metho d . M or e o v er , it can also over co m e  p ossiblerecons t ru ct i on failu re  us ingg l obalmatr i xcompleti on due to  t he cohe re nt mis si nginf ormat i on.
The r est of t his p a per is organiz e d as follows. In s ec t ion\[s ec:Coordina teRe c on\] , we pr opo s e a l ow-ra nk ma t rix completion mode lto rec onstr uct local coo rdinates a n d design a n al g or i thm to solve t he pr oposed con v ex optim izati on probl em basedo n  operato r s pli tti nga n dthe alternati n g  dir ec tion me tho d. Afte r t hat , s ect io n \[sec:S olvePDEs \]  i sde vot ed to discusspo int -w ise ly ap p roxima tingdiff er en t ial  operat o rs b ased  o ntheint ri nsicmeth o dspropose d in [@la i20 1 3loc al ;@Liang: CVPR2012; @li an g2013solvi ng ].We als o discussour two model PDEs, the Laplace -Be ltram i (L B) eigenv alu e prob lem and th e Eikn oal e qu ati o n  fors o lv e g eo desic dist a n ceon ma ni fold s. As a n application of g l oba l information  us ings o lu tio n so f L Be ige n v alue problem, w e proposean ew  method fo r  re co nstruct ing man ifold s  by sti tching lo cal patch es  int h e L B frequenc y spacein sectio n  \[se c :M fdSti tch ing\].  A llnumer ical e x per iment s abou tthe ef fecti ve ness ofour local coordinates r econst ructi onmodel, th e a c cur acy and r obus tness of t hePDE s sol ver s , and  dem o ns tra t ion o f ma n ifold rec o ns tru c t io n are discu s s e d i n sec tio n  \[sec :exp eriments\]. Final l y, we conclude  our w ork  in sect io n \[sec:conclu sio n\ ] . 

Localco ordinates r econstru ct i on vi a matr ix com pletion { #s e c:Coor dina teR econ}
=== === == = ======= == == = ====== ==== == ====== ====== = ==== = = ==

In this sect ion,w e  firs t  re viewth e conce p t of  the class ical multid imensi onal  scal ing (MD S) . Insp ire dby the cla s sical MDS , wepropose  a  met hod  of re cons t r uctin g lo ca l c oordinate s fr o mth e  gi venincom pl eteinter-poi n t distan cem atrix u si ngt h e low- r an k matrix com ple tion[ @ CandesRe20 0 8].A  n u meric al alg orithm  basedo n t he  augmen ted L agrangian  is alsod esi gn ed t o solveth e pr op ose dconvex optimiza tion pro blem. As w e desc rib e d be f ore, o ur  strategy is  to c on d uct loc a l  co ordinat e s r eco nst ruct ion  for the KNN $
 details_about this_method will be discussed_in Section \[sec:SolvePDEs\]_and_will be_verified_in our numerical_experiments. Furthermore, this_approach enables us to_conduct geometric understanding_of_data without global coordinates reconstruction. Examples include global pattern extraction, comparison and classification as_many_existing methods_conducted_using_results of differential equations [@Shi:08a; @Lai:2010CVPR;_@lai2017multi]. As a byproduct of_global information_from PDEs, we also propose a new method_for_reconstructing manifolds by_stitching its local patches. This new method is much_more efficient than the way of_direct global reconstruction_using_matrix_completion method. Moreover, it_can also overcome possible reconstruction failure_using global matrix completion due to_the coherent missing information.

The rest of this_paper is organized as follows. In_section \[sec:CoordinateRecon\], we propose a low-rank_matrix completion_model to reconstruct local coordinates_and design an_algorithm to_solve the proposed_convex optimization problem based on operator_splitting and the_alternating direction method. After that, section \[sec:SolvePDEs\]_is_devoted to discuss_point-wisely_approximating_differential operators_based on the_intrinsic_methods proposed_in [@lai2013local;_@Liang:CVPR2012; @liang2013solving]. We also discuss our_two_model PDEs, the Laplace-Beltrami (LB) eigenvalue problem_and the Eiknoal equation_for_solve geodesic distance on_manifolds. As an application of_global information using solutions of LB_eigenvalue problem,_we propose_a new method for reconstructing manifolds by stitching local patches in_the LB frequency space in section \[sec:MfdStitching\]._All numerical experiments about_the effectiveness_of_our local coordinates_reconstruction_model, the_accuracy and robustness of the PDEs solvers,_and demonstration_of manifold reconstruction are discussed in_section \[sec:experiments\]. Finally, we conclude_our_work in section \[sec:conclusion\].

Local coordinates reconstruction via_matrix completion {#sec:CoordinateRecon}
======================================================

In this section, we_first review the concept of_the_classical_multidimensional scaling (MDS). Inspired by_the classical MDS, we propose a_method of reconstructing_local coordinates from the given incomplete inter-point_distance_matrix using the low-rank matrix completion [@CandesRe2008]._A_numerical algorithm based on the augmented_Lagrangian_is_also designed to solve the_proposed convex optimization problem. As we_described before, our strategy is to conduct local coordinates_reconstruction for the_KNN $