Instruct: Identify the main and secondary category of Arxiv papers based on the titles and abstracts
Query: Fast geodesic shooting for landmark matching using CUDA   Landmark matching via geodesic shooting is a prerequisite task for numerous
registration based applications in biomedicine. Geodesic shooting has been
developed as one solution approach and formulates the diffeomorphic
registration as an optimal control problem under the Hamiltonian framework. In
this framework, with landmark positions q0 fixed, the problem solely depends on
the initial momentum p0 and evolves through time steps according to a set of
constraint equations. Given an initial p0, the algorithm flows q and p forward
through time steps, calculates a loss based on point-set mismatch and kinetic
energy, back-propagate through time to calculate gradient on p0 and update it.
In the forward and backward pass, a pair-wise kernel on landmark points K and
additional intermediate terms have to be calculated and marginalized, leading
to O(N2) computational complexity, N being the number of points to be
registered. For medical image applications, N maybe in the range of thousands,
rendering this operation computationally expensive. In this work we ropose a
CUDA implementation based on shared memory reduction. Our implementation
achieves nearly 2 orders magnitude speed up compared to a naive CPU-based
implementation, in addition to improved numerical accuracy as well as better
registration results.

26
The categories of this paper are: bayes-an
8
The categories of this paper are: hep-ex
16
The categories of this paper are: stat
11
The categories of this paper are: hep-lat
5
The categories of this paper are: gr-qc
18
The categories of this paper are: econ
20
The categories of this paper are: acc-phys
14
The categories of this paper are: math-ph
25
The categories of this paper are: patt-sol
21
The categories of this paper are: adap-org
7
The categories of this paper are: hep-th
12
The categories of this paper are: nucl-th
30
The categories of this paper are: q-alg
19
The categories of this paper are: eess
23
The categories of this paper are: solv-int
29
The categories of this paper are: funct-an