4
4
// Copyright (C) 2009 Benoit Jacob <jacob.benoit.1@gmail.com>
5
5
// Copyright (C) 2009 Gael Guennebaud <gael.guennebaud@inria.fr>
7
// Eigen is free software; you can redistribute it and/or
8
// modify it under the terms of the GNU Lesser General Public
9
// License as published by the Free Software Foundation; either
10
// version 3 of the License, or (at your option) any later version.
12
// Alternatively, you can redistribute it and/or
13
// modify it under the terms of the GNU General Public License as
14
// published by the Free Software Foundation; either version 2 of
15
// the License, or (at your option) any later version.
17
// Eigen is distributed in the hope that it will be useful, but WITHOUT ANY
18
// WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
19
// FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License or the
20
// GNU General Public License for more details.
22
// You should have received a copy of the GNU Lesser General Public
23
// License and a copy of the GNU General Public License along with
24
// Eigen. If not, see <http://www.gnu.org/licenses/>.
7
// This Source Code Form is subject to the terms of the Mozilla
8
// Public License v. 2.0. If a copy of the MPL was not distributed
9
// with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
26
11
#ifndef EIGEN_JACOBI_H
27
12
#define EIGEN_JACOBI_H
29
16
/** \ingroup Jacobi_Module
31
18
* \class JacobiRotation
326
313
// both vectors are sequentially stored in memory => vectorization
327
314
enum { Peeling = 2 };
329
Index alignedStart = first_aligned(y, size);
316
Index alignedStart = internal::first_aligned(y, size);
330
317
Index alignedEnd = alignedStart + ((size-alignedStart)/PacketSize)*PacketSize;
332
319
const Packet pc = pset1<Packet>(j.c());
344
331
Scalar* EIGEN_RESTRICT px = x + alignedStart;
345
332
Scalar* EIGEN_RESTRICT py = y + alignedStart;
347
if(first_aligned(x, size)==alignedStart)
334
if(internal::first_aligned(x, size)==alignedStart)
349
336
for(Index i=alignedStart; i<alignedEnd; i+=PacketSize)