numpy - Fast inverse and transpose matrix in Python -


मेरे पास एक बड़ा मैट्रिक्स ए आकृति का (n, n, 3, 3 ) के साथ n के बारे में 5000 है। अब मुझे मैट्रिक्स a : व्युत्क्रम और स्थानांतरित करना चाहते हैं: numpy as np a = np.random.rand (1000, 1000, 3, 3) पहचान = एनपी। निवासी (3, dtype = A.dtype) Ainv = np.zeros_like (ए) एट्रांस = np.zeros_like (ए) के लिए मैं रेंज (1000) में: सीमा के लिए जी में (1000): ऐन्विन [i, j ] = Np.linalg.solve (ए [i, j], पहचान) एट्रांस [i, j] = np.transpose (ए [i, j])

क्या कोई तेज़ है , यह करने के लिए अधिक कुशल तरीका है?

यह मेरा एक प्रोजेक्ट से लिया जाता है, जहां मैं भी करता हूं Vexized रैखिक बीजगणित कई 3x3 मैट्रिक्स पर।

ध्यान दें कि केवल 3 से अधिक एक लूप है; N पर लूप नहीं है, इसलिए महत्वपूर्ण आयामों में कोड को सदिश किया गया है। मैं यह कैसे C / numba एक्सटेंशन की तुलना में एक ही बात करने के लिए तुलना करने के लिए ज़मानत नहीं करना चाहता, हालांकि, प्रदर्शन के आधार पर। यह काफी तेजी से अभी भी होने की संभावना है, लेकिन कम से कम यह पानी से बाहर n पर छोरों को उड़ाता है। 

  def adjoint (ए): "" "det द्वारा विभाजन के बिना उलटा गणना; ..xv3xc3 इनपुट, या मैट्रिक्स की सरणी के लिए "एआई = एनपी.एम्प्टी_एक्लो (ए)" एक्सरेन्ज (3) में: एआई [..., i ,:] = np.cross (ए [..., I-2 ,:], ए [..., i-1 ,:]) एआई डीएफ़ inverse_transpose (ए): "" "3x3 मैट्रिक्स के ढेर के लिए व्युत्क्रम-पारगमन" कुशलतापूर्वक गणना "" I = adjoint (ए ) Det = डॉट (I, ए)। मीन (अक्ष = -1) रिटर्न I / det [..., कोई नहीं, कोई नहीं] def उलटा (ए): "" "3x3 मैट्रिक्स के एक स्टैक के व्युत्क्रम" "" वापसी Np.swapaxes (inverse_transpose (A), -1, -2) def डॉट (ए, बी): "" "डॉक्स एरेज़ ऑफ वेक्स; अंतिम इंडेक्सेस पर कॉन्ट्रैक्ट" "एनपी.इंसम ('... i। ..i-> ... ..., A, B) A = np.random.rand (2,2,3,3) I = उलटा (ए) प्रिंट np.einsum ('... ij, .. । जेके ', ए, आई)

Comments

Post a Comment

Popular posts from this blog

c - Mpirun hangs when mpi send and recieve is put in a loop -

I was trying to use the program to use mifrewan on 4 node clusters. Distributing node 0 data in nodes 1, 2 and 3. In the program, from '90 'to variable values ​​of' dior 'to be calculated. Therefore the node is distributing data and collecting results in an omitted fashion (for different values ​​of var 'dir') when do {*******} Whereas (dir Loop is given, Mippun is hanging, and no output is found. But when I comment on two {*******} (dir The loop output is obtained for the initial value of variable dir dir = -90 ), and that output is correct. The problem occurs when the loop is given. Can anyone help me in solving this problem. #include "mpi.h" int main (int argc, char * argv []) Float dir = -90; Int rank, nanoprakash; MPI_Status status; MPI_Init (& amp; argc, & argv); MPI_Comm_rank (MPI_COMM_WORLD, & amp; Rank); MPI_Comm_size (MPI_COMM_WORLD, & amp; numprocs); If (rank == 0) {/ / data to begin / / (dest = 1; dest ...
Read more

python - Apply coupon to a customer's subscription based on non-stripe related actions on the site -

I can not find anything in the source about using coupons after subscribing to a user's plan. Is it possible to use a coupon to credit a user's membership after initial subscription? If not, what is the best option? Specifically, I want to give a user 100% coupon for that month, when they take 5 [site specific actions]. This is possible, but may not be the best idea. Coupons are currently applied to the stripe, then apply again at the customer level and each client can have a single coupon. So if you are already offering discounts to the coupon, then the customer who is currently in the present will also give up. It also means that if you are using a recently-added multiple subscription capacity, then all of the user's subscription will be ready for the month. Coupons can be set to "Once", only once (as the only name) can be used by any customer, so if this credit can earn more than once, You will have to maintain a 100% discount coupon continuousl...
Read more

java - Unable to get JDBC connection in Spring application to MySQL -

मैंने स्प्रिंग में डीबी कॉन्फ़िगरेशन का पालन किया है: & lt; bean id = "dataSource" Class = "org.apache.commons.dbcp.BasicDataSource" नष्ट-विधि = "बंद" & gt; & Lt; प्रॉपर्टी नाम = "url" मान = "jdbc: dburl" / & gt; & Lt; संपत्ति नाम = "driverClassName" मान = "com.mysql.jdbc.Driver" / & gt; & Lt; संपत्ति नाम = "उपयोगकर्ता नाम" मान = "उपयोगकर्ता नाम" / & gt; & Lt; संपत्ति नाम = "पासवर्ड" मान = "पासवॉर्फ" / & gt; & Lt; प्रॉपर्टी नाम = "हटाए गए हटाए गए" मान = "सही" / & gt; & Lt; गुण नाम = "प्रारंभिक आकार" मान = "3" / & gt; & Lt; संपत्ति नाम = "अधिकतमअक्टिव" मान = "5" / & gt; & Lt; / सेम ​​& gt; डीबी से गुण मूल्यों को लोड करने के लिए अतिरिक्त कॉन्फ़िगरेशन। & lt; bean id = "config1" class = "org.apache.commons.configuration....
Read more