matlab - Vectorizing in Multiple Dimensions -

-

i have following loop in matlab code i'd rid of:

for i=1:size(thepolygon,3)                 activevalues(:,i) = sum(normalvectors.*thepolygon(:,:,i),2);     checkvalues(:,i) = sign(activevalues(:,i)-sum(normalvectors.*thepolygon3(:,:,i),2)); end

i tried replacing i 1:size(thepolygon,3), dimensions don't line element-by-element multiplication, , i'm not sure else try. in advance tips.

it's difficult test without example data but:

activevalues = squeeze(sum(bsxfun(@times, normalvectors, thepolygon), 2))

ok first test data:

a = rand(3,3,3); b = rand(3,3);

now test yours

for i=1:size(a,3)     c(:,i) = sum(b.*a(:,:,i),2); end

on run got:

c =      0.9773    1.0608    0.3673     0.6670    0.1597    0.7296     0.8372    1.1418    0.9828

and mine:

squeeze(sum(bsxfun(@times, b, a), 2))

i get

ans =      0.9773    1.0608    0.3673     0.6670    0.1597    0.7296     0.8372    1.1418    0.9828

so assuming mine correct.

explanation:

bsxfun(@times, b, a) broadcasting, expands b along it's singleton dimension (in case dim 3) match size of a , applies function (@times .*) element element expanded b , a. identical (but better practice , faster) going repmat(b, [1 1 size(a,3)]).*a.

then sum 3d matrix column wize have, i.e. sum(x, 2) returns 1x3x3 result. want 3x3 result rid of singelton dimension (i.e. dim 1 equals 1) used squeeze

it should trivial vectorize second line in same way.

but

tic; k = 1:10000 c = squeeze(sum(bsxfun(@times, b, a), 2)); end; toc elapsed time 0.394828 seconds.  tic; c = zeros(3,3); k = 1:10000 i=1:size(a,3) c(:,i) = sum(b.*a(:,:,i),2); end; end; toc elapsed time 0.113199 seconds.

the non-vectorized code faster in matlasb 2012b!!!!!!!


Comments

Post a Comment

Popular posts from this blog

java - Jmockit String final length method mocking Issue -

-
using jmockit 1.2 jar, trying mock string's length method getting unexpected invocation exception: failed: test java.lang.illegalstateexception: missing invocation mocked type @ point; please make sure such invocations appear after declaration of suitable mock field or parameter @ stringdemo.testa$1.<init>(testa.java:17) @ stringdemo.testa.test(testa.java:13) i using testng: @test public void test() throws exception { new expectations() { @mocked("length") string astring; { astring.length(); result = 2; } }; system.out.println(a.showa("test")); } } actual class a: public class { public static int showa(string str){ int a= str.length(); return a; } } this wrong way of recording expected results. shouldn't mock , record string's length() method, record showa() instead. here solution @suppresswarnings("unused...
Read more

What is the difference between data design and data model(ERD) -

-
i have google question , find answer in day. know, data design many many table have list column,and mark pk , fk. why er-diagram have find this? data design process of designing database. main output of data design detailed logical data model of database. some people data design includes of needed logical , physical design choices , physical storage parameters needed generate design in data definition language. technically, database analysis, database analysts (dbas) trained do. while person can both data design , database analysis, these 2 different tasks. data design models data. database analysis takes model , applies 1 or more database engines (relational, hierarchical, nosql). a logical data model 1 of main outputs of data design. data model represented entity relationship diagram, or er diagram.
Read more

ios - Can NSManagedObject conform to NSCoding -

-
i need transfer single object across device. right converting nsmanagedobject dictionary , archiving , sending nsdata. upon receiving unarchiving it. transfer nsmanagedobject archiving , unarchiving instead of creating intermediate data object. @interface test : nsmanagedobject<nscoding> @property (nonatomic, retain) nsstring * title; @end @implementation test @dynamic title; - (id)initwithcoder:(nscoder *)coder { self = [super init]; if (self) { self.title = [coder decodeobjectforkey:@"title"]; //<crash } return self; } - (void)encodewithcoder:(nscoder *)coder { [coder encodeobject:self.title forkey:@"title"]; } @end nsdata *archivedobjects = [nskeyedarchiver archiveddatawithrootobject:testobj]; nsdata *objectsdata = archivedobjects; if ([objectsdata length] > 0) { nsarray *objects = [nskeyedunarchiver unarchiveobjectwithdata:objectsdata]; } the problem above code is. crashes @ self.title in initwithcoder sa...
Read more