
We study the noncommutative Poincar\'e duality between the Poisson homology
and cohomology of unimodular Poisson algebras, and show that Kontsevich's
deformation quantization as well as Koszul duality preserve the corresponding
Poincar\'e duality. As a corollary, the BatalinVilkovisky algebra structures
that naturally arise in these cases are all isomorphic.

In this paper we consider a general matrix factorization model which covers a
large class of existing models with many applications in areas such as machine
learning and imaging sciences. To solve this possibly nonconvex, nonsmooth and
nonLipschitz problem, we develop a nonmonotone alternating updating method
based on a potential function. Our method essentially updates two blocks of
variables in turn by inexactly minimizing this potential function, and updates
another auxiliary block of variables using an explicit formula. The special
structure of our potential function allows us to take advantage of efficient
computational strategies for nonnegative matrix factorization to perform the
alternating minimization over the two blocks of variables. A suitable line
search criterion is also incorporated to improve the numerical performance.
Under some mild conditions, we show that the line search criterion is well
defined, and establish that the sequence generated is bounded and any cluster
point of the sequence is a stationary point. Finally, we conduct some numerical
experiments using real datasets to compare our method with some existing
efficient methods for nonnegative matrix factorization and matrix completion.
The numerical results show that our method can outperform these methods for
these specific applications.

Movie recommendation systems provide users with ranked lists of movies based
on individual's preferences and constraints. Two types of models are commonly
used to generate ranking results: longterm models and sessionbased models.
While longterm models represent the interactions between users and movies that
are supposed to change slowly across time, sessionbased models encode the
information of users' interests and changing dynamics of movies' attributes in
short terms. In this paper, we propose an LSIC model, leveraging Long and
Shortterm Information in Contentaware movie recommendation using adversarial
training. In the adversarial process, we train a generator as an agent of
reinforcement learning which recommends the next movie to a user sequentially.
We also train a discriminator which attempts to distinguish the generated list
of movies from the real records. The poster information of movies is integrated
to further improve the performance of movie recommendation, which is
specifically essential when few ratings are available. The experiments
demonstrate that the proposed model has robust superiority over competitors and
sets the stateoftheart. We will release the source code of this work after
publication.

The manual outlining of hepatic metastasis in (US) ultrasound acquisitions
from patients suffering from pancreatic cancer is common practice. However,
such pure manual measurements are often very time consuming, and the results
repeatedly differ between the raters. In this contribution, we study the
indepth assessment of an interactive graphbased approach for the segmentation
for pancreatic metastasis in US images of the liver with two specialists in
Internal Medicine. Thereby, evaluating the approach with over one hundred
different acquisitions of metastases. The two physicians or the algorithm had
never assessed the acquisitions before the evaluation. In summary, the
physicians first performed a pure manual outlining followed by an algorithmic
segmentation over one month later. As a result, the experts satisfied in up to
ninety percent of algorithmic segmentation results. Furthermore, the
algorithmic segmentation was much faster than manual outlining and achieved a
median Dice Similarity Coefficient (DSC) of over eighty percent. Ultimately,
the algorithm enables a fast and accurate segmentation of liver metastasis in
clinical US images, which can support the manual outlining in daily practice.

In this paper we show that for a Koszul CalabiYau algebra, there is a
shifted bisymplectic structure in the sense of
CrawleyBoeveyEtingofGinzburg, on the cobar construction of its counitalized
Koszul dual coalgebra, and hence its DG representation schemes, in the sense of
BerestKhachatryanRamadoss, have a shifted symplectic structure in the sense
of PantevTo\"enVaqui\'eVezzosi.

A solution of twostage stochastic generalized equations is a pair: a first
stage solution which is independent of realization of the random data and a
second stage solution which is a function of random variables.This paper
studies convergence of the sample average approximation of twostage stochastic
nonlinear generalized equations. In particular an exponential rate of the
convergence is shown by using the perturbed partial linearization of functions.
Moreover, sufficient conditions for the existence, uniqueness, continuity and
regularity of solutions of twostage stochastic generalized equations are
presented under an assumption of monotonicity of the involved functions. These
theoretical results are given without assuming relatively complete recourse,
and are illustrated by twostage stochastic noncooperative games of two
players.

In this contribution, a software system for computeraided position planning
of miniplates to treat facial bone defects is proposed. The intraoperatively
used bone plates have to be passively adapted on the underlying bone contours
for adequate bone fragment stabilization. However, this procedure can lead to
frequent intraoperatively performed material readjustments especially in
complex surgical cases. Our approach is able to fit a selection of common
implant models on the surgeon's desired position in a 3D computer model. This
happens with respect to the surrounding anatomical structures, always including
the possibility of adjusting both the direction and the position of the used
osteosynthesis material. By using the proposed software, surgeons are able to
preplan the out coming implant in its form and morphology with the aid of a
computervisualized model within a few minutes. Further, the resulting model
can be stored in STL file format, the commonly used format for 3D printing.
Using this technology, surgeons are able to print the virtual generated
implant, or create an individually designed bending tool. This method leads to
adapted osteosynthesis materials according to the surrounding anatomy and
requires further a minimum amount of money and time.

Patientspecific cranial implants are important and necessary in the surgery
of cranial defect restoration. However, traditional methods of manual design of
cranial implants are complicated and timeconsuming. Our purpose is to develop
a novel software named EasyCrania to design the cranial implants conveniently
and efficiently. The process can be divided into five steps, which are
mirroring model, clipping surface, surface fitting, the generation of the
initial implant and the generation of the final implant. The main concept of
our method is to use the geometry information of the mirrored model as the base
to generate the final implant. The comparative studies demonstrated that the
EasyCrania can improve the efficiency of cranial implant design significantly.
And, the intra and interrater reliability of the software were stable, which
were 87.07+/1.6% and 87.73+/1.4% respectively.

We consider a class of differenceofconvex (DC) optimization problems whose
objective is levelbounded and is the sum of a smooth convex function with
Lipschitz gradient, a proper closed convex function and a continuous concave
function. While this kind of problems can be solved by the classical
differenceofconvex algorithm (DCA) [26], the difficulty of the subproblems of
this algorithm depends heavily on the choice of DC decomposition. Simpler
subproblems can be obtained by using a specific DC decomposition described in
[27]. This decomposition has been proposed in numerous work such as [18], and
we refer to the resulting DCA as the proximal DCA. Although the subproblems are
simpler, the proximal DCA is the same as the proximal gradient algorithm when
the concave part of the objective is void, and hence is potentially slow in
practice. In this paper, motivated by the extrapolation techniques for
accelerating the proximal gradient algorithm in the convex settings, we
consider a proximal differenceofconvex algorithm with extrapolation to
possibly accelerate the proximal DCA. We show that any cluster point of the
sequence generated by our algorithm is a stationary point of the DC
optimization problem for a fairly general choice of extrapolation parameters:
in particular, the parameters can be chosen as in FISTA with fixed restart
[15]. In addition, by assuming the Kurdyka{\L}ojasiewicz property of the
objective and the differentiability of the concave part, we establish global
convergence of the sequence generated by our algorithm and analyze its
convergence rate. Our numerical experiments on two differenceofconvex
regularized least squares models show that our algorithm usually outperforms
the proximal DCA and the general iterative shrinkage and thresholding algorithm
proposed in [17].

In this paper, we propose a discretization scheme for the twostage
stochastic linear complementarity problem (LCP) where the underlying random
data are continuously distributed. Under some moderate conditions, we derive
qualitative and quantitative convergence for the solutions obtained from
solving the discretized twostage stochastic LCP (SLCP). We explain how the
discretized twostage SLCP may be solved by the wellknown progressive hedging
method (PHM). Moreover, we extend the discussion by considering a twostage
distributionally robust LCP (DRLCP) with moment constraints and proposing a
discretization scheme for the DRLCP. As an application, we show how the SLCP
and DRLCP models can be used to study equilibrium arising from twostage
duopoly game where each player plans to set up its optimal capacity at present
with anticipated competition for production in future.

An adaptive regularization algorithm using highorder models is proposed for
partially separable convexly constrained nonlinear optimization problems whose
objective function contains nonLipschitzian $\ell_q$norm regularization terms
for $q\in (0,1)$. It is shown that the algorithm using an $p$th order Taylor
model for $p$ odd needs in general at most $O(\epsilon^{(p+1)/p})$ evaluations
of the objective function and its derivatives (at points where they are
defined) to produce an $\epsilon$approximate firstorder critical point. This
result is obtained either with Taylor models at the price of requiring the
feasible set to be 'kernelcentered' (which includes bound constraints and many
other cases of interest), or for nonLipschitz models, at the price of passing
the difficulty to the computation of the step. Since this complexity bound is
identical in order to that already known for purely Lipschitzian minimization
subject to convex constraints [CartGoulToin2016], the new result shows that
introducing nonLipschitzian singularities in the objective function may not
affect the worstcase evaluation complexity order. The result also shows that
using the problem's partially separable structure (if present) does not affect
complexity order either. A final (worse) complexity bound is derived for the
case where Taylor models are used with a general convex feasible set.

Ultrasound (US) is the most commonly used liver imaging modality worldwide.
Due to its low cost, it is increasingly used in the followup of cancer
patients with metastases localized in the liver. In this contribution, we
present the results of an interactive segmentation approach for liver
metastases in US acquisitions. A (semi) automatic segmentation is still very
challenging because of the low image quality and the low contrast between the
metastasis and the surrounding liver tissue. Thus, the state of the art in
clinical practice is still manual measurement and outlining of the metastases
in the US images. We tackle the problem by providing an interactive
segmentation approach providing realtime feedback of the segmentation results.
The approach has been evaluated with typical US acquisitions from the clinical
routine, and the datasets consisted of pancreatic cancer metastases. Even for
difficult cases, satisfying segmentations results could be achieved because of
the interactive realtime behavior of the approach. In total, 40 clinical
images have been evaluated with our method by comparing the results against
manual ground truth segmentations. This evaluation yielded to an average Dice
Score of 85% and an average Hausdorff Distance of 13 pixels.

Virtual Reality, an immersive technology that replicates an environment via
computersimulated reality, gets a lot of attention in the entertainment
industry. However, VR has also great potential in other areas, like the medical
domain, Examples are intervention planning, training and simulation. This is
especially of use in medical operations, where an aesthetic outcome is
important, like for facial surgeries. Alas, importing medical data into Virtual
Reality devices is not necessarily trivial, in particular, when a direct
connection to a proprietary application is desired. Moreover, most researcher
do not build their medical applications from scratch, but rather leverage
platforms like MeVisLab, MITK, OsiriX or 3D Slicer. These platforms have in
common that they use libraries like ITK and VTK, and provide a convenient
graphical interface. However, ITK and VTK do not support Virtual Reality
directly. In this study, the usage of a Virtual Reality device for medical data
under the MeVisLab platform is presented. The OpenVR library is integrated into
the MeVisLab platform, allowing a direct and uncomplicated usage of the head
mounted display HTC Vive inside the MeVisLab platform. Medical data coming from
other MeVisLab modules can directly be connected per draganddrop to the
Virtual Reality module, rendering the data inside the HTC Vive for immersive
virtual reality inspection.

In the orthognathic surgery, dental splints are important and necessary to
help the surgeon reposition the maxilla or mandible. However, the traditional
methods of manual design of dental splints are difficult and timeconsuming.
The research on computeraided design software for dental splints is rarely
reported. Our purpose is to develop a novel special software named EasySplint
to design the dental splints conveniently and efficiently. The design can be
divided into two steps, which are the generation of initial splint base and the
Boolean operation between it and the maxillamandibular model. The initial
splint base is formed by ruled surfaces reconstructed using the manually picked
points. Then, a method to accomplish Boolean operation based on the distance
filed of two meshes is proposed. The interference elimination can be conducted
on the basis of marching cubes algorithm and Boolean operation. The accuracy of
the dental splint can be guaranteed since the original mesh is utilized to form
the result surface. Using EasySplint, the dental splints can be designed in
about 10 minutes and saved as a stereo lithography (STL) file for 3D printing
in clinical applications. Three phantom experiments were conducted and the
efficiency of our method was demonstrated.

In this publication, the interactive planning and reconstruction of cranial
3D Implants under the medical prototyping platform MeVisLab as alternative to
commercial planning software is introduced. In doing so, a MeVisLab prototype
consisting of a customized dataflow network and an own C++ module was set up.
As a result, the ComputerAided Design (CAD) software prototype guides a user
through the whole workflow to generate an implant. Therefore, the workflow
begins with loading and mirroring the patients head for an initial curvature of
the implant. Then, the user can perform an additional Laplacian smoothing,
followed by a Delaunay triangulation. The result is an aesthetic looking and
wellfitting 3D implant, which can be stored in a CAD file format, e.g.
STereoLithography (STL), for 3D printing. The 3D printed implant can finally be
used for an indepth presurgical evaluation or even as a real implant for the
patient. In a nutshell, our research and development shows that a customized
MeVisLab software prototype can be used as an alternative to complex commercial
planning software, which may also not be available in every clinic. Finally,
not to conform ourselves directly to available commercial software and look for
other options that might improve the workflow.

Let $A$ be a Koszul (or more generally, $N$Koszul) CalabiYau algebra.
Inspired by the works of Kontsevich, Ginzburg and Van den Bergh, we show that
there is a derived noncommutative Poisson structure on $A$, which induces a
graded Lie algebra structure on the cyclic homology of $A$; moreover, we show
that the Hochschild homology of $A$ is a Lie module over the cyclic homology
and the Connes long exact sequence is in fact a sequence of Lie modules.
Finally, we show that the LeibnizLoday bracket associated to the derived
noncommutative Poisson structure on $A$ is naturally mapped to the
Gerstenhaber bracket on the Hochschild cohomology of its Koszul dual algebra
and hence on that of $A$ itself. Relations with some other brackets in
literature are also discussed and several examples are given in detail.

In this paper, we prove that the global version of the ${\L}$ojasiewicz
gradient inequality holds for quadratic sphere constrained optimization problem
with exponent $\theta=\frac{3}{4}$. An example from Ting Kei Pong shows that
$\theta=\frac{3}{4}$ is tight. This is the first ${\L}$ojasiewicz gradient
inequality established for the sphere constrained optimization problem with a
linear term.

In this paper, we study a general optimization model, which covers a large
class of existing models for many applications in imaging sciences. To solve
the resulting possibly nonconvex, nonsmooth and nonLipschitz optimization
problem, we adapt the alternating direction method of multipliers (ADMM) with a
general dual stepsize to solve a reformulation that contains three blocks of
variables, and analyze its convergence. We show that for any dual stepsize
less than the golden ratio, there exists a computable threshold such that if
the penalty parameter is chosen above such a threshold and the sequence thus
generated by our ADMM is bounded, then the cluster point of the sequence gives
a stationary point of the nonconvex optimization problem. We achieve this via a
potential function specifically constructed for our ADMM. Moreover, we
establish the global convergence of the whole sequence if, in addition, this
special potential function is a Kurdyka{\L}ojasiewicz function. Furthermore,
we present a simple strategy for initializing the algorithm to guarantee
boundedness of the sequence. Finally, we perform numerical experiments
comparing our ADMM with the proximal alternating linearized minimization (PALM)
proposed in [5] on the background/foreground extraction problem with real data.
The numerical results show that our ADMM with a nontrivial dual stepsize is
efficient.

In this paper, we study the proximal gradient algorithm with extrapolation
for minimizing the sum of a Lipschitz differentiable function and a proper
closed convex function. Under the error bound condition used in [19] for
analyzing the convergence of the proximal gradient algorithm, we show that
there exists a threshold such that if the extrapolation coefficients are chosen
below this threshold, then the sequence generated converges $R$linearly to a
stationary point of the problem. Moreover, the corresponding sequence of
objective values is also $R$linearly convergent. In addition, the threshold
reduces to $1$ for convex problems and, as a consequence, we obtain the
$R$linear convergence of the sequence generated by FISTA with fixed restart.
Finally, we present some numerical experiments to illustrate our results.

Let $O_\tau(\Gamma)$ be a family of algebras \textit{quantizing} the
coordinate ring of $\mathbb{C}^2 / \Gamma$, where $\Gamma$ is a finite subgroup
of $\mathrm{SL}_2(\mathbb{C})$, and let $G_{\Gamma}$ be the automorphism group
of $O_\tau$. We study the natural action of $G_\Gamma$ on the space of right
ideals of $O_\tau$ (equivalently, finitely generated rank $1$ projective
$O_\tau$modules). It is known that the later can be identified with disjoint
union of algebraic (quiver) varieties, and this identification is
$G_\Gamma$equivariant. In the present paper, when $\Gamma \cong \mathbb{Z}_2$,
we show that the $G_{\Gamma}$action on each quiver variety is transitive. We
also show that the natural embedding of $G_\Gamma$ into $\mathrm{Pic}(O_\tau)$,
the Picard group of $O_\tau$, is an isomorphism. These results are used to
prove that there are countably many nonisomorphic algebras Morita equivalent
to $O_\tau$, and explicit presentation of these algebras are given. Since
algebras $O_\tau(\mathbb{Z}_2)$ are isomorphic to primitive factors of
$U(sl_2)$, we obtain a complete description of algebras Morita equivalent to
primitive factors. A structure of the group $G_{\Gamma}$, where $\Gamma$ is an
arbitrary cyclic group, is also investigated. Our results generalize earlier
results obtained for the (first) Weyl algebra $A_1$.

We consider a class of constrained optimization problems with a possibly
nonconvex nonLipschitz objective and a convex feasible set being the
intersection of a polyhedron and a possibly degenerate ellipsoid. Such problems
have a wide range of applications in data science, where the objective is used
for inducing sparsity in the solutions while the constraint set models the
noise tolerance and incorporates other prior information for data fitting. To
solve this class of constrained optimization problems, a common approach is the
penalty method. However, there is little theory on exact penalization for
problems with nonconvex and nonLipschitz objective functions. In this paper,
we study the existence of exact penalty parameters regarding local minimizers,
stationary points and $\epsilon$minimizers under suitable assumptions.
Moreover, we discuss a penalty method whose subproblems are solved via a
nonmonotone proximal gradient method with a suitable update scheme for the
penalty parameters, and prove the convergence of the algorithm to a KKT point
of the constrained problem. Preliminary numerical results demonstrate the
efficiency of the penalty method for finding sparse solutions of
underdetermined linear systems.

Ultrasound (US) is the most commonly used liver imaging modality worldwide.
It plays an important role in followup of cancer patients with liver
metastases. We present an interactive segmentation approach for liver tumors in
US acquisitions. Due to the low image quality and the low contrast between the
tumors and the surrounding tissue in US images, the segmentation is very
challenging. Thus, the clinical practice still relies on manual measurement and
outlining of the tumors in the US images. We target this problem by applying an
interactive segmentation algorithm to the US data, allowing the user to get
realtime feedback of the segmentation results. The algorithm has been
developed and tested handinhand by physicians and computer scientists to make
sure a future practical usage in a clinical setting is feasible. To cover
typical acquisitions from the clinical routine, the approach has been evaluated
with dozens of datasets where the tumors are hyperechoic (brighter), hypoechoic
(darker) or isoechoic (similar) in comparison to the surrounding liver tissue.
Due to the interactive realtime behavior of the approach, it was possible even
in difficult cases to find satisfying segmentations of the tumors within
seconds and without parameter settings, and the average tumor deviation was
only 1.4mm compared with manual measurements. However, the long term goal is to
ease the volumetric acquisition of liver tumors in order to evaluate for
treatment response. Additional aim is the registration of intraoperative US
images via the interactive segmentations to the patient's preinterventional CT
acquisitions.

The conjugate gradient (CG) method is an efficient iterative method for
solving largescale strongly convex quadratic programming (QP). In this paper
we propose some generalized CG (GCG) methods for solving the
$\ell_1$regularized (possibly not strongly) convex QP that terminate at an
optimal solution in a finite number of iterations. At each iteration, our
methods first identify a face of an orthant and then either perform an exact
line search along the direction of the negative projected minimumnorm
subgradient of the objective function or execute a CG subroutine that conducts
a sequence of CG iterations until a CG iterate crosses the boundary of this
face or an approximate minimizer of over this face or a subface is found. We
determine which type of step should be taken by comparing the magnitude of some
components of the minimumnorm subgradient of the objective function to that of
its rest components. Our analysis on finite convergence of these methods makes
use of an error bound result and some key properties of the aforementioned
exact line search and the CG subroutine. We also show that the proposed methods
are capable of finding an approximate solution of the problem by allowing some
inexactness on the execution of the CG subroutine. The overall arithmetic
operation cost of our GCG methods for finding an $\epsilon$optimal solution
depends on $\epsilon$ in $O(\log(1/\epsilon))$, which is superior to the
accelerated proximal gradient method [2,23] that depends on $\epsilon$ in
$O(1/\sqrt{\epsilon})$. In addition, our GCG methods can be extended
straightforwardly to solve boxconstrained convex QP with finite convergence.
Numerical results demonstrate that our methods are very favorable for solving
illconditioned problems.

This paper presents a generalized integrated framework of semiautomatic
surgical template design. Several algorithms were implemented including the
mesh segmentation, offset surface generation, collision detection, ruled
surface generation, etc., and a special software named TemDesigner was
developed. With a simple user interface, a customized template can be semi
automatically designed according to the preoperative plan. Firstly, mesh
segmentation with signed scalar of vertex is utilized to partition the inner
surface from the input surface mesh based on the indicated point loop. Then,
the offset surface of the inner surface is obtained through contouring the
distance field of the inner surface, and segmented to generate the outer
surface. Ruled surface is employed to connect inner and outer surfaces.
Finally, drilling tubes are generated according to the preoperative plan
through collision detection and merging. It has been applied to the template
design for various kinds of surgeries, including oral implantology, cervical
pedicle screw insertion, iliosacral screw insertion and osteotomy,
demonstrating the efficiency, functionality and generality of our method.

Percutaneous radiofrequency ablation (RFA) is a minimally invasive technique
that destroys cancer cells by heat. The heat results from focusing energy in
the radiofrequency spectrum through a needle. Amongst others, this can enable
the treatment of patients who are not eligible for an open surgery. However,
the possibility of recurrent liver cancer due to incomplete ablation of the
tumor makes postinterventional monitoring via regular followup scans
mandatory. These scans have to be carefully inspected for any conspicuousness.
Within this study, the RF ablation zones from twelve postinterventional CT
acquisitions have been segmented semiautomatically to support the visual
inspection. An interactive, graphbased contouring approach, which prefers
spherically shaped regions, has been applied. For the quantitative and
qualitative analysis of the algorithm's results, manual slicebyslice
segmentations produced by clinical experts have been used as the gold standard
(which have also been compared among each other). As evaluation metric for the
statistical validation, the Dice Similarity Coefficient (DSC) has been
calculated. The results show that the proposed tool provides lesion
segmentation with sufficient accuracy much faster than manual segmentation. The
visual feedback and interactivity make the proposed tool well suitable for the
clinical workflow.