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

wpf - ControlTemplate with changes for every class -

I have the following ControlTemplate that is the same for all my custom controls. & amp; the lift; ControlTemplate x: key = "PssFunctionControlBaseHorizontal" TargetType = "Local: PssFunctionControlBase" & gt; & Lt; Grid & gt; & Lt; Threshold Threshold = "1" & gt; & Lt; Content Control X: Name = "Inner Content" Template = "{Static Resources Inner Contentbase Horizontal}" /> & Lt; / Border & gt; & Lt; Control x: name = "PART_ResizeDecorator" visibility = "short" template = "{static resource resize deccurratematet}" /> & Lt; / Grid & gt; & Lt; ControlTemplate.Triggers & gt; & Lt; Data Trigger Value = "True" Binding = "{Binding Editing Moded, Relative Soros = {ResolveSource Search Engineer, Instant Type = {x: Type Local: PssViewLayoutControl}}}" & gt; & Lt; Setter target name = "PART_ResizeDecorator" ...
Read more

sonata admin - how to add ManyToOne property field in configureListfield in sonataAdmin -

I have a phone unit unit that contains many properties for many unit entity classes I An exception is given during the rendering of a template ("You must define a associated_property option or Create a \ __counter __ \ Examens \ ExamensBundle \ unit \ Grade Proxy :: __ toString service examens.examens.admin.fonctionnaire from the field option idGrade to the method is in SonataDoctrineORMAdminBundle "): CRUD: Line_ on line_list_orm_many_to_one.html.tw Fontionnaire.php class: & lt ;? Php namespace exams \ ExamensBundle \ Entity; Use the theory / ORM \ mapping as an ORM; / ** * PhoneNationLayer * / Class Fictionalair {//...//// ** * @VARAEfimens \ Examinations Bundle \ Unit \ Grade * / Private $ IDGed; //...//} FonctionnaireAdmin.php: & lt ;? Php namespace exams ExamensBundle \ Admin; Use Sonata \ AdminBundle \ Admin \ Admin; Sonata \ AdminBundle \ Datagrid \ ListMapper; Sonata \ AdminBundle \ Datagrid \ DatagridMapper; Sonata \ AdminBundle \ Validat...
Read more

sql - Return Function using Cursor -

My current schedule is: USERID Department NAME Real Estimated 522 Xyz Department John Wayne 2/1/14/14 8:34 522 Investigation Department John Wayne 2/17/1411: 21 78 9Except Department Bill Smith 2/17/14 9:41 78 9 Excise Department Bill Smith 2/17/14 11:16 78 9 Excise Department Bill Smith 2 / 17/14 11:50 78 9 Excise Department Bill Smith 2/17/14 12: 18 856 Explore Department Raul Castill 2/17/14 9:25 I have to create a function which enables me to enter a user ID and that Allows the user to return data from the person. However, I want to take the first line and want to hold the actual date and then put it in the 'column' check and I also carry the last row and store it in 'Check Out' in the column call. Basically, for every day, I want to record the records in the user's time and want to take the time for each day. This is what I have done so far " return @RtnValue table (useridd integer, deptname varchar (max), emplname varchar (100), actualDate da...
Read more