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