1function C = asec (G) 2%ASEC inverse secant. 3% C = asec (G) is the inverse secant of each entry of G. 4% Since asec (0) is nonzero, the result is a full matrix. 5% C is complex if any (abs(G) < 1). 6% 7% See also GrB/sec, GrB/sech, GrB/asech. 8 9% SuiteSparse:GraphBLAS, Timothy A. Davis, (c) 2017-2021, All Rights Reserved. 10% SPDX-License-Identifier: GPL-3.0-or-later 11 12G = G.opaque ; 13type = gbtype (G) ; 14if (~gb_isfloat (type)) 15 type = 'double' ; 16end 17 18C = GrB (gb_trig ('acos', gbapply ('minv', gbfull (G, type)))) ; 19 20