You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
70 lines
2.8 KiB
70 lines
2.8 KiB
/* Copyright (c) 2022, NVIDIA CORPORATION. All rights reserved. |
|
* |
|
* Redistribution and use in source and binary forms, with or without |
|
* modification, are permitted provided that the following conditions |
|
* are met: |
|
* * Redistributions of source code must retain the above copyright |
|
* notice, this list of conditions and the following disclaimer. |
|
* * Redistributions in binary form must reproduce the above copyright |
|
* notice, this list of conditions and the following disclaimer in the |
|
* documentation and/or other materials provided with the distribution. |
|
* * Neither the name of NVIDIA CORPORATION nor the names of its |
|
* contributors may be used to endorse or promote products derived |
|
* from this software without specific prior written permission. |
|
* |
|
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS ``AS IS'' AND ANY |
|
* EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE |
|
* IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR |
|
* PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR |
|
* CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, |
|
* EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, |
|
* PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR |
|
* PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY |
|
* OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT |
|
* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE |
|
* OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. |
|
*/ |
|
|
|
//! includes, system |
|
#include <algorithm> |
|
#include <cassert> |
|
#include <cmath> |
|
#include <fstream> |
|
#include <iomanip> |
|
#include <iostream> |
|
#include <map> |
|
#include <vector> |
|
|
|
// includes, projcet |
|
#include "matlab.h" |
|
|
|
// namespace, unnamed |
|
namespace { |
|
} // namespace |
|
|
|
/////////////////////////////////////////////////////////////////////////////// |
|
//! Write a tridiagonal, symmetric matrix in vector representation and |
|
//! it's eigenvalues |
|
//! @param filename name of output file |
|
//! @param d diagonal entries of the matrix |
|
//! @param s superdiagonal entries of the matrix (len = n - 1) |
|
//! @param eigenvals eigenvalues of the matrix |
|
//! @param indices vector of len n containing the position of the eigenvalues |
|
//! if these are sorted in ascending order |
|
//! @param n size of the matrix |
|
/////////////////////////////////////////////////////////////////////////////// |
|
void writeTridiagSymMatlab(const char *filename, float *d, float *s, float *eigenvals, const unsigned int n) |
|
{ |
|
std::ofstream file(filename, std::ios::out); |
|
|
|
// write diagonal entries |
|
writeVectorMatlab(file, "d", d, n); |
|
|
|
// write superdiagonal entries |
|
writeVectorMatlab(file, "s", s, n - 1); |
|
|
|
// write eigenvalues |
|
writeVectorMatlab(file, "eigvals", eigenvals, n); |
|
|
|
file.close(); |
|
}
|
|
|