\section{Introduction}

Retrieving 3D models from 2D sketches has important applications in computer graphics, information retrieval, and computer vision \cite{eitz2012sbsr, furuya2013ranking, li2014comparison}.
Compared to the early attempts where keywords or 3D shapes are used as queries \cite{shilane2004princeton}, the sketch-based idea is very attractive because sketches by hand provide an easy way to input, yet they are rich enough to specify shapes.

Directly matching 2D sketches to 3D models suffers from significant differences between the 2D and 3D representations.
Thus, in many state of the art methods 3D models are projected to multiple 2D views, and a sketch matches a 3D model if it matches one of its views.
Fig. \ref{fig:examples} shows a few examples of 2D sketches and their corresponding 3D models.
One can immediately see the variations in both the sketch styles and 3D models.

\begin{figure}[t]
\begin{center}
\begin{tabular}{c}
	\includegraphics[width=0.75\columnwidth]{images/sketch_model_sample} \\
\end{tabular}
\end{center}
\caption{Examples of sketch based 3D shape retrieval.\label{fig:examples}}
\end{figure}

In almost all state of the art approaches, sketch based 3D shape retrieval amounts to finding the ``best views'' for 3D models and hand-crafting the right features for matching sketches and views.
First, an automatic procedure is used to select the most representative views of a 3D model.
Ideally, one of the viewpoints is similar to that of the query sketches.
Then, 3D models are projected to 2D planes using a variety of line rendering algorithms.
Subsequently, many 2D matching methods can be used for computing the similarity scores, where features are always manually defined (\textit{e.g.}, Gabor, dense SIFT, and GALIF \cite{eitz2012sbsr}).

This stage-wise methodology appears pragmatic, but it also brings a number of puzzling issues.
To begin with, there is no guarantee that the best views have similar viewpoints with the sketches.
The inherent issue is that identifying the best views is an unsolved problem on its own, partially because the general definition of best views is elusive.
In fact, many best view methods require manually selected viewpoints for training, which makes the view selection by finding ``best views'' a chicken-egg problem.

Further, this viewpoint uncertainty makes it dubious to match samples from two different domains without learning their metrics.
Take Fig. \ref{fig:examples} for example, even when the viewpoints are similar the variations in sketches as well as the different characteristics between sketches and views are beyond the assumptions of many 2D matching methods.

Considering all the above issues arise when we struggle to seek the viewpoints for matching, can we bypass the stage of view selection?
In this paper we demonstrate that by learning cross domain similarities, we no longer require the seemingly indispensable view similarity assumption.

Instead of relying on the elusive concept of ``best views'' and hand-crafted features, we propose to define our views and learn features for views and sketches.
Assuming that the majority of the models are upright, we drastically reduce the number of views to two per object for the whole dataset.
We also make no selections of these two directions as long as they are significantly different.
Therefore, we consider this as the minimalism approach as opposed to multiple best views.

This upright assumption appears to be strong, but it turns out to be sensible for 3D datasets.
Many 3D models are naturally generated upright (\textit{e.g.}, \cite{shilane2004princeton}).
We choose two viewpoints because it is very unlikely to get degenerated views for two significantly different viewpoints.
An immediate advantage is that our matching is more efficient without the need of comparing to more views than necessary.

This seemingly radical approach triumphs only when the features are learned properly.
In principle, this can be regarded as learning representations between sketches and views by specifying similarities, which gives us a semantic level matching.
To achieve this, we need comprehensive shape representations rather than the combination of  shallow features that only capture low level visual information.

We learn the shape representations using Convolutional Neural Network (CNN).
Our model is based on the Siamese network \cite{chopra2005learning}.
Since the two input sources have distinctive intrinsic properties, we use two different CNN models, one for handling the sketches and the other for the views.
This two model strategy can give us more power to capture different properties in different domains.

Most importantly, we define a loss function to ``align'' the results of the two CNN models.
This loss function couples the two input sources into the same target space, which allows us to compare the features directly using a simple distance function.

Our experiments on three large datasets show that our method significantly outperforms state of the art approaches in a number of metrics, including precision-recall and the nearest neighbor.
We further demonstrate the retrievals in each domain are effective.
Since our network is based on filtering, the computation is fast.

Our contributions include
\begin{itemize}
\item We propose to learn feature representations for sketch based shape retrieval, which bypasses the dilemma of best view selection;
\item We adopt two Siamese Convolutional Neural Networks to successfully learn similarities in both the within-domain and the cross domain;
\item We outperform all the state of the art methods on three large datasets significantly.
\end{itemize}

\section{Related work}

Sketch based shape retrieval has received many interests for years \cite{Funkhouser:2003:ASE}.
In this section we review three key components in sketch based shape retrieval: public available datasets, features, and similarity learning.

\vspace{-5pt}
\paragraph*{Datasets}
The effort of building 3D datasets can be traced back to decades ago.
The Princeton Shape Benchmark (PSB) is probably one of the best known sources for 3D models \cite{shilane2004princeton}.
There are some recent advancements for general and special objects, such as the SHREC'14 Benchmark \cite{LLL+14a} and the Bonn Architecture Benchmark \cite{Wessel:2009:SBR:2381128.2381141}.

2D sketches have been adopted as input in many systems \cite{journals/ijcv/DarasA10}.
However, the large scale collections are available only recently.
Eitz \textit{et al.} \cite{eitz2012sbsr} collected sketches based on the PSB dataset. Li \textit{et al.} \cite{li2014comparison} organized the sketches collected by \cite{eitz2012hdhso} in their SBSR challenge.

\vspace{-5pt}
\paragraph*{Features}
Global shape descriptors, such as statistics of shapes \cite{Osada:2002:SD} and distance functions \cite{Kazhdan:2002:RSD:645316.649205}, have been used for 3D shape retrieval \cite{Tangelder:2008:SCB:1395016.1395041}.
Recently, local features is proposed for partial matching \cite{Funkhouser:2006:PMS:1281957.1281974} or used in the bag-of-words model for 3D shape retrieval \cite{Bronstein:2011:SGG:1899404.1899405}.

Boundary information together with internal structures are used for matching sketches against 2D projections.
Therefore, a good representation of line drawing images is a key component for sketch based shape retrieval.
Sketch representation such as shape context \cite{belongie2002shape} was proposed for image based shape retrieval.
Furuya \textit{et al.} proposed BF-DSIFT feature, which is an extended SIFT feature with Bag-of-word method, to represent sketch images \cite{furuya2009dense}.
One recent method is the Gabor local line based feature (GALIF) by Mathias \textit{et al.}, which builds on a bank of Gabor filters followed by a Bag-of-word method \cite{eitz2012sbsr}.

In addition to 2D shape features, some methods also explored geometry features as well as graph-based features to facilitate the 3D shape retrieval \cite{li2015comparison}.
Semantic labeling is also used to bridge the gaps between different domains \cite{gong2013learning}.
In this paper, we focus on view based method and only use 2D shape features.

\vspace{-5pt}
\paragraph*{CNN and Siamese network}
Recently deep learning has achieved great success on many computer vision tasks. Specifically, CNN has set records on standard object recognition benchmarks \cite{Krizhevsky2012}. With a deep structure, the CNN can effectively learn complicated mappings from raw images to the target, which requires less domain knowledge compared to handcrafted features and shallow learning frameworks.

A Siamese network \cite{chopra2005learning} is a particular neural network architecture consisting of two identical sub-convolutional networks, which is used in a weakly supervised metric learning setting.
The goal of the network is to make the output vectors similar if input pairs are labeled as similar, and dissimilar for the input pairs that are labeled as dissimilar.
Recently, the Siamese network has been applied to text classification \cite{Yih:2011:LDP:2018936.2018965} and speech feature classification \cite{NIPS2011_4314}.

\section{Learning feature representations for sketch based 3D shape retrieval}

We first briefly introduce basic concepts in CNNs and Siamese network.
Then, we present our network architecture for cross domain matching, based on the Siamese network.
Given a set of view and sketch pairs, we propose to use two different Siamese networks, one for each domain.
Finally, we revisit the view selection problem, and describe our minimalism approach of viewpoint definition and the line drawing rendering procedure.

\subsection{CNN and Siamese network}

CNN is a multilayer learning framework, which consists of an input layer, a few convolutional layers and fully connected layers, as well as an output layer on which the loss function is defined. The goal of CNN is to learn a hierarchy of feature representations. Signals in each layer are convolved with a number of filters and further downsampled by pooling operations, which aggregate values in a small region by functions including max, min, and average.
The learning of CNN is based on Stochastic Gradient Descent (SGD). Please refer to \cite{LeCun:1998:CNI:303568.303704} for details.

Siamese Convolutional Neural Network has been used successfully for dimension reduction in weakly supervised metric learning.
Instead of taking single sample as input, the network typically takes a pair of samples, and the loss functions are usually defined over pairs.
A typical loss function of a pair has the following form:
\begin{align}
L(s_1, s_2, y) &= (1-y) \alpha D_w^2 + y \beta e^{\gamma D_w},
\label{eq:loss}
\end{align}
where $s_1$ and $s_2$ are two samples, $y$ is the binary similarity label, $D_w = {\Vert f(s_1; w_1)-f(s_2; w_2) \Vert}_1$ is the distance.
Following \cite{chopra2005learning}, we set $\alpha = \frac{1}{C_p}$, $\beta = C_n$, and $\gamma = \frac{-2.77}{C_n}$, where $C_p=0.2$ and $C_n=10$ are two constants.

This can be regarded as a metric learning approach.
Unlike methods that assign binary similarity labels to pairs, the network aims at bring the output feature vectors closer for input pairs that are labeled as similar, or push the feature vectors away if the input pairs are labeled as dissimilar.

The Siamese network is frequently illustrated as two identical networks for two different samples.
In each SGD iteration, pairs of samples are processed using two identical networks, and the error computed by Eq. \ref{eq:loss} is then back-propagated and the gradients are computed individually base on the two sample sets.
The Siamese network is updated by the average of these two gradients.

\subsection{Cross-domain matching using Siamese network}

In this section, we propose a method to match samples from two domains without the heavy assumption of view similarity.
We first provide our motivation using an illustrated sample.
Then, we propose our extension of the basic Siamese network.
Specifically, we use two different networks to handle sources from different domains.

\subsubsection{An illustrated example}

The matching problem in sketch based shape retrieval can be seen as a metric learning paradigm.
In each domain, the samples are mapped to some feature vectors.
The cross domain matching is successful if the features from each domain are ``aligned'' correctly.
\begin{figure}[h]
\begin{center}
\begin{tabular}{c|c}
	\includegraphics[width=0.34\columnwidth]{images/illustration_1} &
	\includegraphics[width=0.34  \columnwidth]{images/illustration_2} \\
	(a) & (b)
\end{tabular}
\end{center}
\caption{An illustrated example, a) the shapes in the original domain may be mixed, and b) after cross-domain metric learning, similar shapes in both domains are grouped together.}
\label{fig:illustration}
\end{figure}

This idea is illustrated in Fig. \ref{fig:illustration}.
Blue denotes samples in the sketch domain, and the orange denotes the ones in the view domain.
Different shapes denote different classes.
Before learning, the feature points from two different domains are initially mixed together
(Fig. \ref{fig:illustration}a).
If we learn the correct mapping using pair similarities in each domain as well as their cross-domain relations jointly, the two point sets may be correctly aligned in the feature space (Fig. \ref{fig:illustration}b).
After this cross domain metric learning, matching can be performed in both the same domain (sketch-sketch and view-view) and cross domain (sketch-view).

Note that, there are no explicit requirements about viewpoint similarity in this perspective ({\it i.e.,} whether the matched pairs are from the same viewpoints is less important).
Instead, the focus is the metric between the two domains and the mapping within the same domain.

\subsubsection{Two networks, one loss}
\begin{figure*}[t]
\begin{center}
\includegraphics[width=0.75\textwidth]{images/siamese-model-2}
\end{center}

\caption{Dimension reduction using Siamese network.\label{fig:Siamese}}
\end{figure*}

The basic Siamese network is commonly used for samples from the same domain. In the cross domain setting, we propose to extend the basic version to two Siamese networks, one for the view domain and the other for the sketch domain.
Then, we define the within-domain loss and the cross domain loss.
This hypothesis is supported in the Sec. \ref{sec:exp}.

Assuming we have two inputs from each domain, \textit{i.e.}, $s_1$ and $s_2$ are two sketches and $v_1$ and $v_2$ are two views.
For simplicity, we assume $s_1$ and $v_1$ are from the same class and $s_2$ and $v_2$ are from the same class as well.
Therefore, one label $y$ is enough to specify their relationships.

As a result, our loss function is composed by three terms: the similarity of sketches, the similarity of views, and the cross domain similarity.
\begin{align}
\mathcal{L}(s_1,& s_2, v_1, v_2, y) \nonumber
\\
                         &= L(s_1,s_2,y) + L(v_1,v_2,y) + L(s_1,v_1,y),
\label{eq:finalloss}
\end{align}
where $L(\cdot, \cdot, \cdot)$ is defined by Eq. \ref{eq:loss}.
Please note that, while the category information available in the dataset can be exploited to improve the performance, we do not use the category labels in the above framework.

\subsection{Network architecture}

Fig. \ref{fig:Siamese} shows the architecture of our network for the inputs being views and sketches, respectively.

We use the same network design for both networks, but they are learned separately.
Our input patch size is $100 \times 100$ for both sources.
The structure of the single CNN has three convolutional layers, each with a max pooling, one fully connected layer to generate the features, and one output layer to compute the loss (Eq. \ref{eq:finalloss}).

The first convolutional layer followed by a $4\times4$ pooling generates 32 response maps, each of size $22\times22$.
The second layer and pooling outputs 64 maps of size $8\times8$.
The third layer layer has 256 response maps, each pooled to a size of $3\times3$.
The 2304 features generated by the final pooling operation are linearly transformed to $64 \times 1$ features in the last layer.
Rectified linear units are used in all layers.

\subsection{View definitions and line drawing rendering}

We present our procedure of generating viewpoints and rendering 3D models.
As opposed to multiple views, we find it sufficient to use two views to characterize a 3D model because the chance that both views are degenerated is little.
Following this observation, we impose the minimal assumptions on choosing views for the whole dataset:
\begin{enumerate}
  \item Most of the 3D models in the dataset are up-right;
  \item Two viewpoints are randomly generated for the whole dataset, provided that the difference in their angles is larger than $45$ degrees.
\end{enumerate}

Fig. \ref{fig:views} shows some of our views in the PSB dataset.
The first row shows that the upright assumption does not require strict alignments of 3D models, because some models may not have well defined orientation.
Further, while the models are upright, they can still has different rotations.
\begin{figure}[h]
\begin{center}
\begin{tabular}{c}
	\includegraphics[width=0.65\columnwidth]{images/view_comp}\\
\end{tabular}
\end{center}
\caption{3D models viewed from predefined viewpoints.}
\label{fig:views}
\end{figure}

We want to stress that our approach does not eliminate the possibility of selecting more (best) views as input, but the comparisons among view selection methods are beyond the scope of this paper.

Once the viewpoints are chosen, we render the 3D models and generate 2D line drawings.
Rendering line drawings that include strong abstraction and stylization effects is a very useful topic in computer graphics, computer vision, and psychology.
Outer edges and internal edges both play an important role in this rendering process.
Therefore, we use the following descriptors: 1) closed boundaries and 2) Suggestive Contours \cite{DeCarlo:2003:SCF} (Fig. \ref{fig:rendering}).
\begin{figure}[h]
\begin{center}
\begin{tabular}{ccc}
	\includegraphics[width=0.2\columnwidth]{images/m110_shade} &
	\includegraphics[width=0.2\columnwidth]{images/m110_sc} &
	\includegraphics[width=0.2\columnwidth]{images/m110_comb} \\
	(a) Shaded & (b) SC & (c) Final
\end{tabular}
\end{center}
\caption{Rendering 3D models.}

\label{fig:rendering}
\end{figure}

\section{Experiments}\label{sec:exp}

We present our experiments on three recent large datasets in this section.
In all experiments our method outperforms the state of the arts in a number of well recognized metrics.
In additional to the cross-domain retrieval, we also present our within-domain retrieval results, which have not been reported in any other comparison methods.
These experiments demonstrate that our Siamese network successfully learns the feature representations for both domains.
The data and the code is available at \url{http://users.cecs.anu.edu.au/~yili/cnnsbsr/}.

\subsection{Datasets}

\paragraph{PSB / SBSR dataset}
The Princeton Shape Benchmark (PSB) \cite{shilane2004princeton} is widely used for 3D shape retrieval system evaluation, which contains 1814 3D models and is equally divided into training set and testing set.

In \cite{eitz2012sbsr}, the Shape Based Shape Retrieval (SBSR) dataset is collected based on the PSB dataset.
The 1814 hand drawn sketches are collected using Amazon Mechanical Turk.
In the collection process, participants are asked to draw sketches given only the name of the categories without any visual clue from the 3D models.

\paragraph{SHREC'13 \& '14 dataset}
Although the PSB dataset is widely used in shape retrieval evaluation, there is a concern that the number of sketches for each class in the SBSR dataset is not enough. Some classes have only very few instances (27 of 90 training classes have no more than 5 instances), while some classes have dominating number of instances, \textit{e.g.}, the ``fighter jet" class and the ``human" class have as many as 50 instances.

To remove the possible bias when evaluating the retrieval algorithms, Li \textit{et al.} \cite{li2014comparison} reorganized the PSB/SBSR dataset, and proposed a SHREC'13 dataset where a subset of PSB with 1258 models is used and the sketches in each classes has 80 instances.
These sketch instances are split in two sets: 50 for training and 30 for testing.
Please note, the number of models in each class still varies.
For example, the largest class has 184 instances but there are 23 classes containing no more than 5 models

Recently, SHREC'14 is proposed to address some above concerns \cite{LLL+14a}, which greatly enlarges the number of 3D models to 8987, and the number of classes is doubled. The large variation of this dataset makes it much more challenging, and the overall performance of all reported methods are very low (\textit{e.g.}, the accuracy for the best algorithm is only 0.16 for the top 1 candidate).
This is probably due to the fact that the models are from various sources and are arbitrarily oriented.
While our performance is still superior (see Fig. \ref{fig:SHRECstat}b and Table. \ref{tab:shreccmp}), we choose to present our results using the SHREC'13 dataset.

\begin{figure*}[htbp]
\begin{center}
\begin{tabular}{c}
	\includegraphics[width=0.85\textwidth]{images/eval_results/shape-retrieval-psb-2.png}
\end{tabular}
\end{center}
\caption{Retrieval examples of PSB/SBSR dataset. Cyan denotes the correct retrievals.}
\label{fig:tub}
\end{figure*}

\vspace{-10pt}
\paragraph{Evaluation criteria}
In our experiment, we use the above datasets and measure the performance using the following criteria:
1) \textit{Precision-recall curve} is calculated for each query and linear interpolated, then the final curve is reported by averaging all precision values for fixed recall rates;
2) \textit{Average precision (mAP)} is the area under the precision-recall curve;
3) \textit{Nearest neighbor (NN)} is used to measure the top 1 retrieval accuracy;
4) \textit{E-Measure (E)} is the harmonic mean of the precision and recall for the top 32 retrieval results;
5) \textit{First/second tier (FT/ST)} and \textit{Discounted cumulated gain (DCG)} as defined in the PSB statistics.

\subsection{Experimental settings}

\paragraph{Stopping criteria}
All three of the datasets had been split into training and testing sets, but no validation set was specified. Therefore, we terminated our algorithm after $50$ epochs for PSB/SBSR and $20$ for SHREC'13 dataset (or until convergence).
Multiple runs were performed and the mean values were reported.

\vspace{-5pt}
\paragraph{Generating pairs for Siamese network}
To make sure we generate reasonable proportion of similar and dissimilar pairs, we use the following approach to generate pair sets.
For each training sketch, we random select $k_p$ view pairs in the same category (matched pairs) and $k_n$ view samples from other categories (unmatched pairs).
Usually, our dissimilar pairs are ten times more than the similar pairs for successful training.
In our experiment, we use $k_p=2$, $k_n=20$.
We perform this random pairing for each training epoch.
To increase the number of training samples, we also used data augmentation for the sketch set.
To be specific, we randomly perform affine transformations on each sketch sample with small scales and angles to generate more variations. We generate two augmentations for each sketch sample in the dataset.

\paragraph{Computational cost}
The implementation of the proposed Siamese CNN is based on the Theano \cite{Theano} library.
We measure the processing time on on a PC with 2.8GHz CPU and GTX 780 GPU.
With preprocessed view features, the retrieval time for each query is approximately 0.002 sec on average on SHREC'13 dataset.

The training time is proportional to the total number of pairs and the number of epochs.
Overall training  takes approximately 2.5 hours for PSB/SBSR, 6 hours for SHREC'13, respectively.
Considering the total number of pairs is large, the training time is sensible.

We test various number of views in our experiments. We find that there was no significant performance gain when we vary the view from two to ten. However, it increased the computational cost significantly when more views are used, and more importantly, the GPU memory. This motivates us to select only two views in the experiments.

\subsection{Shape retrieval on PSB/SBSR dataset}

\subsubsection{Examples}

In this section, we test our method using the PSB/SBSR dataset.
First, we show some retrieval examples in Fig. \ref{fig:tub}.
The first column shows 8 queries from different classes, and each row shows the top 15 retrieval results.
Cyan denotes the correct retrievals, and gray denotes incorrect ones.

Our method performs exceptionally well in popular classes such as human, face, and plane.
We also find that some fine grained categorizations are difficult to distinguish.
For instance, the shelf and the box differ only in a small part of the model.
However, we also want to note that some of the classes only differ in semantics (\textit{e.g.}, barn and house only differ in function).
Certainly, this semantic ambiguity is beyond the scope of this paper.

Finally, we want to stress that the importance of viewpoint is significantly decreased in our metric learning approach.
Some classes may exhibit a high degree of freedom such as the plane, but the retrieval results are also excellent (as shown in Fig. \ref{fig:tub}).

\subsubsection{Analysis}

We further show some statistics on this dataset.
First, we provide the precision-recall values at fixed points in Table \ref{tab:pr}.
Compared to Fig. 9 in \cite{eitz2012sbsr}, our results are approximately $10 \%$ higher.
We then show six standard evaluation metrics in Table \ref{tab:psbmetric}.
Since other methods did not report the results on this dataset, we leave the comprehensive comparison to the next section.
Instead, in this analysis we focus on the effectiveness of metric learning for shape retrieval.

PSB/SBSR is a very imbalanced dataset, where training and testing only partially overlap.
Namely, there are 21 classes appear in both training and testing sets, while 71 classes are used solely for testing.
This makes it an excellent dataset for investigating similarity learning, because the ``unseen'' classes verify the learning is not biased.

We show some examples for these unseen classes in Fig. \ref{fig:tub-test-only} (more statistical curves are available on project website due to the space limitation).
It is interesting to see that our proposed method works well even on failure cases ({\it e.g.,} the flower), where the retrieval returns similar shapes (``potting plant'').
This demonstrates that our method learns the similarity effectively.

\begin{table}[h]
\caption{Precision-recall on fixed points.}
\begin{center}
\begin{tabular}{c|c|c|c|c|c}

   \hline
    5$\%$ &  20$\%$ & 40$\%$ & 60$\%$ & 80$\%$ & 100$\%$   \\ \hline
   0.616 & 0.286 & 0.221 & 0.180 & 0.138 & 0.072   \\ \hline

\end{tabular}
\end{center}
\label{tab:pr}
\end{table}

\begin{table}[h]
\caption{Standard metrics on the PSB/SBSR dataset.}
\begin{center}
\begin{tabular}{c|c|c|c|c|c}

   \hline
    NN &  FT & ST & E & DCG  & mAP  \\ \hline
   0.223 & 0.177 & 0.271 & 0.173 & 0.451 & 0.218   \\ \hline

\end{tabular}
\end{center}
\label{tab:psbmetric}
\end{table}

\begin{figure}[htbp]
\begin{center}
\begin{tabular}{c}
	\includegraphics[width=0.9\columnwidth]{images/eval_results/shape-retrieval-test-only.png}
\end{tabular}
\end{center}
\caption{Retrieval examples of unseen samples in PSB/SBSR dataset. The cyan denotes the correct retrievals.}
\label{fig:tub-test-only}
\end{figure}

\subsection{Shape retrieval on SHREC'13 dataset}
In this section, we use the SHREC'13 benchmark to evaluate our method.
We also show the retrieval results within the same domain.

\vspace{-5pt}
\subsubsection{A visualization of the learned features}
First, we present a visualization of our learned features in Fig. \ref{fig:vis}.
We perform PCA on the learned features and reduce the dimension to two for visualization.
The green dots denote the sketches, and the yellow ones denote views.
For simplicity, we only overlay the views over the point cloud.
Please visit  \url{http://users.cecs.anu.edu.au/~yili/cnnsbsr/} for an interactive demo.

While this is a coarse visualization, we can already see some interesting properties of our method.
First, we can see that classes with similar shapes are grouped together automatically.
On the top right, different animals are mapped to neighboring positions.
On the left, various types of vehicles are grouped autonomously.
Other examples include house and church, which are very similar.
Note that this is an weakly supervised method. This localization suggests that the learned features are very useful for both within-domain and cross domain retrievals.
\begin{figure}[htbp]
\begin{center}
\begin{tabular}{c}
	\includegraphics[width=0.95\columnwidth]{images/view_pca6.png}
\end{tabular}
\end{center}
\caption{Visualization of feature space on SHREC'13. Sketch and view feature points are shown by green \& yellow, respectively.}
\label{fig:vis}
\end{figure}

\vspace{-5pt}
\subsubsection{Statistical results}
We present the statistical results on SHREC'13 in this section.
First, we compare the precision-recall curve against the state of the art methods reported in \cite{li2014comparison}.
\begin{figure*}[htbp]
\begin{center}
\begin{tabular}{cc}
	\includegraphics[height=5cm]{images/eval_results/shrec13-comp.eps}  &
	\includegraphics[height=5cm]{images/eval_results/shrec14-comp.eps} \\
	(a) SHREC'13 comparison  &  (b) SHREC'14 comparison\\
\end{tabular}
\end{center}
\caption{Performance comparison on SHREC'13 \& '14. Please refer to \cite{li2014comparison} and \cite{LLL+14a} for the descriptions of the compared methods.}
\label{fig:SHRECstat}
\end{figure*}

From the Fig. \ref{fig:SHRECstat} we can see that our method significantly outperforms other comparison methods.
On SHREC'13 benchmark, the performance gain of our method is already $10 \%$ when recall is small.
More importantly, the whole curve decreases much slower than other methods when the recall increases, which is desirable because it shows the method is more stable.
Our method has a higher performance gain ($30 \%$) when recall reaches $1$.

We note that there is a noticeable overfitting in the training when a stopping criterion is reached.
It suggests the performance can be even better, if one can fine tune and explore the network structure and training procedure.

We further show the standard metrics for comparison.
These metrics examine the retrieval from different perspectives.
For simplicity, we only select the best method from each research group in  \cite{li2014comparison}.
As shown in Table \ref{tab:shreccmp}, our method performs better in every metric on both benchmarks.
This further demonstrates our method is superior.

\begin{table}
\caption{Comparison on SHREC'13 \& '14 dataset. The best results are shown in red, and the second best results are shown in blue.}
\begin{center}
\begin{tabular}{c|c|c|c|c|c|c}

 \hline

\multicolumn{7}{c}{SHREC'13} \\ \hline

   &  NN    &   FT    &   ST    &   E   &    DCG  &  mAP \\ \hline \hline

Ours         &   \color{red}{0.405} & \color{red}{0.403} & \color{red}{0.548} & \color{red}{0.287} & \color{red}{0.607}  & \color{red}{0.469}  \\ \hline
Identic   &   \color{blue}{0.389} & \color{blue}{0.364} & \color{blue}{0.516} & \color{blue}{0.272} & \color{blue}{0.588}  & \color{blue}{0.434} \\ \hline \hline

\cite{furuya2013ranking}  & 0.279  & 0.203  & 0.296  & 0.166  & 0.458  & 0.250 \\ \hline
\cite{li2014comparison}  & 0.164  & 0.097  & 0.149  & 0.085  & 0.348  & 0.116 \\ \hline
\cite{journals/vlc/SousaF10}  & 0.017  & 0.016  & 0.031  & 0.018  & 0.240  & 0.026 \\ \hline
\cite{Saavedra:2012:SMR:2381194.2381204}   & 0.110  & 0.069  & 0.107  & 0.061  & 0.307  & 0.086 \\ \hline

\hline

\multicolumn{7}{c}{SHREC'14} \\ \hline
   &  NN    &   FT    &   ST    &   E   &    DCG  &  mAP \\ \hline \hline

Ours         &   \color{red}{0.239} & \color{red}{0.212} & \color{red}{0.316} & \color{red}{0.140} & \color{red}{0.496}  & \color{red}{0.228}  \\ \hline
\hline
\cite{tatsuma2012large}  & 0.160  & 0.115  & 0.170  & 0.079  & 0.376  & 0.131  \\ \hline
\cite{furuya2013ranking}  & 0.109  & 0.057  & 0.089  & 0.041  & 0.328  & 0.054  \\ \hline
\cite{li2014comparison}  & 0.095  & 0.050  & 0.081  & 0.037  & 0.319  & 0.050  \\ \hline

\end{tabular}
\end{center}
\label{tab:shreccmp}
\end{table}

We also compare to the case where both networks are identical, \textit{i.e.}, both views and sketches use the same Siamese network.
Fig. \ref{fig:SHRECstat}a suggests that this configuration is inferior than our proposed version, but still it is better than all other methods.
This supports our hypothesis that the variations in two domains are different.
This also sends a message that using the same features (hand-crafted or learned) for both domains may not be ideal.

\subsubsection{Within-domain retrieval}
Finally, we show the retrievals in the same domain.
This interesting experiment shall be straightforward to report because the data is readily available, but was not shown before in any literature.
Since this is a ``by-product'' of our method, we do not tune up any parameter or re-train the system.

Figs. \ref{fig:withindomain} and \ref{fig:withindomain2} visualize some retrieval results in each domain, respectively.
Table \ref{tab:withindomain} further reports the statistics.
The retrieval results demonstrate our method is powerful in learning the features for both within-domain and cross-domain.
From these figures, one can see that the view domain is much more consistent than the sketch domain.
Comparing Table \ref{tab:withindomain} to Table \ref{tab:shreccmp}, we conclude that the inconsistency in sketch is the most challenging issue in the sketch based 3D shape retrieval.

\begin{figure}[htbp]
\begin{center}
\begin{tabular}{c}
	\includegraphics[width=0.9\columnwidth]{images/eval_results/shape-retrieval-sketch-vs-sketch-top8.png} \\
\end{tabular}
\end{center}
\caption{Sketch-sketch retrieval for SHREC'13. The incorrect retrievals are marked as light gray.}
\label{fig:withindomain}
\end{figure}

\begin{figure}[htbp]
\begin{center}
\begin{tabular}{c}
	\includegraphics[width=0.9\columnwidth]{images/eval_results/shape-retrieval-view-vs-view-top8.png}
\end{tabular}
\end{center}
\caption{View-view retrieval for SHREC'13. The cyan denotes the correct retrievals.}
\label{fig:withindomain2}
\end{figure}

\begin{table}
\caption{Standard metrics for the within-domain retrieval on SHREC'13.}
\begin{center}
\begin{tabular}{c|c|c|c|c|c|c}

   \hline
             &   NN &  FT & ST & E & DCG  & mAP  \\ \hline
 view  &   0.965 & 0.877 & 0.982 & 0.536 & 0.971 & 0.909   \\ \hline
 sketch  & 0.431 & 0.352 & 0.514 & 0.298 & 0.679 & 0.373   \\ \hline

\end{tabular}
\end{center}
\label{tab:withindomain}
\end{table}

\section{Conclusion}
In this paper we propose to learn feature presentations for sketch based 3D shape retrieval. Instead of computing ``best views'' and match them against queries, we use predefined viewpoints for the whole dataset and adopt two Siamese CNNs, one for views and one for sketches.
Our experiments on three large datasets demonstrated that our method is superior.

\section{Acknowledgement}
NICTA is funded by the Australian Government
as represented by the Department of Broadband, Communications and the Digital Economy and the Australian Research Council through the ICT Centre of Excellence program.

{\small \bibliographystyle{ieee}


\section{Introduction}

We identify the $d$-dimensional flat torus $\mathbb{T}^{d}=\mathbb{R
^{d}/\mathbb{Z}^{d}$ with the unit cube $\left[  -\frac{1}{2},\frac{1
{2}\right)  ^{d}$ and we recall that a sequence $\{t_{j}\}_{j=1}^{+\infty
}\subset\mathbb{T}^{d}$ is \textit{uniformly distributed} if one of the
following three equivalent conditions is satisfied: (i) for every
$d$-dimensional box $I\subset\left[  -\frac{1}{2},\frac{1}{2}\right)  ^{d}$
with volume $|I|$
\[
\lim_{N\rightarrow+\infty}\frac{1}{N}\,\operatorname*{card}\{t_{j}\in I:1\leq
j\leq N\}=\left\vert I\right\vert \ ;
\]
(ii) for every continuous function $f$ on $\mathbb{T}^{d}$
\[
\lim_{N\rightarrow+\infty}\frac{1}{N}\sum_{j=1}^{N}f(t_{j})=\int
_{\mathbb{T}^{d}}f\left(  t\right)  \ dt\ ;
\]
and (iii) for every $0\neq k\in\mathbb{Z}^{d}$
\[
\lim_{N\rightarrow+\infty}\frac{1}{N}\sum_{j=1}^{N}e^{2\pi ik\cdot t_{j
}=0\ ,
\]
where \textquotedblleft$\cdot$\textquotedblright\ denotes the $d$-dimensional
inner product.

The concept of uniform distribution and the defining properties given above go
back to a fundamental paper written one hundred years ago by H. Weyl \cite{W};
see \cite{KN} for the basic reference on uniformly distributed sequences.
Observe that the above definition does not show the quality of a uniformly
distributed sequence. In the late thirties J. van der Corput coined the term
\textit{discrepancy}: let $\mathfrak{D}_{N}:=\{t_{j}\}_{j=1}^{N}$ be a
sequence of $N$ points in $\mathbb{T}^{d}$, henceforth called a
\textit{distribution} (of $N$ points), and let
\[
D\left(  \mathfrak{D}_{N}\right)  :=\sup_{I\subset\mathbb{T}^{d}}\left\vert
\operatorname*{card}\left(  \{t_{j}\}_{j=1}^{N}\cap I\right)
-N|I|\right\vert
\]
be the (non normalized) discrepancy associated with $\mathfrak{D}_{N}$ with
respect to the $d$-dimensional boxes $I$ in $\mathbb{T}^{d}$. There are
different approaches to define a discrepancy that measures the quality of a
distribution of points; see e.g. \cite{BC88,CST,DT,KN,Mat,DP} for an
introduction of discrepancy theory. See \cite{BCCGST,CF,DG,DHRZ} for the
connections of discrepancy to energy and numerical integration.

Throughout this paper we shall denote by $c,c_{1},\ldots$ positive constants
which may change from step to step.

\bigskip

K. Roth \cite{Roth} proved the following lower estimate: for every
distribution $\mathfrak{D}_{N}$ of $N$ points in $\mathbb{T}^{2}$, we have
\begin{equation}
\int_{\mathbb{T}^{2}}\left\vert \operatorname*{card}(\mathfrak{D}_{N}\cap
I_{x,y})-Nxy\right\vert ^{2}\ dxdy\geq c\ \log N\ , \label{roth
\end{equation}
where $I_{x,y}:=[0,x]\times\lbrack0,y]$ and $0\leq x,y<1$. This yields
$D(\mathfrak{D}_{N})\geq c\ \log^{1/2}N$. H. Davenport \cite{Dav} proved that
the estimate (\ref{roth}) is sharp.

W. Schmidt \cite{Scmdt} investigated the discrepancy with respect to discs.
His results were improved and extended, independently, by J. Beck \cite{Beck}
and H. Montgomery \cite{mont}: for every convex body $C\subset\left[
-\frac{1}{2},\frac{1}{2}\right)  ^{d}$ of diameter less than one and for every
distribution $\mathfrak{D}_{N}$ of $N$ points in $\mathbb{T}^{d}$, one ha
\begin{equation}
\int_{0}^{1}\int_{SO(d)}\int_{\mathbb{T}^{d}}\left\vert \operatorname*{card
(\mathfrak{D}_{N}\cap(\lambda\sigma(C)+t))-\lambda^{d}N|C|\right\vert
^{2}\ dtd\sigma d\lambda\geq c\ N^{(d-1)/d}\ . \label{bm
\end{equation}
This relation implies that for every distribution $\mathfrak{D}_{N}$ there
exists a translated, rotated, and dilated copy $\overline{C}$ of a given
convex body $C\subset\left[  -\frac{1}{2},\frac{1}{2}\right)  ^{d}$ having
diameter less than one, such that
\[
\left\vert \operatorname*{card}(\mathfrak{D}_{N}\cap\overline{C
)-N|\overline{C}|\right\vert \geq c\ N^{(d-1)/(2d)}\ .
\]
J. Beck and W. Chen \cite{BC} proved that (\ref{bm}) is sharp. Indeed, they
showed that for every positive integer $N$ there exists a distribution
$\widetilde{\mathfrak{D}}_{N}\subset\mathbb{T}^{d}$ satisfyin
\begin{equation}
\int_{SO(d)}\int_{\mathbb{T}^{d}}\left\vert \operatorname*{card
(\widetilde{\mathfrak{D}}_{N}\cap C)-N|C|\right\vert ^{2}\ dtd\sigma\leq
c\ N^{(d-1)/d}\ . \label{pallino discr
\end{equation}
This distribution $\widetilde{\mathfrak{D}}_{N}$ can be obtained either by
applying a probabilistic argument or by reduction to a lattice point problem;
see \cite{BGT, CHEN,CT,T} for a comparison of probabilistic and deterministic results.

\bigskip

In the following, we shall consider bounds for the integral in
(\ref{pallino discr}) for distributions of $N$ points that are restrictions of
a shrunk integer lattice to the unit cube $\left[  -\frac{1}{2},\frac{1
{2}\right)  ^{d}$. Due to an argument in \cite[p. 3533]{BIT} that also extends
to higher dimensions, we may assume that $N$ is a $d$th power $N=M^{d}$ for a
positive integer $M$. More precisely, we consider distribution

\[
\mathfrak{D}_{N}:=\left(  \frac{1}{N^{1/d}}\mathbb{Z}^{d}\right)  \cap\left[
-\frac{1}{2},\frac{1}{2}\right)  ^{d}\ .
\]
Given a convex body $C\subset\left[  -\frac{1}{2},\frac{1}{2}\right)  ^{d}$ of
diameter less than one, we then hav
\begin{equation}
\operatorname*{card}\left(  \mathfrak{D}_{N}\cap C\right)  -N\left\vert
C\right\vert =\operatorname*{card}\left(  \mathbb{Z}^{d}\cap N^{1/d}C\right)
-N\left\vert C\right\vert \ . \label{N ro lattice
\end{equation}
Estimation of the RHS in (\ref{N ro lattice}) is a classical lattice point
problem. Results concerning lattice points are extensively used in different
areas of pure and applied mathematics; see, for example, \cite{GL,H,K}.

For the definition of a suitable discrepancy function, we change the discrete
dilation $N^{1/d}$ in (\ref{N ro lattice}) to an arbitrary dilation $\rho
\geq1$ and replace the convex body $C$ in (\ref{N ro lattice}) with a
translated, rotated and then dilated copy $\rho\sigma\left(  C\right)  +t$,
where $\sigma\in SO\left(  d\right)  $ and $t\in\mathbb{T}^{d}$. Thus the
discrepanc
\[
D_{C}^{\rho}(\sigma,t):=\operatorname*{card}(\mathbb{Z}^{d}\cap(\rho
\sigma(C)+t))-\rho^{d}|C|=\sum_{k\in\mathbb{Z}^{d}}\chi_{\rho\sigma
(C)+t}(k)-\rho^{d}|C|
\]
is defined as the difference between the number of integer lattice points in
the set $\rho\sigma\left(  C\right)  +t$ and its volume $\rho^{d}\left\vert
C\right\vert $ (here, $\chi_{A}$ denotes the characteristic function for the
set $A$). It is easy to see (e.g., \cite{BGT}) that the periodic function
$t\mapsto D_{C}^{\rho}(\sigma,t)$ has the Fourier series expansio
\begin{equation}
\rho^{d}\sum_{0\neq m\in\mathbb{Z}^{d}}\widehat{\chi_{\sigma\left(  C\right)
}}\left(  \rho m\right)  \ e^{2\pi im\cdot t}\ . \label{fourier D
\end{equation}
D. Kendall \cite{Kendall} seems to have been the first to realize that
multiple Fourier series expansions can be helpful in certain lattice point
problems. Using our notation, he proved that for every convex body
$C\subset\mathbb{R}^{d}$ and $\rho\geq1
\begin{equation}
\Vert D_{C}^{\rho}\Vert_{L^{2}\left(  SO(d)\times\mathbb{T}^{d}\right)  }\leq
c\ \rho^{(d-1)/2}\ . \label{kend
\end{equation}
This also follows from more recent results in \cite{Pod91} and \cite{BHI} as
demonstrated next. Given a convex body $C\subset\mathbb{R}^{d}$, we define the
(spherical) average decay of $\widehat{\chi_{C}}$ a
\[
\left\Vert \widehat{\chi_{C}}\left(  \rho\cdot\right)  \right\Vert
_{L^{2}\left(  \Sigma_{d-1}\right)  }:=\left\{  \int_{\Sigma_{d-1}}\left\vert
\widehat{\chi_{C}}\left(  \rho\tau\right)  \right\vert ^{2}\ d\tau\right\}
^{1/2}\ ,
\]
where $\Sigma_{d-1}:=\left\{  t\in\mathbb{R}^{d}:\left\vert t\right\vert
=1\right\}  $ and $\tau$ is the rotation invariant normalized measure on
$\Sigma_{d-1}$. Extending an earlier result of A. Podkorytov \cite{Pod91}, L.
Brandolini, S. Hofmann, and A. Iosevich \cite{BHI} proved that
\begin{equation}
\left\Vert \widehat{\chi_{C}}\left(  \rho\cdot\right)  \right\Vert
_{L^{2}\left(  \Sigma_{d-1}\right)  }\leq c\ \rho^{-\left(  d+1\right)  /2}\ .
\label{bhi
\end{equation}
By applying the Parseval identity to the Fourier series (\ref{fourier D}) of
the discrepancy function, we obtain Kendall's result (\ref{kend}); i.e.
\begin{align}
\Vert D_{C}^{\rho}\Vert_{L^{2}\left(  SO(d)\times\mathbb{T}^{d}\right)  }^{2}
&  =\rho^{2d}\sum_{0\neq k\in\mathbb{Z}^{d}}\int_{SO\left(  d\right)
}\left\vert \widehat{\chi_{\sigma\left(  C\right)  }}\left(  \rho k\right)
\right\vert ^{2}\ d\sigma\label{parsev}\\
&  \leq c\ \rho^{2d}\sum_{0\neq k\in\mathbb{Z}^{d}}\left\vert \rho
k\right\vert ^{-\left(  d+1\right)  }\leq c_{1}\ \rho^{d-1}\ .\nonumber
\end{align}
We are interested in the reversed inequalit
\begin{equation}
\Vert D_{C}^{\rho}\Vert_{L^{2}\left(  SO(d)\times\mathbb{T}^{d}\right)
^{2}\geq c_{1}\ \rho^{d-1}\ , \label{reverse
\end{equation}
which, as we shall see, may or may not hold. To understand this, let us assume
that (\ref{bhi})\ can be reverse
\begin{equation}
\left\Vert \widehat{\chi_{C}}\left(  \rho\cdot\right)  \right\Vert
_{L^{2}\left(  \Sigma_{d-1}\right)  }\geq c_{1}\ \rho^{-\left(  d+1\right)
/2}\ . \label{simplex reverse
\end{equation}
This relation (\ref{simplex reverse}) is true for a simplex (see \cite[Theorem
2.3]{BCT}) but it is not true for every convex body (see the next section).

The following result was proved in \cite[Proof of Theorem 3.7]{BCT}.

\begin{proposition}
\label{prop}Let $C$ in $\mathbb{R}^{d}$ be a convex body which satisfies
(\ref{simplex reverse}). Then $C$ satisfies (\ref{reverse}).
\end{proposition}

\begin{proof}
Indeed
\begin{align}
\Vert D_{C}^{\rho}\Vert_{L^{2}\left(  SO(d)\times\mathbb{T}^{d}\right)  }^{2}
&  =\rho^{2d}\sum_{0\neq k\in\mathbb{Z}^{d}}\int_{SO\left(  d\right)
}\left\vert \widehat{\chi_{\sigma\left(  C\right)  }}\left(  \rho k\right)
\right\vert ^{2}\ d\sigma\label{dalbassocondecad}\\
&  \geq c\ \rho^{2d}\int_{SO\left(  d\right)  }\left\vert \widehat
{\chi_{\sigma\left(  C\right)  }}\left(  \rho k^{\prime}\right)  \right\vert
^{2}\ d\sigma\geq c_{1}\ \rho^{d-1}\ ,\nonumber
\end{align}
where $k^{\prime}$ is any non-zero element in $\mathbb{Z}^{d}$.
\end{proof}

We are going to see that (\ref{reverse}) does not imply (\ref{simplex reverse}).

\section{$L^{2}$-regularity of convex bodies}

We say that a convex body $C\subset\mathbb{R}^{d}$ is $L^{2}$\emph{-regular
}if there exists a positive constant $c_{1}$ such that
\begin{equation}
c_{1}\ \rho^{(d-1)/2}\leq\Vert D_{C}^{\rho}\Vert_{L^{2}\left(  SO(d)\times
\mathbb{T}^{d}\right)  } \label{L2reg
\end{equation}
(by (\ref{kend}) we already know that $\Vert D_{C}^{\rho}\Vert_{L^{2}\left(
SO(d)\times\mathbb{T}^{d}\right)  }\leq c_{2}$ $\rho^{(d-1)/2}$ for some
$c_{2}>0$). If (\ref{L2reg}) fails we say that $C$ is $L^{2}$\emph{-irregular}.

Let $d>1$. L. Parnovski and A. Sobolev \cite{PS} proved that the
$d$-dimensional ball $B_{d}:=\left\{  t\in\mathbb{R}^{d}:\left\vert
t\right\vert \leq1\right\}  $ is $L^{2}$-regular if and only if $d\not \equiv
1\ (\operatorname{mod}4)$.

More generally, it was proved \cite{BCGT} that if $C\subset\mathbb{R}^{d}$
($d>1$) is a convex body with smooth boundary, having everywhere positive
Gaussian curvature, then (i) if $C$ is not symmetric about a point, or if
$d\not \equiv 1\left(  \operatorname{mod}4\right)  $, then $C$ is $L^{2
$-regular; (ii) if $C$ is symmetric about a point and if $d\equiv1\,\left(
\operatorname{mod}4\right)  $ then $C$ is $L^{2}$-irregular.

L. Parnovski and N. Sidorova \cite{psi} studied the above problem for the
non-convex case of a $d$-dimensional annulus ($d>1$). They provided a complete
answer in terms of the width of the annulus.

In the case of a polyhedron $P$, inequality (\ref{kend}) was extended to
$L^{p}$ norms in \cite{BCT}: for any $p>1$ and $\rho\geq1$ we have
\[
\Vert D_{P}^{\rho}\Vert_{L^{p}\left(  SO(d)\times\mathbb{T}^{d}\right)  }\leq
c_{p}\ \rho^{(d-1)(1-1/p)
\]
and, specifically for simplices $S$, one ha
\[
c_{p}^{\prime}\ \rho^{(d-1)(1-1/p)}\leq\Vert D_{S}^{\rho}\Vert_{L^{p}\left(
SO(d)\times\mathbb{T}^{d}\right)  }\leq c_{p}\ \rho^{(d-1)(1-1/p)}\ .
\]
In particular, this implies that the $d$-dimensional simplices are $L^{2}$-regular.

For the planar case it was proved in \cite[Theorem 6.2]{BRT} that every convex
body with piecewise $C^{\infty}$ boundary that is not a polygon is $L^{2}$-regular.

Related results can be found in \cite{BCT,CT,ko}.

Until now no example of a $L^{2}$-irregular polyhedron has been found.

We are interested in identifying the $L^{2}$-regular convex polyhedrons. In
this paper we give a complete answer for the planar case.

\bigskip

Let us first compare the $L^{2}$-regularity for a disc $B\subset\mathbb{R
^{2}$ and a square $Q\subset\mathbb{R}^{2}$. Their characteristic functions
$\chi_{B}$ and $\chi_{Q}$ do not satisfy (\ref{simplex reverse}). Indeed,
$\widehat{\chi_{B}}\left(  \xi\right)  =\left\vert \xi\right\vert ^{-1
J_{1}\left(  2\pi\left\vert \xi\right\vert \right)  $, where $J_{1}$ is the
Bessel function (see e.g. \cite{T}). Then the zeroes of $J_{1}$ yield an
increasing diverging sequence $\left\{  \rho_{u}\right\}  _{u=1}^{\infty}$
such that
\[
\left\Vert \widehat{\chi_{C}}\left(  \rho_{u}\cdot\right)  \right\Vert
_{L^{2}\left(  \Sigma_{1}\right)  }=0\ .
\]

Less obvious is the fact that the inequality $\left\Vert \widehat{\chi_{C
}\left(  \rho\cdot\right)  \right\Vert _{L^{2}\left(  \Sigma_{1}\right)  }\geq
c\ \rho^{-3/2}$ fails for a square $Q$: it was observed in \cite{BCT} the
existence of a positive constant $c$ such that, for every positive integer
$n$, one ha
\begin{equation}
\left\Vert \widehat{\chi_{Q}}\left(  n\cdot\right)  \right\Vert _{L^{2}\left(
\Sigma_{1}\right)  }\leq c\ n^{-7/4}\ . \label{7quarti
\end{equation}
For completeness we write the short proof of (\ref{7quarti}). Indeed, let
$\ Q=\left[  -\frac{1}{2},\frac{1}{2}\right]  ^{2}$ and let $n$ be a positive
integer. Let $\Theta:=\left(  \cos\theta,\sin\theta\right)  $. Then an
explicit computation of $\widehat{\chi_{Q}}$ yield
\begin{align*}
&  \int_{0}^{2\pi}\left\vert \widehat{\chi_{Q}}(n\Theta)\right\vert
^{2}~d\theta=8\int_{0}^{\pi/4}\left\vert \dfrac{\sin(\pi n\cos\theta)}{\pi
n\cos\theta}\frac{\sin(\pi n\sin\theta)}{\pi n\sin\theta}\right\vert
^{2}~d\theta\\
&  \leq c\,\frac{1}{n^{4}}\int_{0}^{\pi/4}\left\vert \dfrac{\sin(\pi
n\cos\theta)}{\sin\theta}\right\vert ^{2}~d\theta=c\,\frac{1}{n^{4}}\int
_{0}^{\pi/4}\left\vert \dfrac{\sin(\pi n\left(  1-2\sin^{2}\left(
\theta/2\right)  \right)  )}{\sin\theta}\right\vert ^{2}~d\theta\\
&  \leq c^{\prime}\,\frac{1}{n^{4}}\int_{0}^{\pi/4}\left\vert \sin(2\pi
n\sin^{2}\left(  \theta/2\right)  )\right\vert ^{2}\theta^{-2}~d\theta\\
&  \leq c^{\prime\prime}\,\frac{1}{n^{4}}\int_{0}^{n^{-1/2}}n^{2}\theta
^{2}~d\theta+c^{\prime\prime}\ \frac{1}{n^{4}}\int_{n^{-1/2}}^{\pi/4
\theta^{-2}~d\theta\leq c^{\prime\prime\prime}\ n^{-7/2}\;.
\end{align*}

Then $B$ and $Q$ \textit{may} be $L^{2}$-irregular.

On the one hand it is known that a disc $B$ is $L^{2}$-regular (see \cite{PS}
or \cite[Theorem 6.2]{BRT}), so that (\ref{reverse}) does not imply
(\ref{simplex reverse}). On the other hand we shall prove in this paper that
$Q$ is $L^{2}$-irregular.

The $L^{2}$-irregularity of the square $Q$ is shared by each member of the
family of polygons described in the following definition.

\begin{definition}
\label{def Gotic B}Let $\mathfrak{P}$ be the family of all convex polygons in
$\mathbb{R}^{2}$ which can be inscribed in a circle and are symmetric about
the centre.
\end{definition}

\section{Statements of the results}

We now state our main result.

\begin{theorem}
\label{maintheorem}A convex polygon $P$ is $L^{2}$-regular if and only if
$P\notin\mathfrak{P}$.
\end{theorem}

The \textquotedblleft only if\textquotedblright\ part is a consequence of the
following more precise result.

\begin{proposition}
\label{propeps}If $P\in\mathfrak{P}$, then for every $\varepsilon>0$ there is
an increasing diverging sequence $\left\{  \rho_{u}\right\}  _{u=1}^{\infty}$
such tha
\[
\Vert D_{P}^{\rho_{u}}\Vert_{L^{2}\left(  SO(2)\times\mathbb{T}^{2}\right)
}\leq c_{\varepsilon}\ \rho_{u}^{1/2}\log^{-1/\left(  32+\varepsilon\right)
}(\rho_{u})\ .
\]

\end{proposition}

Theorem \ref{maintheorem} above and \cite[Theorem 6.2]{BRT} yield the
following more general result.

\begin{corollary}
Let $C$ be a convex body in $\mathbb{R}^{2}$ having piecewise smooth boundary.
Then $C$ is not $L^{2}$-regular if and only if it belongs to $\mathfrak{P}$.
\end{corollary}

The following result shows that Theorem \ref{maintheorem} is essentially sharp.

\begin{proposition}
\label{stimaepsilon} For every $P\in\mathfrak{P}$, for $\varepsilon>0$
arbitrary small, and for any $\rho$ large enough
\[
\Vert D_{P}^{\rho}\Vert_{L^{2}(SO(2)\times\mathbb{T}^{2})}\geq c_{\varepsilon
}\ \rho^{1/2-\varepsilon}\ ,
\]
where $c_{\varepsilon}$ is independent of $\rho$.
\end{proposition}

The \textquotedblleft if\textquotedblright\ part of Theorem \ref{maintheorem}
is a consequences of the following three lemmas.

\begin{lemma}
\label{lemma1}Let $P$ in $\mathbb{R}^{2}$ be a polygon having a side not
parallel to any other side. Then $P$ is $L^{2}$-regular.
\end{lemma}

\begin{lemma}
\label{lemma2}Let $P$ in $\mathbb{R}^{2}$ be a convex polygon with a pair of
parallel sides having different lengths. Then $P$ is $L^{2}$-regular.
\end{lemma}

\begin{lemma}
\label{lemma3} Let $P$ in $\mathbb{R}^{2}$ be a convex polygon which cannot be
inscribed in a circle. Then $P$ is $L^{2}$-regular.
\end{lemma}

\section{ Notation and preliminary arguments}

In the remainder of the paper, a polygon $P$ is given by its vertex set
$\{P_{h}\}_{h=1}^{s}$, where it is assumed that the numbering indicates
counterclockwise ordering of the vertices;\ we write $P\sim\{P_{h}\}_{h=1
^{s}$. For convenience we use periodic labeling; i.e., $P_{h+s},$ $P_{h+2s},$
$\ldots$ refer to the same point $P_{h}$ for $1\leq h\leq s$. For every $h$ le

\[
\tau_{h}:=\frac{P_{h+1}-P_{h}}{\left\vert P_{h+1}-P_{h}\right\vert
\]
be the direction of the oriented side $P_{h}P_{h+1}$ and $\ell_{h}:=\left\vert
P_{h+1}-P_{h}\right\vert $ its length.

For every $h$ let $\nu_{h}$ be the outward unit normal vector corresponding to
the side $P_{h}P_{h+1}$. Let
\[
\mathcal{L}_{h}:=\left\vert P_{h}+P_{h+1}\right\vert
\]
be the length of the vector $P_{h}+P_{h+1}$. Observe that if $\left\vert
P_{h}\right\vert =\left\vert P_{h+1}\right\vert $ (in particular if the
polygon $P$ is inscribed in a circle centred at the origin) then
\[
P_{h}+P_{h+1}=\mathcal{L}_{h}\nu_{h}\ .
\]
We shall always assume $\ell_{h}\geq1$ and $\mathcal{L}_{h}\geq1$.

Let $\nu(s)$ be the outward unit normal vector at a point $s\in\partial P$
which is not a vertex of $P$. By applying Green's formula we see that, for any
$\rho\geq1$, we hav
\begin{align}
&  \widehat{\chi_{P}}(\rho\Theta)\nonumber\\
&  =\int_{P}e^{-2\pi i\rho\Theta\cdot t}\ dt=-\frac{1}{2\pi i\rho
\int_{\partial P}e^{-2\pi i\rho\Theta\cdot s}\ \left(  \Theta\cdot
\nu(s)\right)  \ ds\nonumber\\
&  =-\frac{1}{2\pi i\rho}\sum_{h=1}^{s}\ell_{h}\left(  \Theta\cdot\nu
_{h}\right)  \int_{0}^{1}e^{-2\pi i\rho\Theta\cdot\left(  P_{h}+\lambda\left(
P_{h+1}-P_{h}\right)  \right)  }\ d\lambda\nonumber\\
&  =-\frac{1}{4\pi^{2}\rho^{2}}\sum_{h=1}^{s}\frac{\Theta\cdot\nu_{h}
{\Theta\cdot\tau_{h}}\ \left[  e^{-2\pi i\rho\Theta\cdot P_{h+1}}-e^{-2\pi
i\rho\Theta\cdot P_{h}}\right] \nonumber\\
&  =-\frac{1}{4\pi^{2}\rho^{2}}\sum_{h=1}^{s}\frac{\Theta\cdot\nu_{h}
{\Theta\cdot\tau_{h}}\ e^{-\pi i\rho\Theta\cdot(P_{h+1}+P_{h})}\left[  e^{-\pi
i\rho\Theta\cdot(P_{h+1}-P_{h})}-e^{\pi i\rho\Theta\cdot(P_{h+1}-P_{h
)}\right] \nonumber\\
&  =\frac{i}{2\pi^{2}\rho^{2}}\sum_{h=1}^{s}\frac{\Theta\cdot\nu_{h}
{\Theta\cdot\tau_{h}}\ e^{-\pi i\rho\mathcal{L}_{h}\Theta\cdot\nu_{h}}\sin
(\pi\rho\ell_{h}\Theta\cdot\tau_{h})\ . \label{divergenza
\end{align}

For any $1\leq h\leq s$, let $\theta_{h}\in\lbrack0,2\pi)$ be the angle
defined by
\begin{equation}
\tau_{h}=:(\cos\theta_{h},\sin\theta_{h})\ . \label{thetaj
\end{equation}
Hence
\begin{equation}
\nu_{h}=(\sin\theta_{h},-\cos\theta_{h}) \label{nu j
\end{equation}
and, if $\Theta:=\left(  \cos\theta,\sin\theta\right)  $,
\[
\Theta\cdot\tau_{h}=\cos(\theta-\theta_{h})\ ,\ \ \ \ \ \Theta\cdot\nu
_{h}=-\sin(\theta-\theta_{h}).
\]
Then (\ref{divergenza}) can be written a
\[
\widehat{\chi_{P}}(\rho\Theta)=-\frac{i}{2\pi^{2}\rho^{2}}\sum_{h=1}^{s
\frac{\sin\left(  \theta-\theta_{h}\right)  }{\cos\left(  \theta-\theta
_{h}\right)  }\ e^{\pi i\rho\mathcal{L}_{h}\sin\left(  \theta-\theta
_{h}\right)  }\sin\left(  \pi\rho\ell_{h}\cos\left(  \theta-\theta_{h}\right)
\right)
\]
and the equality in (\ref{parsev}) yields
\begin{align}
&  \Vert D_{P}^{\rho}\Vert_{L^{2}\left(  SO(2)\times\mathbb{T}^{2}\right)
}^{2}\label{generale}\\
&  =\rho^{4}\sum_{0\neq k\in\mathbb{Z}^{2}}\int_{0}^{2\pi}\left\vert
\widehat{\chi_{P}}(\rho\left\vert k\right\vert \Theta)\right\vert
^{2}\ d\theta\nonumber\\
&  =c\sum_{0\neq k\in\mathbb{Z}^{2}}\frac{1}{|k|^{4}}\nonumber\\
&  \times\int_{0}^{2\pi}\left\vert \sum_{h=1}^{s}\frac{\sin(\theta-\theta
_{h})}{\cos(\theta-\theta_{h})}\ e^{-\pi i\rho\left\vert k\right\vert
\mathcal{L}_{h}\sin\left(  \theta-\theta_{h}\right)  }\sin(\pi\rho\left\vert
k\right\vert \ell_{h}\cos(\theta-\theta_{h}))\right\vert ^{2}d\theta
\ .\nonumber
\end{align}

For $P\in\mathfrak{P}$, relation (\ref{generale}) can be further simplified.
Let $P\in\mathfrak{P}$ have $s=2n$ sides (i.e. $P\sim\{P_{h}\}_{h=1}^{2n}$)
and be inscribed in a circle centered at the origin. Then $P_{h
P_{h+1}=-P_{n+h}P_{n+h+1}$ for any $1\leq h\leq n$ and $P_{h+1}+P_{h
=\mathcal{L}_{h}\nu_{h}$. Therefore, for every $1\leq h\leq n$,
\[
\tau_{h}=-\tau_{n+h}\ ,\ \ \ \ \ \nu_{h}=-\nu_{n+h}\ ,\ \ \ \ \ \ell_{h
=\ell_{n+h}\ ,\ \ \ \ \ \mathcal{L}_{h}=\mathcal{L}_{n+h}\ .
\]
Then the relation (\ref{divergenza}) become
\begin{equation}
\widehat{\chi}_{P}(\rho\Theta)=\frac{1}{\pi^{2}\rho^{2}}\sum_{h=1}^{n
\frac{\sin(\theta-\theta_{h})}{\cos(\theta-\theta_{h})}\sin(\pi\rho
\mathcal{L}_{h}\sin(\theta-\theta_{h}))\sin(\pi\rho\ell_{h}\cos(\theta
-\theta_{h}))\ . \label{divergenza2
\end{equation}
and the equality in (\ref{parsev}) yields
\begin{align}
&  \Vert D_{P}^{\rho}\Vert_{L^{2}\left(  SO(2)\times\mathbb{T}^{2}\right)
}^{2}\label{final}\\
&  =c\sum_{0\neq k\in\mathbb{Z}^{2}}\frac{1}{|k|^{4}}\nonumber\\
&  \times\int_{0}^{2\pi}\left\vert \sum_{h=1}^{n}\frac{\sin(\theta-\theta
_{h})}{\cos(\theta-\theta_{h})}\sin(\pi\rho|k|\ell_{h}\cos(\theta-\theta
_{h}))\sin(\pi\rho|k|\mathcal{L}_{h}\sin(\theta-\theta_{h}))\right\vert
^{2}d\theta\nonumber\\
&  \leq c\ \sum_{0\neq k\in\mathbb{Z}^{2}}\frac{1}{|k|^{4}}\sum_{h=1}^{n
\int_{0}^{\pi/2}\left\vert \frac{\sin(\pi\rho|k|\ell_{h}\sin\theta)
{\sin\theta}\sin(\pi\rho|k|\mathcal{L}_{h}\cos\theta)\right\vert ^{2
d\theta\ .\nonumber
\end{align}
The last relation holds for every $P\in\mathfrak{P}$ with $2n$ sides.

\section{Proofs}

\begin{proof}
[Proof of Lemma \ref{lemma1}]The proof of Lemma \ref{lemma1} is essentially
the proof of \cite[Theorem 3.7]{BCT}, which is stated for a simplex but the
argument also works for every polyhedron having a face not parallel to any
other face.
\end{proof}

\begin{proof}
[Proof of Lemma \ref{lemma2}]By Lemma \ref{lemma1} we can assume that
$P\sim\{P_{h}\}_{h=1}^{2n}$ is a convex polygon with an even number of sides,
and that for every $h=1,\ldots,n$ the sides $P_{h}P_{h+1}$ and $P_{h+n
P_{h+n+1}$ are parallel. Suppose that the length $\ell_{j}$ of the $j$th side
$P_{j}P_{j+1}$ is longer than the length $\ell_{j+n}$ of the opposite side
$P_{j+n}P_{j+n+1}$. Then there exist $0<\varepsilon<1$ and $0<\alpha<1$ such
tha
\begin{equation}
(1+\varepsilon)\frac{\ell_{j+n}}{\ell_{j}}<\alpha\ . \label{alfa
\end{equation}
Let $H>1$ be a large constant satisfying
\begin{equation}
\sin\left(  \theta-\theta_{j}\right)  \geq\sqrt{\alpha}\left(  \theta
-\theta_{j}\right)  \ \ \ \ \ \text{if \ }0\leq\theta-\theta_{j}\leq
\frac{1+\varepsilon}{H}\ . \label{alfa A
\end{equation}
We further assume (recall $\rho\geq1$
\[
\frac{1}{H\pi\rho\ell_{j}}\leq\theta-\theta_{j}\leq\frac{1+\varepsilon
{H\pi\rho\ell_{j}}\ .
\]
Observe that (\ref{alfa}) and (\ref{alfa A}) yiel
\begin{align*}
&  |\sin(\pi\rho\ell_{j}\sin\left(  \theta-\theta_{j}\right)  )|-|\sin(\pi
\rho\ell_{j+n}\sin\left(  \theta-\theta_{j+n}\right)  )|\\
&  \geq\sin(\pi\rho\ell_{j}\sqrt{\alpha}\left(  \theta-\theta_{j}\right)
)-\sin(\pi\rho\ell_{j+n}\left(  \theta-\theta_{j+n}\right)  )\geq\frac{\alpha
}{H}-\frac{1+\varepsilon}{H}\frac{\ell_{j+n}}{\ell_{j}}=:a_{j}>0\ .
\end{align*}
Henc
\begin{align}
&  \left\vert \frac{\sin(\pi\rho\ell_{j}\sin(\theta-\theta_{j}))}{\sin
(\theta-\theta_{j})}\cos(\theta-\theta_{j})e^{-\pi i\rho\Theta\cdot\left(
P_{j+1}+P_{j}\right)  }\right. \label{19}\\
&  \left.  +\frac{\sin(\pi\rho\ell_{j+n}\sin(\theta-\theta_{j}))}{\sin
(\theta-\theta_{j})}\cos(\theta-\theta_{j+n})e^{-\pi i\rho\Theta\cdot\left(
P_{j+n+1}+P_{j+n}\right)  }\right\vert \nonumber\\
&  \geq\frac{|\cos(\theta-\theta_{j})|}{|\sin(\theta-\theta_{j})|}\left(
|\sin(\pi\rho\ell_{j}\sin(\theta-\theta_{j}))|-|\sin(\pi\rho\ell_{j+n
\sin(\theta-\theta_{j+n}))|\right) \nonumber\\
&  \geq a_{j}\ \frac{|\cos(\theta-\theta_{j})|}{|\sin(\theta-\theta_{j
)|}\ .\nonumber
\end{align}
We use the previous estimates to evaluate the last integral in (\ref{generale
) in a neighborhood of $\theta_{j}$ and therefore obtain an estimate from
below of $\Vert D_{P}^{\rho}\Vert_{L^{2}\left(  SO(2)\times\mathbb{T
^{2}\right)  }$. By the arguments in \cite[Theorem 2.3]{BCT} or \cite[Lemma
10.6]{T}, the contribution of all the sides $P_{h}P_{h+1}$ (with $h\neq j$ and
$h\neq j+n$) to the term $\Vert D_{P}^{\rho}\Vert_{L^{2}\left(  SO(2)\times
\mathbb{T}^{2}\right)  }$ is $\mathcal{O}\left(  1\right)  $. Then
(\ref{dalbassocondecad}), (\ref{divergenza2}) and (\ref{19}) yiel
\[
\Vert D_{P}^{\rho}\Vert_{L^{2}\left(  SO(2)\times\mathbb{T}^{2}\right)
^{2}\geq c\ \int_{\frac{1}{H\pi\rho\ell_{j}}}^{\frac{1+\varepsilon}{H\pi
\rho\ell_{j}}}\frac{\cos^{2}\theta}{\sin^{2}\theta}\ d\theta+c_{1}\geq
c\int_{\frac{1}{H\pi\rho\ell_{j}}}^{\frac{1+\varepsilon}{H\pi\rho\ell_{j}
}\frac{d\theta}{\theta^{2}}+c_{1}\geq c_{2}\ \rho\ .
\]

\end{proof}

\begin{proof}
[Proof of Lemma \ref{lemma3}]We can assume that $P\sim\{P_{h}\}_{h=1}^{2n}$ is
a convex polygon such that for every $h=1,\ldots,n$ the sides $P_{h}P_{h+1}$
and $P_{h+n}P_{h+n+1}$ are parallel and of the same length (that is, $\ell
_{h}=\ell_{h+n}$, $\ \tau_{h}=-\tau_{h+n}$, \ $\nu_{h}=-\nu_{h+n}$). Then we
may assume that $P$ is symmetric about the origin. As $P$ cannot be inscribed
in a circle, there exists an index $1\leq j\leq n$ such that the two opposite
equal and parallel sides $P_{j}P_{j+1}$ and $P_{j+n}P_{j+n+1}$ are not the
sides of a rectangle. Then $P_{j}+P_{j+1}$ is not orthogonal to $P_{j+1
-P_{j}$. Let $\phi_{j}\in\lbrack\theta_{j}-\pi,\theta_{j}]$ be defined by
\[
P_{j+1}+P_{j}=\mathcal{L}_{j}(\cos\phi_{j},\sin\phi_{j})\ .
\]
Since $\tau_{j}=(\cos\theta_{j},\sin\theta_{j})$ and $\nu(j)=(\cos(\theta
_{j}-\frac{\pi}{2}),\sin(\theta_{j}-\frac{\pi}{2}))$, see (\ref{thetaj}) and
(\ref{nu j}), we have $\phi_{j}-\theta_{j}\neq-\frac{\pi}{2}$. We put
$\varphi_{j}:=\phi_{j}-\theta_{j}$. Then
\[
\varphi_{j}\in\lbrack-\pi,0]\setminus\{-\frac{\pi}{2}\}\ .
\]
Again we need to find a lower bound for the last integral in (\ref{generale}).
As in the previous proof it is enough to conside
\begin{align*}
F_{j}(\theta)  &  :=\sum_{h\in\{j,j+n\}}\frac{\sin(\theta-\theta_{h})
{\cos(\theta-\theta_{h})}\sin(\pi\rho\ell_{j}\cos(\theta-\theta_{h}))e^{-\pi
i\rho\Theta\cdot(P_{h+1}+P_{h})}\\
&  =\frac{\sin(\theta-\theta_{j})}{\cos(\theta-\theta_{j})}\sin(\pi\rho
\ell_{j}\cos(\theta-\theta_{j}))\left[  e^{-\pi i\rho\mathcal{L}_{j
\cos(\theta-\phi_{j})}-e^{\pi i\rho\mathcal{L}_{j}\cos(\theta-\phi_{j
)}\right] \\
&  =-2i\ \frac{\sin(\theta-\theta_{j})}{\cos(\theta-\theta_{j})}\sin(\pi
\rho\ell_{j}\cos(\theta-\theta_{j}))\sin(\pi\rho\mathcal{L}_{j}\cos
(\theta-\phi_{j}))\ .
\end{align*}
We writ
\begin{align*}
&  \int_{0}^{2\pi}\left\vert \frac{\sin(\theta-\theta_{j})}{\cos(\theta
-\theta_{j})}\sin(\pi\rho\ell_{j}\cos(\theta-\theta_{j}))\sin(\pi
\rho\mathcal{L}_{j}\cos(\theta-\phi_{j}))\right\vert ^{2}\ d\theta\\
&  =\int_{0}^{2\pi}\left\vert \frac{\sin\left(  \pi\rho\ell_{j}\sin
\theta\right)  }{\sin\theta}\cos\theta\sin(\pi\rho\mathcal{L}_{j}\sin
(\theta-\varphi_{j}))\right\vert ^{2}\ d\theta\ .
\end{align*}
We shall integrate $\theta$ in a neighborhood of $0$ (actually $0\leq
\theta\leq1$ suffices). As for $\varphi_{j}$ we first assume $\varphi_{j
\in(-\frac{\pi}{2},0]$. Then $\cos\varphi_{j}>0$ and $\sin\varphi_{j}\leq0$.
Let $0<\gamma<1$ satisfy $\cos\varphi_{j}>\gamma$. In order to prove that
$\left\vert \sin(\pi\rho\mathcal{L}_{j}\sin(\theta-\varphi_{j})\right\vert
\geq c$ we consider two cases.

\noindent\emph{Case 1:} $|\sin(\pi\rho\mathcal{L}_{j}\sin\varphi_{j
)|>\gamma/2$.

\noindent We need to bound $\sin(\theta-\varphi_{j})-\left\vert \sin
\varphi_{j}\right\vert $. Since $\sin\varphi_{j}\leq0$ one ha
\[
\frac{\theta}{2}\cos\varphi_{j}+\left[  1-\frac{\theta^{2}}{2}\right]
|\sin\varphi_{j}|\leq\sin\theta\cos\varphi_{j}-\cos\theta\sin\varphi_{j
\leq\theta\cos\varphi_{j}+|\sin\varphi_{j}|\ .
\]
Therefore
\begin{equation}
\frac{\theta}{2}\gamma-\frac{\theta^{2}}{2}\leq\sin(\theta-\varphi
_{j})-\left\vert \sin\varphi_{j}\right\vert \leq\theta\ . \label{newast
\end{equation}
Let $\rho\geq1$ and assume
\[
\frac{\gamma}{8\pi\rho\mathcal{L}_{j}}\leq\theta\leq\frac{\gamma}{4\pi
\rho\mathcal{L}_{j}}\ .
\]
We recall that $\mathcal{L}_{j}\geq1$. Again we have to estimate $\sin
(\theta-\varphi_{j})-\left\vert \sin\varphi_{j}\right\vert $. By
(\ref{newast}) we hav
\[
0<\frac{\gamma}{16\pi\rho\mathcal{L}_{j}}-\frac{\gamma^{2}}{32(\pi
\rho\mathcal{L}_{j})^{2}}<\sin(\theta-\varphi_{j})-|\sin\varphi_{j}|\leq
\frac{\gamma}{4\pi\rho\mathcal{L}_{j}}\ .
\]
Therefore
\begin{equation}
0<\pi\rho\mathcal{L}_{j}\sin(\theta-\varphi_{j})-\pi\rho\mathcal{L}_{j
|\sin\varphi_{j}|\leq\frac{\gamma}{4}\ . \label{astetrig1
\end{equation}
Hence the assumption of \emph{Case 1} and (\ref{astetrig1}) yield
\begin{align*}
&  \left\vert \sin(\pi\rho\mathcal{L}_{j}\sin(\theta-\varphi_{j}))\right\vert
\\
&  =\left\vert \sin(\pi\rho\mathcal{L}_{j}\left[  \sin(\theta-\varphi
_{j})+\sin\varphi_{j}\right]  -\pi\rho\mathcal{L}_{j}\sin\varphi
_{j})\right\vert \\
&  =\left\vert \sin\left(  \pi\rho\mathcal{L}_{j}\left[  \sin(\theta
-\varphi_{j})-\left\vert \sin\varphi_{j}\right\vert \right]  \right)
\cos\left(  \pi\rho\mathcal{L}_{j}\sin\varphi_{j}\right)  \right. \\
&  \left.  -\cos\left(  \pi\rho\mathcal{L}_{j}\left[  \sin(\theta-\varphi
_{j})-\left\vert \sin\varphi_{j}\right\vert \right]  \right)  \sin\left(
\pi\rho\mathcal{L}_{j}\sin\varphi_{j}\right)  \right\vert \\
&  \geq\left\vert \sin(\pi\rho\mathcal{L}_{j}\sin\varphi_{j})\right\vert
\left\vert \cos(\pi\rho\mathcal{L}_{j}\sin(\theta-\varphi_{j})-\pi
\rho\mathcal{L}_{j}|\sin\varphi_{j}|)\right\vert \\
&  -\left\vert \sin(\pi\rho\mathcal{L}_{j}\sin(\theta-\varphi_{j})-\pi
\rho\mathcal{L}_{j}|\sin\varphi_{j}|)\right\vert \\
&  >\frac{\gamma}{2}\left[  1-\frac{\gamma^{2}}{32}\right]  -\frac{\gamma
{4}\\
&  >\frac{\gamma}{5}\ .
\end{align*}

\noindent\emph{Case 2:} $|\sin(\pi\rho\mathcal{L}_{j}\sin\varphi_{j
)|\leq\gamma/2$.

\noindent Let $\rho$ be large so that $0\leq\theta\leq\frac{3}{2\pi
\rho\mathcal{L}_{j}}$ implies $\sin\theta\geq(1-\delta)\theta$, with
$\delta<1/20$. Then for
\begin{equation}
\frac{1}{\pi\rho\mathcal{L}_{j}}\leq\theta\leq\frac{3}{2\pi\rho\mathcal{L
_{j}} \label{newpo
\end{equation}
we hav
\[
\theta(1-\delta)\gamma+\left[  1-\frac{\theta^{2}}{2}\right]  |\sin\varphi
_{j}|\leq\sin\theta\cos\varphi_{j}-\cos\theta\sin\varphi_{j}\leq\theta
+|\sin\varphi_{j}|
\]
and
\begin{equation}
\theta\gamma(1-\delta)-\frac{\theta^{2}}{2}\leq\sin(\theta-\varphi_{j
)-|\sin\varphi_{j}|\leq\theta\ . \label{newsta
\end{equation}
For $\rho$ large enough we have $\frac{9}{8(\pi\rho\mathcal{L}_{j})^{2}
<\frac{\gamma\delta}{\pi\rho\mathcal{L}_{j}}$. Then (\ref{newpo}) and
(\ref{newsta}) yiel
\begin{align*}
\frac{\gamma(1-2\delta)}{\pi\rho\mathcal{L}_{j}}  &  <\frac{\gamma(1-\delta
)}{\pi\rho\mathcal{L}_{j}}-\frac{9}{8(\pi\rho\mathcal{L}_{j})^{2}}\leq
\theta\gamma\left(  1-\delta\right)  -\frac{\theta^{2}}{2}\\
&  <\sin(\theta-\varphi_{j})-|\sin\varphi_{j}|\leq\theta\leq\frac{3}{2\pi
\rho\mathcal{L}_{j}
\end{align*}
and
\begin{equation}
\gamma(1-2\delta)<\pi\rho\mathcal{L}_{j}\sin(\theta-\varphi_{j})-\pi
\rho\mathcal{L}_{j}|\sin\varphi_{j}|\leq\frac{3}{2}\ . \label{caso2tr
\end{equation}
We choose $\gamma$ small enough so that
\[
\sin(\gamma(1-2\delta))\geq(1-2\delta)^{2}\gamma\text{ \ \ \ and\ \ \ \
\gamma^{2}/4<2\delta\ .
\]
Then (\ref{caso2tr}) and the assumption of \emph{Case 2} yield
\begin{align}
&  \left\vert \sin(\pi\rho\mathcal{L}_{j}\sin(\theta-\varphi_{j}))\right\vert
=\left\vert \sin(\pi\rho\mathcal{L}_{j}\left[  \sin(\theta-\varphi_{j
)-\sin\varphi_{j}\right]  +\pi\rho\mathcal{L}_{j}\sin\varphi_{j})\right\vert
\nonumber\\
&  \geq\left\vert \cos(\pi\rho\mathcal{L}_{j}\sin\varphi_{j})\right\vert
\left\vert \sin(\pi\rho\mathcal{L}_{j}\sin(\theta-\varphi_{j})-\pi
\rho\mathcal{L}_{j}|\sin\varphi_{j}|)\right\vert -\left\vert \sin(\pi
\rho\mathcal{L}_{j}\sin\varphi_{j})\right\vert \nonumber\\
&  \geq\gamma(1-2\delta)^{2}\sqrt{1-\frac{\gamma^{2}}{4}}-\frac{\gamma
{2}>\gamma\left[  (1-2\delta)^{5/2}-\frac{1}{2}\right]  >\frac{\gamma
{4}\ .\nonumber
\end{align}
\emph{Case 1} and \emph{Case 2} prove that for a suitable choice of
$0<\gamma<1$, such that $\cos\varphi_{j}>\gamma$, there exist $0<\alpha<\beta$
such that for $\frac{\alpha}{\pi\rho\mathcal{L}_{j}}\leq\theta\leq\frac{\beta
}{\pi\rho\mathcal{L}_{j}}$ and $\rho$ large enough we have
\begin{equation}
\left\vert \sin(\pi\rho\mathcal{L}_{j}\sin(\theta-\varphi_{j}))\right\vert
>\frac{\gamma}{5}\ . \label{stimaseno
\end{equation}
If $\varphi_{j}\in\lbrack-\pi,-\frac{\pi}{2})$ we have $\cos\varphi_{j}<0$ and
$\sin\varphi_{j}\leq0$. Then for $0\leq\theta<1$ we have
\[
-\theta+\left[  1-\frac{\theta^{2}}{2}\right]  |\sin\varphi_{j}|\leq\sin
\theta\cos\varphi_{j}-\cos\theta\sin\varphi_{j}\leq-\sin\theta|\cos\varphi
_{j}|+|\sin\varphi_{j}|\ .
\]
Hence, for a positive constant $K$,
\[
\sin\theta|\cos\varphi_{j}|\leq|\sin\varphi_{j}|-\sin(\theta-\varphi_{j})\leq
K\ \theta\ .
\]
If we choose a suitable constant $\gamma>0$ such that $|\cos\varphi
_{j}|>\gamma$, we can prove as for the case $\varphi_{j}\in(-\frac{\pi}{2},0]$
that (\ref{stimaseno}) still holds for $\frac{\alpha}{\pi\rho\mathcal{L}_{j
}\leq\theta\leq\frac{\beta}{\pi\rho\mathcal{L}_{j}}$, with $0<\alpha<\beta$
and $\rho$ large enough. Then (\ref{stimaseno}) yield
\begin{align*}
&  \int_{0}^{2\pi}\left\vert F_{j}(\theta)\right\vert ^{2}d\theta\geq
\int_{\frac{\alpha}{\pi\rho\mathcal{L}_{j}}}^{\frac{\beta}{\pi\rho
\mathcal{L}_{j}}}\left\vert \frac{\sin(\pi\rho\ell_{j}\sin\theta)}{\sin\theta
}\cos\theta\sin(\pi\rho\mathcal{L}_{j}\sin(\theta-\varphi_{j}))\right\vert
^{2}d\theta\\
&  \geq c\ \gamma^{2}\int_{\frac{\alpha}{\pi\rho\mathcal{L}_{j}}}^{\frac
{\beta}{\pi\rho\mathcal{L}_{j}}}\left\vert \frac{\sin(\pi\rho\ell_{j
\sin\theta)}{\sin\theta}\right\vert ^{2}d\theta\geq c_{1}\ \rho\int_{c_{1
}^{c_{2}}\left\vert \frac{\sin(t)}{t}\right\vert ^{2}dt\geq c_{2}\ \rho.
\end{align*}
This ends the proof.
\end{proof}

The proof of Theorem \ref{maintheorem} will be complete after the proof of
Proposition \ref{propeps}. We need a simultaneous approximation lemma from
\cite{PS}.

\begin{lemma}
\label{diri} Let $r_{1},r_{2},...,r_{n}\in\mathbb{R}$. For every positive
integer $j$ there exists $j\leq q\leq j^{n+1}$ such that $\left\Vert
r_{s}q\right\Vert <j^{-1}$ for any $1\leq s\leq n$, where $\left\Vert
x\right\Vert $ denotes the distance of a real number $x$ from the integers.
\end{lemma}

\begin{proof}
[Proof of Proposition \ref{propeps}]Let $P\sim\{P_{j}\}_{j=1}^{2n}$ be a
polygon in $\mathfrak{P}$. For every positive integer $u$ let
\[
A_{u}^{j}:=\{k\in\mathbb{Z}^{2}:0<\mathcal{L}_{j}|k|\leq u^{2
\}\ \ \ \text{for }j=1,\ldots,n\ ,\ \ \ \ \ \ \ \ \text{\ }A_{u
:=\bigcup_{j=1}^{n}A_{u}^{j}\ .
\]
Observe that $\operatorname*{card}(A_{u}^{j})\leq4u^{4}$ and therefore
$\operatorname*{card}(A_{u})\leq4nu^{4}$. By Lemma \ref{diri} there exists a
sequence $\left\{  \rho_{u}\right\}  _{u=1}^{+\infty}$ of positive integers
such that, for$\ $every $k\in A_{u}$ and every $j=1,\ldots,n$,
\begin{equation}
u\leq\rho_{u}\leq u^{4nu^{4}+1}\ ,\ \ \ \ \ |\sin(\pi\rho_{u}|k|\mathcal{L
_{j})|<1/u\ . \label{i e ii
\end{equation}
Observe that (\ref{i e ii}) implie
\begin{equation}
u\geq c_{\varepsilon}\log^{\frac{1}{4+\varepsilon}}(\rho_{u}) \label{iii
\end{equation}
for every $\varepsilon>0$. For any $1\leq j\leq n$ and $k\in A_{u}^{j}$ we
split the integral in (\ref{final}) into several parts
\[
E_{1,j,|k|}^{\rho}:=\int_{0}^{\left(  8\rho_{u}|k|\right)  ^{-1}}\left\vert
\frac{\sin(\pi\rho_{u}|k|\ell_{j}\sin\theta)}{\sin\theta}\sin(\pi\rho
_{u}|k|\mathcal{L}_{j}\cos\theta)\right\vert ^{2}d\theta\ .
\]
For $0\leq\theta\leq\left(  8\rho_{u}|k|\right)  ^{-1}$ we have $0\leq
1-\cos\theta\leq\left(  128\rho_{u}^{2}|k|^{2}\right)  ^{-1}$. Then
(\ref{i e ii}) yiel
\begin{align}
&  \left\vert \sin(\pi\rho_{u}|k|\mathcal{L}_{j}\cos\theta)\right\vert
\label{33}\\
&  =\left\vert \sin(\pi\rho_{u}|k|\mathcal{L}_{j}\left[  \cos\theta
-1+1\right]  )\right\vert \nonumber\\
&  \leq\left\vert \sin(\pi\rho_{u}|k|\mathcal{L}_{j}(\cos\theta-1))\cos
(\pi\rho_{u}|k|\mathcal{L}_{j})\right\vert \nonumber\\
&  +\left\vert \sin(\pi\rho_{u}|k|\mathcal{L}_{j})\cos(\pi\rho_{u
|k|\mathcal{L}_{j}(\cos\theta-1)\right\vert \nonumber\\
&  \leq\left\vert \sin(\pi\rho_{u}|k|\mathcal{L}_{j}(1-\cos\theta))\right\vert
+\left\vert \sin(\pi\rho_{u}|k|\mathcal{L}_{j})\right\vert \nonumber\\
&  \leq\frac{\pi\mathcal{L}_{j}}{128\rho_{u}|k|}+\left\vert \sin(\pi\rho
_{u}|k|\mathcal{L}_{j})\right\vert \nonumber\\
&  \leq c\ \frac{1}{u}\ .\nonumber
\end{align}
By (\ref{33}) we obtain
\begin{align*}
E_{1,j,|k|}^{\rho}  &  \leq c\ \frac{1}{u^{2}}\int_{0}^{\left(  8\rho
_{u}|k|\right)  ^{-1}}\left\vert \frac{\sin(\pi\rho_{u}|k|\ell_{j}\sin\theta
)}{\sin\theta}\right\vert ^{2}d\theta\\
&  \leq c_{1}\ \frac{|k|\rho_{u}}{u^{2}}\int_{0}^{1}\left\vert \frac{\sin
(t)}{t}\right\vert ^{2}dt\leq c_{2}\ \frac{|k|\rho_{u}}{u^{2}}\ .
\end{align*}
Let
\[
E_{2,j,|k|}^{\rho}:=\int_{\left(  8\rho_{u}|k|\right)  ^{-1}}^{\left(
8u^{1/4}\rho_{u}^{1/2}\left\vert k\right\vert ^{1/2}\right)  ^{-1}}\left\vert
\frac{\sin(\pi\rho_{u}|k|\ell_{j}\sin\theta)}{\sin\theta}\sin(\pi\rho
_{u}|k|\mathcal{L}_{j}\cos\theta)\right\vert ^{2}d\theta\ .
\]
For $\left(  8\rho_{u}|k|\right)  ^{-1}\leq\theta\leq\left(  8\ell^{1/4
\rho_{u}^{1/2}\left\vert k\right\vert ^{1/2}\right)  ^{-1}$ we have
\[
\frac{1}{2000\rho_{u}^{2}|k|^{2}}\leq2\sin^{2}\left(  \theta/2\right)
=1-\cos\theta\leq\frac{1}{128u^{1/2}\rho_{u}|k|}\ .
\]
As in (\ref{33}) we obtain
\begin{align*}
\left\vert \sin(\pi\rho_{u}|k|\mathcal{L}_{j}\cos\theta)\right\vert  &
\leq\left\vert \sin(\pi\rho_{u}|k|\mathcal{L}_{j}(1-\cos\theta))\right\vert
+\left\vert \sin(\pi\rho_{u}|k|\mathcal{L}_{j})\right\vert \\
&  \leq\frac{\pi\mathcal{L}_{j}}{128u^{1/2}}+\frac{1}{u}\leq c\ u^{-1/2
\end{align*}
and then
\begin{align*}
E_{2,j,|k|}^{\rho}  &  \leq c\ \frac{1}{u}\int_{\left(  8\rho_{u}|k|\right)
^{-1}}^{\left(  8u^{1/4}\rho_{u}^{1/2}\left\vert k\right\vert ^{1/2}\right)
^{-1}}\left\vert \frac{\sin(\pi\rho_{u}|k|\ell_{j}\sin\theta)}{\sin\theta
}\right\vert ^{2}d\theta\\
&  \leq c_{1}\ \frac{1}{u}\int_{\left(  8\rho_{u}|k|\right)  ^{-1}}^{\left(
8u^{1/4}\rho_{u}^{1/2}\left\vert k\right\vert ^{1/2}\right)  ^{-1}
\frac{d\theta}{\theta^{2}}\leq c_{2}\ \frac{\rho_{u}|k|}{u}\ .
\end{align*}
Let $1/4<\lambda<1/2$ and let
\[
E_{3,j,|k|}^{\rho}:=\int_{\left(  8u^{1/4}\rho_{u}^{1/2}\left\vert
k\right\vert ^{1/2}\right)  ^{-1}}^{\lambda}\left\vert \frac{\sin(\pi\rho
_{u}|k|\ell_{j}\sin\theta)}{\sin\theta}\sin(\pi\rho_{u}|k|\mathcal{L}_{j
\cos\theta)\right\vert ^{2}d\theta\ .
\]
We hav
\[
E_{3,j,|k|}^{\rho}\leq\int_{\left(  8u^{1/4}\rho_{u}^{1/2}\left\vert
k\right\vert ^{1/2}\right)  ^{-1}}^{\lambda}\frac{d\theta}{\theta^{2}
\leq8u^{1/4}\rho_{u}^{1/2}\left\vert k\right\vert ^{1/2}\ .
\]
Finally we have
\[
E_{4,j,|k|}^{\rho}:=\int_{\lambda}^{\frac{\pi}{2}}\left\vert \frac{\sin
(\pi\rho_{u}|k|\ell_{j}\sin\theta)}{\sin\theta}\sin(\pi\rho_{u}|k|\mathcal{L
_{j}\cos\theta)\right\vert ^{2}d\theta\leq c\ .
\]
By the above estimates, (\ref{final}), (\ref{i e ii})\ and (\ref{iii}) we hav
\begin{align*}
&  \Vert D_{P}^{\rho_{u}}\Vert_{L^{2}\left(  SO(2)\times\mathbb{T}^{2}\right)
}^{2}\\
&  \leq c\ \rho_{u}\sum_{k\in A_{u}}\frac{1}{|k|^{3}}\left(  \frac{1}{u^{2
}+\frac{1}{u}+u^{1/4}\rho_{u}^{-1/2}\left\vert k\right\vert ^{-1/2}+\rho
_{u}^{-1}\left\vert k\right\vert ^{-1}\right) \\
&  +c_{1}\ \sum_{k\notin A_{u}}\frac{1}{|k|^{4}}\int_{0}^{\pi/2}\left\vert
\frac{\sin(\pi\rho_{u}|k|\ell_{j}\sin\theta)}{\sin\theta}\right\vert
^{2}\ d\theta\\
&  \leq c\ \rho_{u}\sum_{0\neq k\in A_{u}}\frac{1}{|k|^{3}}u^{-1/4}\\
&  +c_{1}\ \sum_{\left\vert k\right\vert >c_{1}u^{2}}\frac{1}{|k|^{4}}\left(
\int_{0}^{\left(  \rho_{u}\left\vert k\right\vert \right)  ^{-1/2}}\left(
\rho_{u}\left\vert k\right\vert \right)  \ d\theta+\int_{\left(  \rho
_{u}\left\vert k\right\vert \right)  ^{-1/2}}^{\pi/2}\frac{1}{\theta^{2
}\ d\theta\right) \\
&  \leq c_{\varepsilon}\ \rho_{u}\sum_{0\neq k\in\mathbb{Z}^{2}}\frac
{1}{|k|^{3}}\log^{-\frac{1}{16+\varepsilon}}\left(  \rho_{u}\right)
+c\ \sum_{\left\vert k\right\vert >c_{1}u^{2}}\frac{1}{|k|^{4}}\left(
\rho_{u}\left\vert k\right\vert \right)  ^{1/2}\\
&  \leq c_{\varepsilon}\ \rho_{u}\log^{-\frac{1}{16+\varepsilon}}\left(
\rho_{u}\right)  +c\ \rho_{u}^{1/2}\int_{\{t\in\mathbb{R}^{2}:|t|>c_{1
u^{2}\}}\frac{1}{\left\vert t\right\vert ^{7/2}}\ dt\\
&  \leq c_{\varepsilon}\ \rho_{u}\log^{-\frac{1}{16+\varepsilon}}\left(
\rho_{u}\right)  \ .
\end{align*}

\end{proof}

We now turn to the proof of Proposition \ref{stimaepsilon}, which depends on
the following lemma proved by L. Parnovski and A. Sobolev \cite{PS}.

\begin{lemma}
\label{PSlemma} For any $\varepsilon>0$ there exist $\rho_{0}\geq1$ and
$0<\alpha<1/2$ such that for every $\rho\geq\rho_{0}$ there exists
$k\in\mathbb{Z}^{d}$ such that $|k|\leq\rho^{\varepsilon}$ and $\Vert
\rho|k|\Vert\geq\alpha$, where $\Vert x\Vert$ is the distance of a real number
$x$ from the integers.
\end{lemma}

\begin{proof}
[Proof of Proposition \ref{stimaepsilon}]Let $P\sim\{P_{j}\}_{j=1}^{2n}$ be a
polygon in $\mathfrak{P}$. Let $\varepsilon>0$ and let $j\in\left\{
1,2,\ldots,n\right\}  $. By Lemma \ref{PSlemma} there exist $\rho_{0}\geq1$
and $0<a<1/2$ such that for any $\rho\geq\rho_{0}$ there is $\widetilde{k
\in\mathbb{Z}^{2}$ such that $|\widetilde{k}|\leq\rho^{\varepsilon/3}$ and
$|\sin(\pi\rho|\widetilde{k}|\mathcal{L}_{j})|>a$. Then we consider the
interva
\begin{equation}
\theta_{j}\leq\theta\leq\theta_{j}+\frac{1}{\pi\rho|\widetilde{k}|}\ .
\label{interval j
\end{equation}
We have
\[
0\leq1-\cos(\theta-\theta_{j})\leq\frac{1}{2(\pi\rho|\widetilde{k}|)^{2}}\ .
\]
Then for large $\rho$ we have
\begin{align}
&  |\sin(\pi\rho|\widetilde{k}|\mathcal{L}_{j}\cos(\theta-\theta
_{j}))|\label{epsi1}\\
&  \geq|\sin(\pi\rho|\widetilde{k}|\mathcal{L}_{j})||\cos(\pi\rho
|\widetilde{k}|\mathcal{L}_{j}(1-\cos(\theta-\theta_{j})))|\nonumber\\
&  -|\sin(\pi\rho|\widetilde{k}|\mathcal{L}_{j}(1-\cos(\theta-\theta
_{j})))|\nonumber\\
&  \geq c\ \left[  1-\frac{\mathcal{L}_{j}^{2}}{8(\pi\rho|\widetilde{k}|)^{2
}\right]  -\frac{\mathcal{L}_{j}}{2\pi\rho|\widetilde{k}|}>c_{1}\ .\nonumber
\end{align}
As before the sides non parallel to $P_{j}P_{j+1}$ give a bounded contribution
to the integration of $\left\vert D_{P_{n}}^{\rho}\right\vert ^{2}$ over the
interval in (\ref{interval j}). Finally (\ref{epsi1}) yield
\begin{align*}
&  \Vert D_{P_{n}}^{\rho}\Vert_{L^{2}\left(  SO(2)\times\mathbb{T}^{2}\right)
}^{2}\\
&  \geq c+c_{1}\ \frac{1}{|\widetilde{k}|^{4}}\int_{\theta_{j}}^{\theta
_{j}+1/(\pi\rho|\widetilde{k}|)}\left\vert \frac{\sin(\pi\rho|\widetilde
{k}|\ell_{j}\sin(\theta-\theta_{j}))}{\sin(\theta-\theta_{j})}\right\vert
^{2}\\
&  \times|\sin(\pi\rho|\widetilde{k}|\mathcal{L}_{j}\cos(\theta-\theta
_{j}))|^{2}|\cos(\theta-\theta_{j})|^{2}\ d\theta\\
&  \geq c+c_{2}\ \frac{1}{|\widetilde{k}|^{4}}\int_{0}^{1/(\pi\rho
|\widetilde{k}|)}\left\vert \frac{\sin(\pi\rho|\widetilde{k}|\ell_{j
\sin\theta)}{\sin\theta}\right\vert ^{2}\ d\theta\\
&  \geq c+c_{3}\frac{\rho}{|\widetilde{k}|^{3}}\\
&  \geq c_{4}\ \rho^{1-\varepsilon}\ .
\end{align*}

\end{proof}

The proofs of Lemmas \ref{lemma1}, \ref{lemma2}, and \ref{lemma3} actually
show that $\left\Vert \widehat{\chi_{P}}\left(  \rho\cdot\right)  \right\Vert
_{L^{2}\left(  \Sigma_{1}\right)  }\geq c\ \rho^{-3/2}$ whenever
$P\notin\mathfrak{P}$. Hence Theorem \ref{maintheorem} and Proposition
\ref{prop} readily yield the following result.

\begin{corollary}
Let $P$ be a polygon in $\mathbb{R}^{2}$. Then $P$ satisfie
\[
\left\Vert \widehat{\chi_{P}}\left(  \rho\cdot\right)  \right\Vert
_{L^{2}\left(  \Sigma_{1}\right)  }\geq c\ \rho^{-3/2
\]
if and only if $P\notin\mathfrak{P}$.
\end{corollary}

\bigskip

The results in this paper (apart from Lemma \ref{lemma1}) seem to be tailored
for the planar case. A different (perhaps simpler) approach might be necessary
in order to deal with the multi-dimensional cases.