vectorization - Accelerating Iterations- MATLAB -

-

consider 2 vectors a = [20000000 x 1] , b = [20000000 x 1 ]

i need find sum of corresponding every unique element of b.

although looks easy, taking forever in matlab.

currently, using 

u = unique(b); length_u = length(u); c = zeros(length_u,1);  = 1:length_u    c(i,1) = sum(a(b==u(i))); end

is there anyway make run faster? tried splitting loop , running 2 parfor loops using parallel computing toolbox(because have 2 cores). still takes hours.

p.s: yes, should better computer.

you must see this answer first.
if must, can use combination of histc , accumarray

a = randi( 500, 1, 100000 ); b = randi( 500, 1, 100000 );  ub = unique( b );  [ignore idx] = histc( b, [ub-.5 ub(end)+.5] ); c = accumarray( idx', a' )';

see toy comparison naive for-loop implementation on ideone.

how work?

we use second outout of histc map elements of b (and later a) bins defined elements of ub (the unique elements of b).
accumarray used sum entries of a accorind mapping defined idx.
note: assume unique elements of b @ least 0.5 apart.


Comments

Post a Comment

Popular posts from this blog

sql server - MSSQL Text and Varchar(MAX) fields shown (MEMO) in DBGrid -

-
Image
i have encountered problem in mssql & delphi xe3. text , varchar(max) fields seemed (memo) or (widememo) in tdbgrid & tdbadvgrid (tms dbgrid). in first screenshot can see table fiels, in second screentshot can see sample data in database: in screenshot, see how in dbgrid or in query results. i figured out changing field type char , varchar not max solves issue don't know why having problem. can please let me know?
Read more

php - Changing the visibility scope of parent methods in child classes -

-
i've got validator class , uservalidator class extends it. my validator has public method setrule(...) public visibility. when extend want change visibility of setrule(...) parent method private/protected within child it's visible child , no outsiders can call method from child. is possible? if so, how achieve it? from architectural point of view not recommended. stated in comments clean way set method protected children can access it. i cannot think of single use case put me in need call public method on parent class not allowed call on child class. that's against open/closed principle. classes should open extension, not modification. since not question i'll provide way how can achieved though. note: this method makes use of class responsible instantiation it's hack. solution not make use of php's native language features when throwing accessibility errors. first let's define classes had <?php class validator { ...
Read more

qml - Is it possible to implement SystemTrayIcon functionality in Qt Quick application -

-
i writing qtquick desktop application. use both c++ (for functionality) , qml (for ui) in it. use qquickview show interface written in qml. i want application reside in system tray when minimised. i mean functionality similar example. http://qt-project.org/doc/qt-4.8/desktop-systray.html . trying implement feature not find way in qt quick application. here main.cpp code: #include <qguiapplication> #include <qqmlengine> #include <qqmlcontext> #include <qqmlfileselector> #include <qquickview> #include "myapp.h" int main(int argc, char* argv[]) { qguiapplication app(argc,argv); app.setapplicationname(qfileinfo(app.applicationfilepath()).basename()); qdir::setcurrent(qapp->applicationdirpath()); myapp myappobject; qquickview view; view.connect(view.engine(), signal(quit()), &app, slot(quit())); view.rootcontext()->setcontextproperty("myappobject", &myappobject); new qqmlfilese...
Read more