小编Yos*_*ian的帖子

如何在python中绘制时间序列

我一直试图从CSV文件中绘制时间序列图.我已经设法读取文件并使用字符串将数据从字符串转换为日期strptime.当我尝试使用包含日期信息的列表在matplotlib中绘制测试图时,它将日期绘制为一系列点; 也就是说,对于2012年5月31日19:00的日期,我在y轴上获得了2012,05,19,31,00处的点的图,其中x = 1的值等等.我知道这不是传递绘图日期信息的正确方法.有人能告诉我如何正确传递这些信息.

python plot matplotlib

43
推荐指数
2
解决办法
14万
查看次数

如何使用magit中的git mv?

有没有一种很好的方法来调用git mvmagit中的文件?我知道可以运行任何git命令:,但这不会自动完成文件名.

git emacs magit

23
推荐指数
3
解决办法
3690
查看次数

将函数应用于3D numpy数组的每个2D切片的有效方法

我想应用一个将2D数组(并返回相同形状之一)的函数应用于3D数组的每个2D切片.这样做的有效方法是什么?numpy.fromiter返回一维数组,numpy.fromfunction需要单独应用于每个坐标.

目前我在做

foo = np.array([func(arg, bar2D) for bar2D in bar3D])
Run Code Online (Sandbox Code Playgroud)

这给了我想要的东西,但列表理解非常慢.此外,func是具有特定边界条件的一维导数.numpy.gradient似乎只有N的数组维度的ND衍生物,但也许有另一个例程可以为我做整件事?

编辑:列表理解有效,但我正在寻找一种更快的方法.bar3D可以很大,最多(500,500,1000).numpy我发现的将函数应用于数组的所有例程似乎都假设函数或数组都是1D.

python arrays numpy

14
推荐指数
1
解决办法
2031
查看次数

就地,mpi_reduce与OpenMPI崩溃

每当我尝试mpi_reduce使用mpi_in_place发送缓冲区进行调用时,它都会崩溃.谷歌的一个拖网显示这对Mac OS OMPI 1.3.3来说是一个问题 - 但我在CentOS上使用OMPI 1.6.3(使用gfortran 4.4.6).

以下程序崩溃:

PROGRAM reduce

  USE mpi

  IMPLICIT NONE

  REAL, DIMENSION(2, 3) :: buffer, gbuffer

  INTEGER :: ierr, me_world
  INTEGER :: buf_shape(2), counts

  CALL mpi_init(ierr)
  CALL mpi_comm_rank(mpi_comm_world, me_world, ierr)

  buffer = 1.
  IF (me_world .EQ. 0) PRINT*, "buffer: ", buffer

  buf_shape = SHAPE(buffer)
  counts = buf_shape(1)*buf_shape(2)

  CALL mpi_reduce(MPI_IN_PLACE, buffer, counts, mpi_real, mpi_sum, 0, mpi_comm_world, ierr)
  IF (me_world .EQ. 0) PRINT*, "buffer: ", buffer

  CALL mpi_finalize(ierr)

END PROGRAM reduce
Run Code Online (Sandbox Code Playgroud)

MPI错误是:

MPI_ERR_ARG: …
Run Code Online (Sandbox Code Playgroud)

fortran mpi

8
推荐指数
1
解决办法
4359
查看次数

TRAMP忽略了tramp-remote-path

我的.emacs中有以下内容:

(require 'tramp)
(add-to-list 'tramp-remote-path "/some/path")
(add-to-list 'tramp-remote-path 'tramp-default-remote-path)
Run Code Online (Sandbox Code Playgroud)

所以当我使用tramp over ssh打开文件时,我希望我PATH能包含/some/path.相反,运行M-! echo $PATH返回

/bin:/usr/bin:/usr/sbin:/usr/local/bin
Run Code Online (Sandbox Code Playgroud)

即使我设置export PATH=/hwdisks/data/modules/pkg/git/1.8.4.1/bin/git:$PATH.bashrc.profile,PATH设置不正确.

在tramp日志中*debug tramp/ssh remotehost*,我可以看到tramp明确设置PATH:

12:28:34.202135 tramp-send-command (6) # PATH=/bin:/usr/bin:/usr/sbin:/usr/local/bin; export PATH
Run Code Online (Sandbox Code Playgroud)

如果我包括echo "in .bashrc"在我的.bashrc,它出现在tramp日志中,所以我知道tramp正在阅读它.

如何让流浪汉使用正确PATH

Emacs版本:24.2.1
Tramp版本:2.2.3-24.1

emacs tramp

8
推荐指数
2
解决办法
1358
查看次数

重命名所有'.' 除了扩展名之外的文件名中的'_'

我正在尝试创建一个替换所有"."的脚本.用"_"出现在文件名中.例如,当我尝试替换所有""符号时,我使用此:

rename 'y/ /_/' '{}' file
# test 1.2.jpg -> test_1.2.jpg
Run Code Online (Sandbox Code Playgroud)

哪个工作正常,但当我尝试使用"."时.符号扩展名也会改变:

rename 'y/./_/' '{}' file
# test 1.2.jpg -> test 1_2_jpg
Run Code Online (Sandbox Code Playgroud)

如何在不更改扩展名的情况下重命名文件(如果有的话)?

regex linux bash rename file-rename

6
推荐指数
1
解决办法
1348
查看次数

停止变量中的代码?

是否可以使用变量中的错误代码/消息来停止或中止fortran程序?这似乎不可能与内在STOP:

integer :: status = 1
character(len=3) :: err_msg = "err"

stop status  
stop err_msg 
Run Code Online (Sandbox Code Playgroud)

两个stop调用都会在编译时抛出语法错误.我错过了什么,或者我是否必须stop 1直接打电话?或者写我自己的包装?

fortran gfortran

5
推荐指数
2
解决办法
711
查看次数

奇异向量 ind matlabsvd 函数的符号确定

有谁知道Matlab函数产生的奇异向量的符号svd是如何确定的?

让:

B = U*S*V'

svd是实数或复数 2×2 矩阵 B 的有效分解,则:

B = (U*c)*S *(V*c)'

也是有效的,其中c是一个改变一个或两个奇异向量符号的矩阵:

c = diag([1 -1]), diag([-1 1]) or diag([-1 -1])

我想知道Matlab的svd算法如何确定U和V中奇异向量的符号。

matlab svd

5
推荐指数
1
解决办法
693
查看次数

在gdb中打印漂亮的Fortran动态类型

在gdb中打印Fortran可分配多态变量的值非常痛苦.鉴于下面的程序,为了看到它的价值alloc_ext,我必须做以下事情:

(gdb) p alloc_ext
$1 = ( _data = 0x606260, _vptr = 0x400ce0 <__foo_MOD___vtab_foo_My_extended_type> )
(gdb) ptype alloc_ext
type = Type __class_foo_My_base_type_a
PTR TO -> (     Type my_base_type :: _data)
PTR TO -> (     Type __vtype_foo_My_base_type :: _vptr)
End Type __class_foo_My_base_type_a
(gdb) ptype alloc_ext%_data
type = PTR TO -> ( Type my_base_type
character*4 :: base_char
End Type my_base_type )
(gdb) p alloc_ext%_data
$2 = (PTR TO -> ( Type my_base_type )) 0x606260
(gdb) p *(my_extended_type*)(alloc_ext%_data)
$3 = ( …
Run Code Online (Sandbox Code Playgroud)

fortran gdb

5
推荐指数
1
解决办法
348
查看次数

np.multiply如何工作?

我试图在Java中实现np.multiply,我对它实际上在做什么很困惑.文档简单地说它是元素乘法.它与我能找到的任何数学矩阵产品都不匹配.它部分匹配元素方面的Hadamard产品,但不需要相同数量的行和列.有谁知道什么数学产品np.multiply执行,并有任何更多的信息到它是如何工作的?

以下是我得到的不同产出.这些似乎是非常不同的功能.

a = np.array([[1,1,1],[2,2,2],[2,3]])
b = np.array([[2,3,4]])
print(np.multiply(a,b))
#output
[[1, 1, 1, 1, 1, 1] [2, 2, 2, 2, 2, 2, 2, 2, 2] [2, 3, 2, 3, 2, 3, 2, 3]]]
Run Code Online (Sandbox Code Playgroud)

a = np.array([[1,1,1],[2,2,2]])
b = np.array([[2,3,4]])
print(np.multiply(a,b))
#output
[[2 3 4]
 [4 6 8]]
Run Code Online (Sandbox Code Playgroud)

python arrays numpy matrix matrix-multiplication

4
推荐指数
1
解决办法
5945
查看次数

如何退出嵌套的 Fortran 循环?

我正在尝试编写一个程序(在 Fortran 95 中),该程序将最多 N 的自然数最小分解为最多 4 个正整数的总和。

一段时间以来,我一直在尝试添加和删除语句,以使其仅停留在最小的分解上,但我一无所获。如何在找到最小分解后立即停止程序?

PROGRAM SummeQuadrat
IMPLICIT NONE
 real:: start,finish
 integer:: a,b,c,d,g,x,y

 write(*,*) "Max n"
 read (*,*) y
 call cpu_time(start)

do x=1,y,1

  do a=0,x,1
   do b=a,x-a,1
    do c=b,x-b,1 
     do d=c,x-c,1
      if (a**2+b**2+c**2+d**2 .eq. x) then 
      write(*,*) "x=",x,d,c,b,a

      end if
     end do 

    end do
   end do
  end do
end do
call cpu_time(finish)
    write(*,*)finish-start
 end program SummeQuadrat
Run Code Online (Sandbox Code Playgroud)

fortran

4
推荐指数
1
解决办法
4391
查看次数

如何在emacs中安装melpa?

我需要将MELPA安装到我的emacs 23.3版本中。我已经搜索了安装它的方法,但是找不到/理解其中的任何一个。

谁能帮我安装MELPA。我正在使用ubuntu 12.04

我需要MELPA安装OmniSharp,这将有助于在C#中自动完成

emacs emacs23 melpa

3
推荐指数
1
解决办法
9124
查看次数

汇编 x86 EDI、SIL 和 DL 寄存器的功能是什么?

lk:
movl %edi , %r9d
jmp .L8

L8:
movzbl %sil , %eax
movzbl %dl , %ecx
leal -1(%rcx) , %r8d
cmpl %r8d, %eax
jl .L4 //if (%eax <%r8d) goto .L4
cltq // extend %eax int %rax
movl %esi, %eax
cmovle %edx, %eax // if (ecran(%rax)<=%dil) %eax=%edx
ret

main:
...
movzbl %bl,%eax
movsbl clavier(%rax), %ebp
movsbl le(%rip),%edx
movl $0, %esi
movl %ebl, %edi
call lk
movzbl %al, %esi
...
ret
Run Code Online (Sandbox Code Playgroud)

什么%edi%sil%dl代表在功能lk和他们有什么类型?

x86 assembly

-4
推荐指数
1
解决办法
8934
查看次数