added the given libraries
parent
31843b8448
commit
223de5d8fd
|
|
@ -0,0 +1,41 @@
|
|||
/*
|
||||
'lrmp.c' Obfuscated by COBF (Version 1.06 2006-01-07 by BB) at Thu Nov 24 09:52:47 2011
|
||||
*/
|
||||
#include<stdio.h>
|
||||
#include<malloc.h>
|
||||
#include"lrmp.h"
|
||||
#define lt static
|
||||
#define lh double
|
||||
#define lg int
|
||||
#define lv LRmP
|
||||
#define lf for
|
||||
#define lm fabs
|
||||
#define lo if
|
||||
#define ln printf
|
||||
#define ll return
|
||||
#define ly Solve
|
||||
#define l0 malloc
|
||||
#define l1 sizeof
|
||||
#define ls vwsubs
|
||||
#define lw rwsubs
|
||||
#define lz free
|
||||
#define lx else
|
||||
#define l2 LRoP
|
||||
#define lu void
|
||||
lt lh lp=1.0e-8;lg lv(lg ld,lh* *lc,lg*lk){lg la,lb,le;lh lj;lf(la=0;
|
||||
la<ld;la++)lk[la]=la;lf(lb=0;lb<ld-1;lb++){lj=lm(lc[lb][lb]);le=lb;lf
|
||||
(la=lb+1;la<ld;la++){lo(lj<lm(lc[la][lb])){le=la;lj=lm(lc[la][lb]);}}
|
||||
la=lk[le];lk[le]=lk[lb];lk[lb]=la;lf(la=0;la<ld;la++){lj=lc[le][la];
|
||||
lc[le][la]=lc[lb][la];lc[lb][la]=lj;}lo(lm(lc[lb][lb])<lp){ln("\n\x20"
|
||||
"\x44\x61\x73\x20\x44\x69\x61\x67\x6f\x6e\x61\x6c\x65\x6c\x65\x6d\x65"
|
||||
"\x6e\x74\x20\x69\x73\x74\x20\x6d\x69\x74\x20\x25\x32\x30\x2e\x31\x36"
|
||||
"\x65\x20\x7a\x75\x20\x6b\x6c\x65\x69\x6e\x2e",lc[lb][lb]);ln("\n\x20"
|
||||
"\x64\x65\x72\x20\x25\x33\x64\x2e\x20\x53\x63\x68\x72\x69\x74\x74\x20"
|
||||
"\x64\x65\x72\x20\x4c\x52\x2d\x5a\x65\x72\x6c\x65\x67\x75\x6e\x67\x20"
|
||||
"\x6b\x61\x6e\x6e\x20\x6e\x69\x63\x68\x74\x20\x61\x75\x73\x67\x65\x66"
|
||||
"\x75\x65\x68\x72\x74\x20\x77\x65\x72\x64\x65\x6e",lb+1);ll 1;}lf(la=
|
||||
lb+1;la<ld;la++){lc[la][lb]=lc[la][lb]/lc[lb][lb];lf(le=lb+1;le<ld;le
|
||||
++)lc[la][le]=lc[la][le]-lc[la][lb] *lc[lb][le];}}ll 0;}lg ly(lg ld,
|
||||
lh* *lc,lh*li,lg*lk){lg la;lh*lj;lo(!(lv(ld,lc,lk))){lj=(lh* )l0(ld*
|
||||
l1(lh));lf(la=0;la<ld;la++)lj[la]=li[lk[la]];ls(ld,lc,lj);lw(ld,lc,lj
|
||||
);lf(la=0;la<ld;la++)li[la]=lj[la];lz(lj);ll 0;}lx ll 1;}
|
||||
|
|
@ -0,0 +1,42 @@
|
|||
#ifndef _LRMP_H_
|
||||
#define _LRMP_H_ 1
|
||||
|
||||
#include "vorrueckwaertsub.h"
|
||||
#include <math.h>
|
||||
|
||||
int LRmP (int n, double ** a, int * p);
|
||||
/* LR-Zerlegung mit halbmaximaler Pivotsuche */
|
||||
/* Eingabe: */
|
||||
/* n : Groesse der Matrix, */
|
||||
/* a : Matrix, */
|
||||
/* p : Vektor der Laenge n fuer Permutationen. */
|
||||
/* Rueckgabe: */
|
||||
/* a : LR-zerlegte Matrix, */
|
||||
/* unterhalb der Diagonale ist L gespeichert, */
|
||||
/* in der Diagonale und oberhalb R, */
|
||||
/* p : Permutationen bei Pivotsuche. */
|
||||
/* Rueckgabewert: */
|
||||
/* bei Fehler */
|
||||
/* (Matrix 'fast' singulaer, d.h. */
|
||||
/* Pivotelement kleiner als eps = 1.0e-8) */
|
||||
/* wird 1 zurueckgegeben, */
|
||||
/* sonst 0. */
|
||||
|
||||
int Solve (int n, double ** a, double * b, int * p);
|
||||
/* Berechnet Loesung von Ax=b mittels : */
|
||||
/* Routine LRmP(n,a,p) fuer LR-Zerlegung, */
|
||||
/* dann Substitutionen mit 'vorrueckwaertsub'. */
|
||||
/* Eingabe: */
|
||||
/* n : Groesse des Systems, */
|
||||
/* a : Matrix, */
|
||||
/* b : rechte Seite des LGS, */
|
||||
/* p : Vektor, um die Permutationen zu speichern. */
|
||||
/* Rueckgabe: */
|
||||
/* b : Loesung, */
|
||||
/* p : Permutationen von LRmP, */
|
||||
/* a : LR-zerlegte Matrix. */
|
||||
/* Rueckgabewert: */
|
||||
/* falls die LR-Zerlegung nicht funktioniert 1, */
|
||||
/* sonst 0. */
|
||||
|
||||
#endif
|
||||
|
|
@ -0,0 +1,26 @@
|
|||
/*
|
||||
'vorrueckwaertsub.c' Obfuscated by COBF (Version 1.06 2006-01-07 by BB) at Thu Nov 24 09:52:47 2011
|
||||
*/
|
||||
#include"vorrueckwaertsub.h"
|
||||
#define lt static
|
||||
#define lh double
|
||||
#define lg int
|
||||
#define lv LRmP
|
||||
#define lf for
|
||||
#define lm fabs
|
||||
#define lo if
|
||||
#define ln printf
|
||||
#define ll return
|
||||
#define ly Solve
|
||||
#define l0 malloc
|
||||
#define l1 sizeof
|
||||
#define ls vwsubs
|
||||
#define lw rwsubs
|
||||
#define lz free
|
||||
#define lx else
|
||||
#define l2 LRoP
|
||||
#define lu void
|
||||
lu ls(lg ld,lh* *lr,lh*li){lg la,lb;lf(la=0;la<ld;la++){lf(lb=0;lb<la
|
||||
;lb++)li[la]=li[la]-lr[la][lb] *li[lb];}ll;}lu lw(lg ld,lh* *lq,lh*li
|
||||
){lg la,lb;lf(la=ld-1;la>=0;la--){lf(lb=la+1;lb<ld;lb++)li[la]=li[la]
|
||||
-lq[la][lb] *li[lb];li[la]=li[la]/lq[la][la];}ll;}
|
||||
|
|
@ -0,0 +1,40 @@
|
|||
/**********************************************************************/
|
||||
/*** ***/
|
||||
/*** vorrueckwaertsubs.h ***/
|
||||
/*** ***/
|
||||
/*** Vorwaerts und Rueckwaertssubstitution ***/
|
||||
/*** ***/
|
||||
/**********************************************************************/
|
||||
/*** ***/
|
||||
/*** Autorin: Luise Blank ***/
|
||||
/*** ***/
|
||||
/*** Datum: 2.11.2011 ***/
|
||||
/*** ***/
|
||||
/**********************************************************************/
|
||||
|
||||
/*
|
||||
vwsubs
|
||||
Diese Funktion loest das Gleichungssystem Lx=b,
|
||||
wobei L eine normierte untere Dreiecksmatrix ist mittels
|
||||
Vorwaertssubstitution. Dabei wird die rechte Seite b mit
|
||||
der Loesung ueberschrieben.
|
||||
input: int n Dimension der Matrix
|
||||
double** L Zeiger auf normierte untere Dreiecksmatrix
|
||||
double* b Rechte Seite und Ausgabe
|
||||
output: void
|
||||
*/
|
||||
void vwsubs(int n, double ** L, double * b);
|
||||
|
||||
|
||||
/*
|
||||
rwsubs
|
||||
Diese Funktion loest das Gleichungssystem Rx=b,
|
||||
wobei R eine regulaere obere Dreiecksmatrix ist mittels
|
||||
Rueckwaertssubstitution. Dabei wird die rechte Seite b
|
||||
mit der Loesung ueberschrieben.
|
||||
input: int n Dimension der Matrix
|
||||
double** R Zeiger auf obere Dreiecksmatrix
|
||||
double* b Rechte Seite und Ausgabe
|
||||
output: void
|
||||
*/
|
||||
void rwsubs(int n, double ** R, double * b);
|
||||
Loading…
Reference in New Issue