Jared Frazier's questions

问题正是标题：Fortran 允许在SUBROUTINE即使有的情况下更改派生类型数据成员INTENT(IN)，如何保护数据？我对 Fortran 还很陌生，尽管INTENT(IN)派生类型伪参数具有属性，但这种允许改变派生类型数据成员的行为让我非常担心。有没有办法用编译器捕获这些不安全的突变？例如，gfortran -Wall似乎没有这样做，但也许这是遗留甚至现代代码库中一个相当常见的问题，以至于可以通过某种方式对问题进行某种静态检测？

示例程序：

! @file main.f90
!
! @purpose: Show data members of derived types are mutable in SUBROUTINE 
!
! @compile: gfortran -Wall main.f90
PROGRAM test
    IMPLICIT NONE

    INTEGER :: n_states = 1
    INTEGER :: len_state = 5 
    INTEGER :: i, j

    TYPE t_states 
        INTEGER, POINTER, CONTIGUOUS :: state(:) ! (len_state,)
    END TYPE t_states

    TYPE(t_states), ALLOCATABLE :: states(:) ! (n_states,)

    ALLOCATE(states(1:n_states)) 

    ! Populate allocated array of states
    DO i = 1, n_states 
        ALLOCATE(states(i)%state(len_state))
        DO j = 1, len_state 
            states(i)%state(j) = j
            PRINT *, states(i)%state(j)
        END DO
        PRINT *
    END DO
 
    CALL unsafe_mutation(states(1))

    ! Print state array after modification
    DO j = 1, len_state
        PRINT *, states(1)%state(j)
    END DO 

CONTAINS

    ! Arbitrary update of state data member even though INTENT(IN)!!!
    SUBROUTINE unsafe_mutation(p_states)
        TYPE(t_states), INTENT(in) :: p_states
        p_states%state(1) = -1 
    END SUBROUTINE

END PROGRAM test

输出：

           1
           2
           3
           4
           5

          -1
           2
           3
           4
           5

我正在尝试将DRY原则应用于玩具绘图类，作为一种智力练习，以提高我对 OOP 的理解（目前正在阅读Python 面向对象编程），但直观地看，使用潜在的越来越深的继承层次结构可能会导致更复杂的代码库（特别是因为我已经读过组合似乎比继承更受青睐，请参阅组合优于继承 wiki）。似乎这样的玩具库最终可能会有太多的抽象类，如AbstractMonthlyMultiPanelPlot等等AbstractSeasonalPlot，用于任意绘图类型以适应不同的输入数据。

有没有一种更符合 Python 风格的方法来处理下面我可能遗漏的问题？我是否违反了某种我误解或完全忽略的设计原则？

from abc import abstractmethod, ABC
from numpy import ndarray
from typing import List, Tuple
import matplotlib.pyplot as plt 

class AbstractPlot(ABC)
    @abstractmethod
    def plot(self): 
        raise NotImplementedError

class AbstractMonthlyPlot(AbstractPlot):
    @abstractmethod
    def plot_for_month(ax, data):
        raise NotImplementedError

    @property
    def n_months(self): 
        """number of months in a year"""
        return 12 

    def plot(self, month_to_data: List[Tuple[ndarray]]):
        fig, axs = plt.subplots(self.n_months, 1)
        for month in range(self.n_months):
            self._plot_for_month(ax=axs[month], data=month_to_data[month])

class Contour(AbstractMonthlyPlot):
    def _plot_for_month(self, ax, data):
        ax.contourf(*data)

class Linear(AbstractMonthlyPlot):
    def _plot_for_month(self, ax, data):
        ax.plot(*data)

我正在使用shdoc为我的 bash 函数生成 markdown 文件；但是，我还想将相关的 bash 函数主体以及该主体中的注释附加到由 生成的 markdown 文件中shdoc。如何打印 bash 函数主体，同时还包括函数主体中的注释？

举一个简单的例子，假设我的所有 bash 文件都包含在一个名为的目录中src，并且所有这些文件都以 结尾.sh。文件中的以下函数src/foo.sh

# @description A function that does something
foo() {
  # oh wow, a comment
  echo "something cool"
}

现在，如果我想获取没有注释的函数主体，我可以轻松地做到这一点，declare -f如下所示，在名为的文件中src/doc.sh：

# assuming you are in the `src` dir and both .sh files are executable
. foo.sh
foo_doc=$(declare -f foo)
echo "$foo_doc"

但问题是foo_doc只有函数体，因为echo "$foo_doc"

foo () 
{ 
    echo "something cool"
}

并且没有任何注释。因此，如果我将此信息附加到由 生成的 markdown 中shdoc，它就不是源代码的理想文档。

我认为一种方法是使用sed来确定与感兴趣的函数相对应的行（我的 bash 函数由 和 “分隔” {，}但根据评论我理解还有其他方法可以声明有效的 bash 函数），然后我可以使用来仅获取这里sed完成的那些行。

我了解如何分别使用 1、2 和 3for循环来定义 1、2 和 3 维的求和，但是如何创建一个对多索引求和的通用函数（如下面定义的函数）？

为简单起见，假设u仅采样一个随机整数，即

from random import randint
def u(a=0, b=100):
    return randint(a, b)

举个例子，假设多索引向量js = [0, 0, 0]和上限向量Ms = [5, 5, 5]，这意味着d = 3（即3维多索引），并且求和可以简单地写成如下

mysum = 0
for j1 in range(Ms[0]):
    for j2 in range(Ms[1]):
        for j3 in range(Ms[2]):
            mysum += u()

但是编写d尺寸的一般方法是什么（如图所示）？

编辑：我知道在这个简单的例子中，人们可以简单地做类似的事情

from numpy import prod
M = prod(Ms)
mysum = sum([u() for i in range(M)]) 

但同样，我对一种执行d嵌套for循环的通用方法特别感兴趣，因为该函数u可能更复杂并且取决于索引 j_1、j_2、...、j_d。