Add vectorized variant for functions of complex variable, especially exp
维护者通常 1 天内回复
还没有人认领这个 Issue。
评估
调研方向
该 issue 没有指明文件或测试。先从复值基本函数指南和现有的向量化 pow 行为开始，然后追踪复标量函数如何与向量或矩阵 autodiff 交互。完成内容应涵盖所请求的 exp 向量化行为，并澄清 log 和三角函数更广泛的一致性问题。
由索引模型根据 Issue 内容生成。
描述
Currently, according to the guide https://mc-stan.org/docs/functions-reference/complex-valued_basic_functions.html it seems that basic arithmetic operators and pow are the only functions that accept complex vectors.
I am modelling some periodicities in data using Fourier series, and I find it the most compact (and computationally efficient) to express what I do in terms of complex numbers. However, it's not efficient to run a function over a matrix using a loop when there's autodiff.
The most striking case is the exp function, which is ubiquitous in complex analysis and I want to be able to run it elementwise over a vector or a matrix.
I figured out a possible workaround is to use pow(exp(1), z) which does work on vectorized z, but still I believe that exp is more fundamental than pow and it should be implemented, because the current state is confusing at best.
Of course then there are all the other functions which accept a real or complex scalar and accept vectorized real, but don't accept vectorized complex, which is inconsistent. This includes log and the trigonometric functions.
- 主要语言
- C++
- 星标
- 845
- 派生
- 225
- 平均合并
- 4 天 21 小时
- 30 天内合并 PR
- 13
环境准备
- 没有 Dockerfile 或 Docker Compose 文件
- 有 Pull Request 模板
- 阅读贡献指南
从这里开始
- 先读完整个 Issue，再读项目的贡献指南。
- 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
- Fork 仓库，在一个分支上完成修改。
- 提交 Pull Request，并在描述里引用这个 Issue 编号。
stan-dev/math 的其他 Issue
-
难度 3/5 1-2 天 新手友好度 54/100
维护者通常 1 天内回复
-
难度 4/5 3-5 天 新手友好度 54/100
维护者通常 1 天内回复
-
难度 5/5 一周以上 新手友好度 35/100
维护者通常 1 天内回复
-
numerics performance
难度 5/5 一周以上 新手友好度 35/100
stan-dev/math#3481 · 1 个 reaction ·
维护者通常 1 天内回复
-
`eigendecompose_sym`, `eigenvalues_sym` and `eigenvectors_sym` do not compile for complex matrices of `fvar<var>` or `fvar<fvar<double>>`可能重新可做 关联的 PR 已关闭且未合并。 未关闭bug
难度 3/5 1-2 天 新手友好度 55/100
维护者通常 1 天内回复
相似的 Issue
-
bug
难度 2/5 1-3 小时 新手友好度 76/100
pocoproject/poco#5664 ·
维护者通常 2 天内回复
-
难度 2/5 1-3 小时 新手友好度 78/100
维护者通常 1 天内回复
-
难度 2/5 1-3 小时 新手友好度 72/100
charmplusplus/charm#4210 ·
维护者通常 1 天内回复
-
难度 1/5 1 小时以内 新手友好度 88/100
cp-algorithms/cp-algorithms#1748 ·
维护者通常 1 天内回复
-
难度 2/5 1-3 小时 新手友好度 66/100
维护者通常 1 天内回复