From 2 Numpy vectors(length=m,n) , create crossed matrix(shape=(m, n,2))

As a reminder for MYSELF,

To create a matrix from two 1D vectors, like 

v1 = np.arange(4)           #[0,1,2,3]
v2 = np.arange(4) *10    #[0,10,20,30]

to create crossed matrix 

[[[0,0],[0,10],[0,20],[0,30]],
  [[1,0],[1,10],[1,20],[1,30]],
  [[2,0],[2,10],[2,20],[2,30]],
  [[3,0],[3,10],[3,20],[3,30]]]
I do 
crossed = np.dstack(( np.zeros((len(v1),len(v2)), dtype=v1.dtype) + v1[...,np.newaxis], 
                                    np.zeros((len(v1),len(v2)),dtype=v2.dtype) + v2))
Is there any better way to do it? Any suggestion will be appreciated.

Comments

Post a Comment

Popular posts from this blog

Calling OpenCV functions via Cython from Python 3.X.

I have managed to run OpenCV function from Python3 via Cython today, so I want to write about it. Before you start something like this, you should know you have to read at least header files of C library to write right definition in Cython,because Cython does not do that for you. I am using virtualenv, python from pythonz, and home-brewed OpenCV. Almost everything of C++ works just fine with Cython, but still, you should do some Voodoo things like for integer template argument. ctypedef void* int_parameter ctypedef int_parameter two "2" ctypedef Point_[float, two] Point2f Now, main sample Cython code is following. #===== testopencv.pyx ======== import numpy as np cimport numpy as np # for np.ndarray from libcpp.string cimport string from libc.string cimport memcpy cdef extern from "opencv2/core/core.hpp": cdef int CV_WINDOW_AUTOSIZE cdef int CV_8UC3 cdef extern from "opencv2/core/core.hpp" namespace "cv": cdef cppclass Mat: Ma...
Read more

How to update lookup fileld of Kintone using pykintone

 When you try to update kintone record's lookup field, you will fail see the error saying you don't have enough privilege to do that. "フィールド「XXXX」の値「YYYY」が、ルックアップの参照先のフィールドにないか、またはアプリやフィールドの閲覧権限がありません。" This is because you need to have both writing access to looking-up application and reading access to looked-up application. You can include two api tokens in header for REST API, you can also do that by pykintone just connect two api tokens by comma. import pykintone app = pykintone.app('tr3xr', "30", "tokenXXX...XXX,tokenYYY....YYY") Isn't that easy?
Read more

Showing CPU/Memory usage on tmux status bar(tmuxのステータスバーにCPUとMemoryの使用状況を表示する)

Image
On OSX, to show CPU usage, script is like this. #!/bin/bash top -l 1 | head -n 4 | tail -n 1|awk '{printf "%3d%%\n", $3}' to get Memory Usage, script is like this. #!/bin/bash vm_stat | awk 'BEGIN{FS="[:]+"}{if(NR<7&&NR>1)sum+=$2; if(NR==2||NR==4||NR==5)free+=$2} END{printf "%3d%%\n",100*((sum - free)/sum)}' Name thsese scripts like getCpuUsage.sh/getMemUsage.sh and put them in folder where the $PATH is set. Change their mode to u+x. To show them on tmux status line, .tmux.conf is like set -g status-right '#[fg=blue,bold]#H#[default] #[fg=blue,bold][CPU=#(getCpuUsage.sh) MEM=#(getMemUsage.sh)]#[default]'
Read more