Performa indeks di MATLAB bagian 2

Pada kesempatan kali ini sekali lagi saya menekankan pentingnya menghindari looping dalam pemrograman MATLAB. Kebetulan ini merupakan soal yang saya berikan dalam pelatihan MATLAB kemarin kepada teman-teman. Saya g tau apa ini bisa dikerjakan atau tidak, yang jelas saya kasikan aja jawabannya, sekalian berbagi pengetahuan. Semoga bisa diambil manfaatnya.

Soalnya adalah diberikan sebuah matriks yang ukurannya 1000 x 1000. Kenapa 1000 x 1000, karena biar kita terbiasa dengan perhitungan sesungguhnya di lapangan. Pada matriks ini,  sederhana saja yang dilakukan yakni, jika penjumlahan baris dan kolom genap, isi matriks tersebut dengan angka dua, sementara jika ganjil, isi dengan 1. Adapaun hasil codingnya adalah sebagai berikut.

	% nomor 1
	function nomor1
	% draw a matrix 1000 x 1000
	% if row + col is even, set to 2, else set to 1
	clear all ;
	clc;
	dim = 1000;
	% first way, iterate trough looping
	Y = zeros(dim,dim) ;
	tic ;
	for i=1:length(Y),
	for j=1:length(Y(1,:)),
	if mod(i+j, 2) == 0
	Y(i,j) = 2;
	else
	Y(i,j) = 1;
	end
	end
	end
	disp(['lama cara 1 = ' num2str(toc)]);
	% disp(Y);
	% second way, using MATLAB index
	X = zeros(dim,dim) ;
	tic;
	[a,b] = ind2sub(size(X), 1:numel(X)) ;
	indeks_genap = mod(a+b,2) == 0;
	indeks_ganjil = mod(a+b,2) > 0;
	XX = sub2ind(size(X), a(indeks_genap),b(indeks_genap));
	YY = sub2ind(size(X), a(indeks_ganjil),b(indeks_ganjil));
	X(YY) = 1;
	X(XX) = 2;
	disp(['lama cara 2 = ' num2str(toc)]);
	% disp(X);
	disp('Is the first get the same output with the second?')
	indeks = find(Y(:,:) ~= X(:,:));
	if(numel(indeks) == 0),
	disp('oh... yes.. ');
	end
	% contoh hasil eksekusi
	% lama cara 1 = 2.72
	% lama cara 2 = 0.2636
	% Is the first get the same output with the second?
	% oh... yes..

view raw nomor1.m hosted with ❤ by GitHub

Jadi terlihat jelas bahwa performa indeks 10 kali lipat lebih cepat ketimbang pake looping.