INNER CODE UNIT · C

js

lh3/yak · 6gjdn.c:9

	js = malloc(n * sizeof(int));
	l = 1;
	for (k = 0; k <= n - 1; k++) {
		d = 0.0;
		for (i = k; i <= n - 1; i++)
			for (j = k; j <= n - 1; j++) {
				t = fabs(a[i * n + j]);
				if (t > d) {
					d = t;
					js[k] = j;
					is = i;
				}
			}
		if (d + 1.0 == 1.0)
			l = 0;
		else {
			if (js[k] != k)
				for (i = 0; i <= n - 1; i++) {

View source record →