[Numpy-discussion] another seemingly unNumpyizeable routine

Rob europax at home.com
Sat Sep 1 11:58:08 EDT 2001

I hate to bother you guys, but this is how I learn.  I have had some
success using the take() and add.reduce() routines all by myself. 
However this one seems intractable:

for IEdge in range(0,a.HybrdBoundEdgeNum):  
    for JEdge in range(0,a.HybrdBoundEdgeNum):  
        AVec[IEdge+a.TotInnerEdgeNum] +=

## below is my attempt at this, but it doesn't
## h=a.TotInnerEdgeNum
## for IEdge in range(0,a.HybrdBoundEdgeNum):
##     AVec[IEdge+h] = add.reduce(a.Cdd[IEdge,0:a.HybrdBoundEdgeNum] *
I think the problem lies in that add.reduce is using the slice
0:a.HybridBoundEdgeNum for Cdd, but then has to contend with XVec slice
which is offset by h.  Is there any elegant way to deal with this.  This
routine is part of a sparse matrix multiply routine.  If I could find
SparsePy, maybe I could eliminate it altogether.

Thanks,  Rob.

The Numeric Python EM Project


More information about the NumPy-Discussion mailing list