R Under development (unstable) (2024-01-07 r85787 ucrt) -- "Unsuffered Consequences" Copyright (C) 2024 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > #****************************************************************************** > # Copyright (c) 2016, College of William & Mary > # 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 the College of William & Mary 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 AND CONTRIBUTORS "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 COLLEGE OF WILLIAM & MARY 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. > # > # PRIMME: https://github.com/primme/primme > # Contact: Andreas Stathopoulos, a n d r e a s _at_ c s . w m . e d u > #****************************************************************************** > # File: tests.R > # > # Purpose - Test functions provided by the library > # > #***************************************************************************** > > require(PRIMME) Loading required package: PRIMME > > # Test for dense and sparse matrices > > As <- list(diag(1:100), diag(as.complex(1:100))); > if (requireNamespace("Matrix", quietly = TRUE)) + As[[3]] <- Matrix::Matrix(diag(1:100), sparse=TRUE); > > for (i in 1:length(As)) { + A = As[[i]]; + d <- eigs_sym(A, 3); + stopifnot(all.equal(c(100,99,98), d$values, tolerance=1e-7)); + + d <- svds(A, 3); + stopifnot(all.equal(c(100,99,98), d$d, tolerance=1e-7)); + } > > proc.time() user system elapsed 0.89 0.17 1.00