c - Unable to find determinant using partial pivoting -


whats problem. executes doesnt show correct answer. gets compiled well. i've used calling function in many places. assumed matrix square , give input througgh terminal. ex=3 random 3x3 matrix value seems incorrect

#include<stdio.h> #include<math.h> #include<stdlib.h>  void row(int r1,int r2,int n,float a[n][n]){ int c; float temp[n][n]; for(c=0;c<n;c++){ temp[r1][c]=a[r1][c]; a[r1][c]=a[r2][c]; a[r2][c]=temp[r1][c]; } }  void maximum(int i, int n,float a[n][n]){ int j; float max=fabs(a[i][i]); for(j=i+1;j<n;j++){     if(fabs(a[j][i])>max){         max=fabs(a[j][i]);         row(i,j,n,a);     } } }  void op(int k, int n,float a[n][n]){ int i,j; float f;f for(i=1;i<n-1-k;i++){      f=-(a[k+i][k]/a[k][k]);     for(j=0;j<n-1-k;j++){            a[k+i][k+j]=a[k+i][k+j]+f*(a[k][k+j]);     } } }     int main(){ int i,j,n; printf("enter order of matrix:"); scanf("%d",&n); float a[n][n]; for(i=0;i<n;i++){     for(j=0;j<n;j++)     scanf("%f",&a[i][j]);     }  float det=1; for(i=0;i<n-1;i++){     maximum(i,n,a);     op(i,n,a);         det*=a[i][i];     } det*=a[n-1][n-1]; printf("%f\n",det); return 0; }

you reading uninitialized values

for(j=i+1;j<n;j++){     if(fabs(a[j][i])>max){         max=fabs(a[j][i]);         row(i,j,n,a);     } }

should be

for(j=i+1;j<n;j++){     if(fabs(a[i][j])>max){         max=fabs(a[i][j]);         row(i,j,n,a);     } }

Comments

Post a Comment

Popular posts from this blog

python Tkinter Capturing keyboard events save as one single string -

is possible store keyboard events 1 string? code below stores 1 char , prints it. card reader or bar code reader, contains collection of character/ string not 1 character @ time. goal save char pressed text variable. from tkinter import * root = tk() def key(event): text= event.char text+= event.char print ("pressed", text) def callback(event): frame.focus_set() print ("clicked at", event.x, event.y) frame = frame(root, width=100, height=100) frame.bind("<key>", key) frame.bind("<button-1>", callback) frame.pack() root.mainloop() currently, creating text variable , everytime key() function gets called, , text stores last character typed. you can define text module level variable , , use module level text inside key function - text = '' def key(event): global text text+= event.char print("pressed", text)
Read more

android - InAppBilling registering BroadcastReceiver in AndroidManifest -

i'm implementing inapp billing in android app. works well, however, i'm trying decouple broadcast receiver activity manifest. suggestion in android's trivialdrive sample: // important: dynamically register broadcast messages updated purchases. // register receiver here instead of <receiver> in manifest // because call getpurchases() @ startup, therefore can ignore // broadcasts sent while app isn't running. // note: registering listener in activity bad idea, done here // because sample. regardless, receiver must registered after // iabhelper setup, before first call getpurchases(). currently there's class extends broadcastreceiver : public class iabbroadcastreceiver extends broadcastreceiver { /** * intent action receiver should filter for. */ public static final string action = "com.android.vending.billing.purchases_updated"; private final iabbroadcastlistener mlistener; public iabbroadcastreceiver(iabbroadcastlistener listener) { m
Read more

javascript - Z-index in d3.js -

i have map , when hover effect on region of map below yellow line , circle, last ones disappear. how can solve problem? http://jsfiddle.net/ot21n9qx/ var arc = g.append("path") .style("fill", "none") .style("stroke", "yellow") .style("stroke-width", 2) .attr("d", "m" + pathorigin[0] + "," + pathorigin[1] + " q" + svgpoint[0] + "," + pathorigin[1] + " " + svgpoint[0] + "," + svgpoint[1]); var circlesize = d3.scale.linear() .domain([0, 0.5, 1]) .range([4, 10, 4]); svg not have z order draws objects in order created. need create path , circle after create highlight regions. here forked fiddle: http://jsfiddle.net/o00fzgaf/1/
Read more