全部问题(coding)

我有几个 python 包需要安装在不同的操作系统/环境中。这些包有依赖关系，其中一些像 Polars 一样，根据操作系统的不同需要不同的包，例如：MacOS（Darwin）上的 polars-lts-cpu 和所有其他操作系统上的 polars。

我用来setuptools创建whl文件，但安装的依赖项取决于创建 wheel 文件的操作系统。这是我的代码：

import platform
from setuptools import find_packages, setup

setup(
    ...
    install_requires=["glob2>=0.7",
                       "numpy>=1.26.4",
                       "polars>=1.12.0" if platform.system() != "Darwin" else "polars-lts-cpu>=1.12.0"]
    ...)

如上所述，此代码根据创建 wheel 文件的操作系统来安装 Polars 版本，而不是根据包的安装位置来安装。

我该如何修复此问题？

我正在对我的流量数据进行对数皮尔逊 III 分布拟合。拟合后，我想绘制观测值和拟合分布。

然而，这个数字并不是我所期望的：

1. x 轴标签没有显示 0.99 0.02 和 0.01，尽管我已将其设置breaks为 c(0.99、0.9、0.5、0.2、0.1、0.02、0.01)。

2.如何在图中绘制拟合的 Log Pearson III 分布？我提供了一个示例图供参考，谢谢您的帮助。

library(fitdistrplus)
library(PearsonDS)
library(ggplot2)

low_flows <- data.frame(Year = 1981:2010, Flow = c(0.03357143, 0.01328571, 0.02285714, 0.02657143, 0.04957143, 0.04085714, 0.16900000, 0.01057143, 
                                                   0.04128571, 0.10871429, 0.08771429, 0.09585714, 0.22057143, 0.11571428, 0.08300000, 0.11257143, 
                                                   0.13614286, 0.07742857, 0.09785714, 0.04728571, 0.04300000, 0.08385714, 0.02828571, 0.07271429, 
                                                   0.21428571, 0.10142857, 0.04400000, 0.10928571, 0.12471429, 0.31500000))


Log_mydata <- log10(low_flows$Flow)

m <- mean(Log_mydata)
v <- stats::var(Log_mydata)
g <- e1071::skewness(Log_mydata, type = 2)
my_shape <- (2/g)^2
my_scale <- sqrt(v)/sqrt(my_shape) * sign(g)
my_location <- m - my_scale * my_shape
start <- list(shape = my_shape, location = my_location, scale = my_scale)

dPIII <<- function(x, shape, location, scale) PearsonDS::dpearsonIII(x, 
                                                                     shape, location, scale, log = FALSE)
pPIII <<- function(q, shape, location, scale) PearsonDS::ppearsonIII(q, 
                                                                     shape, location, scale, lower.tail = TRUE, log.p = FALSE)
qPIII <<- function(p, shape, location, scale) PearsonDS::qpearsonIII(p, 
                                                                     shape, location, scale, lower.tail = TRUE, log.p = FALSE)

fit_lp3 <- fitdist(Log_mydata, distr = "PIII", start = start)

plotdata = low_flows
plotdata <- plotdata[order(plotdata$Flow), ]
plotdata$prob <- seq_len(length(plotdata$Flow))/(length(plotdata$Flow) + 1)

prob_scale_points = c(0.99, 0.9, 0.5, 0.2, 0.1, 0.02, 0.01)

ggplot(data = plotdata, aes(x = prob, y = Flow)) + 
  geom_point(size = 2) + 
  xlab("Probability") + 
  scale_x_continuous(transform = scales::probability_trans("norm", lower.tail = FALSE), 
                     breaks = prob_scale_points, 
                     sec.axis = dup_axis(name = "Return Period", 
                                         labels = function(x) {ifelse(1/x < 2, round(1/x, 2), round(1/x, 0))})) +
  theme_bw()  + 
  theme(panel.border = element_rect(colour = "black", fill = NA, linewidth = 1), 
        panel.grid = element_line(linewidth = 0.2), 
        axis.title = element_text(size = 12), axis.text = element_text(size = 10), 
        axis.title.x.top = element_text(size = 12), 
        legend.text = element_text(size = 10), legend.title = element_text(size = 10))

以下是我从上面的代码中得到的结果：

我的目标如下：

按照 Allan 的建议，Log Pearson III 可以正确绘制了！为了简化我的原始代码，我尝试使用scale_x_log10而不是scale_x_continuous。但是，Log Pearson III 的拟合线没有正确绘制。我想知道如何解决这个问题？

ggplot(data = plotdata, aes(x = prob, y = Flow)) + 
  geom_point(size = 2, color="blue") + 
  geom_function(fun = ~ 10^(qpearsonIII(.x, params = fit_lp3$estimate)),
                color = "black", linewidth = 1.5, alpha = 0.7) + 
  scale_x_log10(name = "Probability", 
                trans = "reverse", 
                breaks = prob_scale_points, 
                limits = rev(range(prob_scale_points)), 
                sec.axis = dup_axis(name = "Return Period", 
                                    labels = function(x) {ifelse(1/x < 2, round(1/x, 2), round(1/x, 0))})) + 
  scale_y_log10(breaks = seq(0.05, 0.4, 0.05)) + 
  theme_bw()  + 
  theme(panel.border = element_rect(colour = "black", fill = NA, linewidth = 1), 
        panel.grid = element_line(linewidth = 0.2), 
        axis.title = element_text(size = 12), axis.text = element_text(size = 10), 
        axis.title.x.top = element_text(size = 12), 
        legend.text = element_text(size = 10), legend.title = element_text(size = 10))

这是图。我们可以看到拟合的线不正确。

我正在尝试为我的 testrunner 帮助脚本设置一个 python wheel，以便从我的 Windows 机器上的任何地方轻松访问它。因此，我在 setup.py 中配置了一个控制台入口点。我可以在 entry_points.txt 中看到生成的入口点，但如果我尝试调用我的脚本，我会收到错误消息。

No module named testrunner.__main__; 'testrunner' is a package and cannot be directly executed

我的安装程序文件夹树如下所示

setup.py
README.md
LICENSE.txt
testrunner/
    __init__.py
    testrunner.py
    templates/
        testDesc.json

testrunner.py如下所示

def runTests():
    print("Hello World")
    #More not relevant code here

if __name__ == '__main__':
    runTests()

setup.py内容​

from setuptools import find_packages, setup
from pathlib import Path

# read the contents of your README file
this_directory = Path(__file__).parent
long_description = (this_directory / "README.md").read_text()
        
setup(
    name='testrunner',
    version='0.3.0',
    packages=find_packages(include=['testrunner']),

    description='C/C++ test runner',
    long_description=long_description,
    long_description_content_type='text/markdown',

    author='Me',

    license="Proprietary",
    license_files = ('LICENSE.txt',),
    
    entry_points={
        'console_scripts': ['testrunner = testrunner:main']
    },

    classifiers=[
        'Topic :: Software Development :: Build Tools',
        'Topic :: Software Development :: Compilers',

        'Private :: Do Not Upload',
        'Operating System :: Microsoft :: Windows',
        'Intended Audience :: Developers',
        'Intended Audience :: Science/Research ',
        'Intended Audience :: Education ',
        'Natural Language :: English',
        'License :: Other/Proprietary License',
        'Programming Language :: Python :: 3',
        'Programming Language :: Python :: 3.7',
        'Programming Language :: Python :: 3.8',
        'Programming Language :: Python :: 3.8',
        'Programming Language :: Python :: 3.9',
        'Programming Language :: Python :: 3.10',
    ],
    install_requires=[],
    package_data={'':['templates\\*.json']},
)

最后是__init__.py

from .testrunner import main, runTests

我使用以下命令构建轮子：py -m pip wheel --no-deps -w dist . 使用 pip 安装后，检查 site-packages 目录中的内容，执行py -m testrunner -h结果为No module named testrunner.__main__; 'testrunner' is a package and cannot be directly executed

我需要用零替换valuea 列中的最后两个值pl.DataFrame，从而我需要group_by该symbol列。

import polars as pl

df = pl.DataFrame(
    {"symbol": [*["A"] * 4, *["B"] * 4], "value": range(8)}
)

shape: (8, 2)
┌────────┬───────┐
│ symbol ┆ value │
│ ---    ┆ ---   │
│ str    ┆ i64   │
╞════════╪═══════╡
│ A      ┆ 0     │
│ A      ┆ 1     │
│ A      ┆ 2     │
│ A      ┆ 3     │
│ B      ┆ 4     │
│ B      ┆ 5     │
│ B      ┆ 6     │
│ B      ┆ 7     │
└────────┴───────┘

这是我的预期结果：

shape: (8, 2)
┌────────┬───────┐
│ symbol ┆ value │
│ ---    ┆ ---   │
│ str    ┆ i64   │
╞════════╪═══════╡
│ A      ┆ 0     │
│ A      ┆ 1     │
│ A      ┆ 0     │<-- replaced
│ A      ┆ 0     │<-- replaced
│ B      ┆ 4     │
│ B      ┆ 5     │
│ B      ┆ 0     │<-- replaced
│ B      ┆ 0     │<-- replaced
└────────┴───────┘

我的组织使用 Excel 2016，但计划很快取消使用 VBA 宏或用户定义公式的功能。我们正在尝试寻找适应这一变化的方法。

问题：

我们可以用其他公式替换下面的用户定义公式吗？

用户定义公式：FYFinder

FYFinder 接受对日期的单元格引用（或尝试将数字引用转换为 Excel 日期）并提供该日期的财政年度。新的财政年度从每年的 10 月 1 日开始。因此，对于 2025 年，从 2025 年 1 月 1 日到 2025 年 9 月 30 日的所有日子都属于 2025 财政年度，而从 2025 年 10 月 1 日到 2025 年 12 月 31 日的所有日子都属于 2026 财政年度。日期时间值的时间部分将被忽略，非数字输入会产生错误

例如，给定以下输入，输出为：

• A1:= 11/20/2024 ; FYFinder(A1) = 2025
• A2：= 2017 年 9 月 30 日； FYFinder(A2) = 2017
• A3:= 2017 年 10 月 1 日； FYFinder(A3) = 2018
• A4:= 1 ；FYFinder（A4）= 1900
• A5:= 11/20/2024 5:15 下午 ; FYFinder(A5) = 2025
• A6:= "11/20/2024" ; FYFinder(A6) = #VALUE！
• A7：=“StackOverflow；FYFinder（A7）=#VALUE！

请注意，A4 演示了该函数将所有数字视为日期，而非数字会产生错误

可以使用以下公式轻松复制此功能：

=IF(MONTH(A1)>=10,1+YEAR(A1),YEAR(A1))

或者简洁版本

=--(MONTH(A1)>=10)+YEAR(A1)

那么，鉴于这些信息，我们可以做些什么来在没有 VBA 的情况下复制该功能？我尝试使用自定义命名范围，但我看到的示例无法引用特定单元格。例如：

=--(MONTH(TODAY())>=10)+YEAR(TODAY()) 仅适用于今天。或者=--(MONTH(A1)>=10)+YEAR(A1)仅引用单元格 A1

（我的组织计划很快升级到 Excel 365，因此如果答案来自该版本，那么它可以帮助我们更快地采用它。

假设我们有这样一个用于聚类的数据框架。

# df
dput(df)
structure(c(1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 
0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 
1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 0L, 0L, 0L, 
0L, 1L, 0L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 0L, 
0L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 
1L, 1L, 1L, 1L, 0L, 0L, 0L, 0L, 0L, 0L, 1L, 1L, 1L, 0L, 0L, 0L, 
0L, 0L, 0L, 0L, 1L, 1L), dim = c(9L, 11L), dimnames = list(c("1", 
"2", "3", "4", "5", "6", "7", "8", "9"), c("C", "D", "E", "F", 
"G", "H", "K", "L", "M", "N", "P")))

dist_matrix <- dist(df, method = "manhattan")
clust <- hclust(dist_matrix, method = "complete")
plot(clust)

我们可以得到以下聚类结果，其中斜体字母是附加标记。

请允许我使用伪代码来说明我想要实现的功能。

1. 获取特定节点的左节点和右节点。

FUNC GET_RIGHT_NODE(D)
RETUEN E,F,G
FUNC GET_LEFT_NODE(ROOT)
RETUEN A,B,C

2. 获取特定节点的叶元素

FUNC GET_LEAF(D)
RETURN 5,3,1,2,4

相关主题：

R：帮助分析层次聚类中的聚类内容

如何在 R 中打印 hclust 对象的行？

R 中的层次聚类最近邻算法

我正在尝试将方法转换为表达式树以进行测试，但收到错误

Exception thrown: 'System.ArgumentException' in System.Linq.Expressions.dll Expression of type 'System.Void' cannot be used for return type 'TestDatabase.Models.TestEnum'

同时，我想知道如何获取异常测试输出的输入值。

我尝试重新创建表达式树的原始方法是：

private TEnum ConvertStringToEnum(TDbValueType dbValue)
{
    if (Enum.IsDefined(typeof(TEnum), dbValue))
    {
        return (TEnum)Enum.Parse(typeof(TEnum), (string)(object)dbValue);
    }

    throw new Exception("");
}

表达式树构建方法是

public static Func<TString, TEnum> BuildStringToEnumConvertor<TString, TEnum>()
{
    var enumType = typeof(TEnum);
    var enumTypeExp = Expression.Constant(enumType);

    //var labeleExp = Expression.Label(enumType);

    var dbValueExp = Expression.Parameter(typeof(TString), "dbValue");
    var m = typeof(Enum).GetMethod("IsDefined", [typeof(Type), typeof(object)]);

    //test
    var testExp = Expression.Call(null, m, enumTypeExp, dbValueExp);

    //true
    var objExp = Expression.Convert(dbValueExp, typeof(object));
    var strExp = Expression.Convert(objExp, typeof(string));

    m = typeof(Enum).GetMethod("Parse", [typeof(Type), typeof(string)]);
    var retExp = Expression.Call(null, m, enumTypeExp, strExp);
    objExp = Expression.Convert(retExp, enumType);

    //var returnExp = Expression.Return(labeleExp, objExp);

    //false
    //var x = Expression.Lambda<Func<string>>(dbValueExp).CompileFast();
    //var x = Expression.Lambda<Func<string>>(dbValueExp).Compile();

    var falseExp = Expression.Throw(Expression.Constant(new Exception($"could not convert \" {{how to get the input parameter value ??}} \" to type \"{typeof(TEnum).Name}\".")));

    var xxExp = Expression.IfThenElse(testExp, objExp, falseExp);

    //var block = Expression.Block(objExp, falseExp);

    //return Expression.Lambda<Func<TString, TEnum>>(block, dbValueExp).CompileFast();
    return Expression.Lambda<Func<TString, TEnum>>(xxExp, dbValueExp).Compile(); <-- error here 

}


internal enum TestEnum
{
    Age = 1,
    Name,
    Address,
    CellPhone,
    Max = 500
}

当我打电话时var func = BuildStringToEnumConvertor<string, TestEnum>();，出现以下异常：

Unhandled exception. System.ArgumentException: Expression of type 'System.Void' cannot be used for return type 'Program+TestEnum'
   at System.Linq.Expressions.Expression.ValidateLambdaArgs(Type delegateType, Expression& body, ReadOnlyCollection`1 parameters, String paramName)
   at System.Linq.Expressions.Expression.Lambda[TDelegate](Expression body, String name, Boolean tailCall, IEnumerable`1 parameters)
   at Program.BuildStringToEnumConvertor[TString,TEnum]()

在此处运行演示： https: //dotnetfiddle.net/IkYEpJ

演示程序

#include<stdio.h>
int main(){
static int ia[2];
 static volatile _Atomic (int *) a = (int *) (&ia[1]);
 if ((a += (1)) != (int *) ((int *) (&ia[1]) + (1)))
         printf("1");
}

编译失败

<source>:5:9: error: invalid operands to binary expression ('volatile _Atomic(int *)' and 'int')
    5 |  if ((a += (1)) != (int *) ((int *) (&ia[1]) + (1)))
      |       ~ ^  ~~~
1 error generated.
Compiler returned: 1

Gcc 没问题！

https://godbolt.org/z/fMqzoe4b5

clang 是否支持 _Atomic 指针算法

std::atomic::fetch_add()（或atomic_fetch_add()在 C11 中）基于二进制补码对有符号整数进行运算。这与通常的有符号整数算法不同，后者允许使用其他表示方案。

来自N3337 [atomics.types.operations/30]：

备注：对于有符号整数类型，算术定义为使用二进制补码表示。没有未定义的结果。对于地址类型，结果可能是未定义的地址，但操作没有未定义的行为。

我的问题是：

(1) 在使用除二进制补码以外的有符号整数表示的平台上如何std::atomic::fetch_add()实现？ 是否只有使用互斥锁等的非无锁实现？

(2) 当在一个使用二进制补码以外的表示方法的平台上读取包含负值的原子变量的值时，它是否会自动转换为平台的本机格式？还是用户有责任将其转换为本机表示，因为二进制补码表示的位串是从中读取的？

std::atomic<int> a{0};
a.fetch_sub(1);
int x = a.load();  // The value of a is -1 (0xffffffff), so what is the value of x?

（3）我想知道在采用除二进制补码之外的表示方案的平台上的实现示例std::atomic::fetch_add()（或atomic_fetch_add()）。（如果这样的平台不存在，我不会感到惊讶！）

请注意，N2427对上述内容作出了如下规定：

正常的有符号整数加法运算在发生溢出时具有未定义的行为。对于原子变量，这种未定义的行为会带来很大的负担，因为通常没有办法预先测试溢出。我们认识到现代系统的实际行为，并定义原子提取和添加（减法）以使用二进制补码算法，而不是让溢出保持未定义状态。我们知道没有一个 C++ 实现会因为这个定义而出现问题。

（这个问题是基于学术兴趣。我问这个问题并不是因为我在面对任何实际问题时遇到了困难。）

我正在寻找在二维中平铺 Excel 数组的方法，如另一个数组所定义，理想情况下不使用 VBA。

为了更好地解释，我想要一个 excel 公式f(A, B)

其中 A 是布尔值且 =

[[1, 0],
 [0, 1]

其中 B =

[[1, 2],
 [3, 4]

并f(A, B)根据 A 中的 bool 值对 B 进行平铺，在本例中，f(A, B) =

[[1, 2, 0, 0],
 [3, 4, 0, 0],
 [0, 0, 1, 2],
 [0, 0, 3, 4]]

我最初的想法是=IF(MUNIT(2), {1, 2;3, 4}, {0, 0;0, 0})，但是输出扩展得不太好。