begin process at 2010 03 18 23:42:55
  Trouver un code source :
 
dans
 
Accueil > 

Code

 > 

Maths et Algorithmes

 > CALCUL DE PI

CALCUL DE PI


 Information sur la source

Note :
Aucune note
Catégorie :Maths et Algorithmes Classé sous :calculer, pi, mathématique, calcul, nombre Niveau :Débutant Date de création :06/02/2005 Date de mise à jour :06/02/2005 12:00:33 Vu :6 287

Auteur : cyberfastfred

Ecrire un message privé
Commentaire sur cette source (1)
Ajouter un commentaire et/ou une note

 Description

Ce petit programme calule Pi avec un grand nombre de décimales.
J'utilise une série simple qui converge assez rapidement (regarder le code de la classe calculpi)
L'implémentation du calcul sur beaucoup de décimales est minimaliste. J'ai juste fait ce qui servait à la série.
L'objectif de ce prog était de mesurer les perfs de Java sur des algo très itératif et de comparer avec le C. (Le prog C est à venir).

Source

  • /* fichier calcpi.java
  • * Created on 22 janv. 2005
  • * @author cyberfastfred
  • *
  • * Calcul de pi avec une serie simple
  • * objectif comparer les perf entre Java et C sur des algo
  • * tres itératifs
  • */
  • import java.util.Date;
  • public class calcpi {
  • public static void main(String[] args) {
  • int i;
  • UDChar pi, frac;
  • java.lang.String StrPI="PI=";
  • Date now = new Date();
  • System.out.println(now);
  • pi = new UDChar(2);
  • frac = new UDChar(2);
  • for (i=1; i <3000; i++){
  • frac.mulEntiere(i);
  • frac.divEntiere(2*i +1);
  • pi.add(frac);
  • }
  • now = new Date();
  • System.out.println(now);
  • System.out.println("PI="+ pi);
  • System.out.println("frac="+ frac);
  • }
  • }
  • /**********************fichier UDChar.java**************************/
  • /**
  • * @author cyberfastfred
  • * Created on 23 janv. 2005
  • * implémentation basique orientée pour le calcul de pi à 100 décimale
  • * La série calculant pi n'ajoutant que des nombre fractionnaires de plus en plus petit
  • * Les grands nombres sont uniquement positifs
  • */
  • public class UDChar {
  • /*nombre de chiffres*/
  • public static final int PRECISION = 1000;
  • protected char [] chiffres = new char [PRECISION];
  • /* 10 chiffres avant la virgule PRECISION-10 chiffres après */
  • /*créé un grand décimal à partir d'un int
  • *
  • */
  • public UDChar(int entier){
  • int i;
  • /*
  • * Initiailisation de la partie entière
  • */
  • for (i=9; i>=0; i--){
  • int chiffre;
  • chiffre = entier % 10;
  • chiffres[i] = intToChiffre(chiffre);
  • entier /=10;
  • }
  • /*mise à 0 des décimales*/
  • for (i=10; i<PRECISION; i++) chiffres[i]='0';
  • }
  • /*
  • * Ajoute l'operande Ultradecmal
  • */
  • public void add(UDChar operande) {
  • int i, somme, retenue;
  • retenue = 0;
  • for(i = PRECISION - 1; i >= 0; i--){
  • somme = chiffreToInt(chiffres[i])
  • + chiffreToInt(operande.chiffres[i])
  • + retenue;
  • retenue = somme / 10;
  • somme %= 10;
  • chiffres[i] = intToChiffre(somme);
  • }
  • }
  • /*
  • * Soustrait l'opérande UltraDecimal
  • */
  • public void sub(UDChar operande) {
  • int i, retrait, retenue;
  • retenue = 0;
  • for(i = PRECISION - 1; i >= 0; i--){
  • retrait = 10 + chiffreToInt(chiffres[i])
  • - chiffreToInt(operande.chiffres[i])
  • - retenue;
  • retenue = 1 - (retrait / 10);
  • retrait %= 10;
  • chiffres[i] = intToChiffre(retrait);
  • }
  • }
  • /*
  • * multiplication d'un UD par un entier
  • * suffisant pour le calcul de séries simples
  • */
  • public void mulEntiere(int operande){
  • UDChar partiel, original;
  • int i, chiffre, operandeSuiv;
  • original = this.dupl();
  • /* remise à zero de l'UD */
  • for (i=0; i<PRECISION; i++) chiffres[i]='0';
  • operandeSuiv = operande;
  • for (i=0; i <10 ;i++){
  • partiel = original.dupl();
  • chiffre = operandeSuiv % 10;
  • operandeSuiv /=10;
  • partiel.mulUnit(chiffre);
  • this.add(partiel);
  • original.shiftLeft(1);
  • }
  • }
  • /*
  • * multiplication d'un UD par un entier
  • * suffisant pour le calcul de séries simples
  • */
  • public void divEntiere(int operande){
  • int i, reste, result;
  • reste = 0;
  • result = 0;
  • for (i = 0; i <PRECISION; i++){
  • reste = reste * 10 + chiffreToInt(chiffres[i]);
  • result = reste / operande;
  • reste = reste % operande;
  • chiffres[i] = intToChiffre(result);
  • }
  • }
  • /*
  • * Renvoie une chaine cotenant le nombre représenté
  • */
  • public String toString() {
  • String str= new String("");
  • int i;
  • for(i=0; i<10; i++){
  • str = str + chiffres[i];
  • }
  • str = str +",";
  • for(i=10; i<PRECISION; i++){
  • if (i%80==0)str=str+"\n";
  • str = str + chiffres[i];
  • }
  • return str;
  • }
  • /*
  • * clone un UD
  • */
  • private UDChar dupl(){
  • UDChar ud = new UDChar(0);
  • ud.add(this); /* pas optimal mais facile */
  • return ud;
  • }
  • /*
  • * decalage à droite
  • * equivalent à multiplication par 10 puissance pas
  • */
  • private void shiftRight(int pas){
  • int i;
  • if (pas <= 0) return;
  • if (pas > PRECISION) pas = PRECISION;
  • for (i = pas; i <PRECISION; i++){
  • chiffres[i] = chiffres[i-pas];
  • }
  • for (i = 0; i < pas; i++){
  • chiffres[i] = '0';
  • }
  • }
  • /*
  • * decalage à gauche
  • * equivalent à multiplication par 10 puissance pas
  • */
  • private void shiftLeft(int pas){
  • int i;
  • if (pas <= 0) return;
  • if (pas > PRECISION) pas = PRECISION;
  • for (i = 0; i < PRECISION-pas; i++){
  • chiffres[i] = chiffres[i+pas];
  • }
  • for (i = PRECISION-pas; i <PRECISION; i++){
  • chiffres[i] = '0';
  • }
  • }
  • /*
  • * convertie un chiffre sous forme de char en entier
  • */
  • private int chiffreToInt(char chiffre){
  • switch(chiffre){
  • case '0':
  • return(0);
  • case '1':
  • return(1);
  • case '2':
  • return(2);
  • case '3':
  • return(3);
  • case '4':
  • return(4);
  • case '5':
  • return(5);
  • case '6':
  • return(6);
  • case '7':
  • return(7);
  • case '8':
  • return(8);
  • case '9':
  • return(9);
  • default:
  • return(0);
  • }
  • }
  • /*
  • * Convertie un entier de 0 à 9 en char
  • */
  • private char intToChiffre(int chiffre){
  • switch(chiffre){
  • case 0:
  • return('0');
  • case 1:
  • return('1');
  • case 2:
  • return('2');
  • case 3:
  • return('3');
  • case 4:
  • return('4');
  • case 5:
  • return('5');
  • case 6:
  • return('6');
  • case 7:
  • return('7');
  • case 8:
  • return('8');
  • case 9:
  • return('9');
  • default:
  • return('0');
  • }
  • }
  • /*
  • * multiplie par un chiffre de 0 à 9
  • */
  • private void mulUnit(int operande){
  • int chiffre, i, resultat, retenue;
  • chiffre = operande%10;
  • retenue = 0;
  • for(i = PRECISION - 1; i >= 0; i--){
  • resultat = chiffreToInt(chiffres[i])
  • * chiffre
  • + retenue;
  • retenue = resultat / 10;
  • resultat %= 10;
  • chiffres[i] = intToChiffre(resultat);
  • }
  • }
  • }
/* fichier calcpi.java
 * Created on 22 janv. 2005
 * @author cyberfastfred
 * 
 * Calcul de pi avec une serie simple 
 * objectif comparer les perf entre Java et C sur des algo
 * tres itératifs
 */
import java.util.Date;

public class calcpi {

	public static void main(String[] args) {
		int i;
		UDChar pi, frac;
		java.lang.String StrPI="PI="; 
		Date now = new Date();
		System.out.println(now);
		 
		pi = new UDChar(2);
		frac = new UDChar(2);
		for (i=1; i <3000; i++){
			frac.mulEntiere(i);
			frac.divEntiere(2*i +1);
			pi.add(frac);
		}
        now = new Date();
		System.out.println(now);
		System.out.println("PI="+ pi);
		System.out.println("frac="+ frac);	
	}
}
/**********************fichier UDChar.java**************************/
/**
 * @author cyberfastfred
 * Created on 23 janv. 2005
 * implémentation basique orientée pour le calcul de pi à 100 décimale
 * La série calculant pi n'ajoutant que des nombre fractionnaires de plus en plus petit
 * Les grands nombres sont uniquement positifs
 */
public class UDChar  {
	/*nombre de chiffres*/
	public static final int PRECISION = 1000;

	protected char [] chiffres = new char [PRECISION];
	/* 10 chiffres avant la virgule PRECISION-10 chiffres après */
	
	/*créé un grand décimal à partir d'un int
	 * 
	 */   	
    public UDChar(int entier){
    	int i;
        /*
         * Initiailisation de la partie entière
         */
         for (i=9; i>=0; i--){
         	int chiffre;
         	chiffre = entier % 10;
         	chiffres[i] = intToChiffre(chiffre);

         	entier /=10;
         }
         /*mise à 0 des décimales*/
         for (i=10; i<PRECISION; i++) chiffres[i]='0';
    }
    
	/*
	 * Ajoute l'operande Ultradecmal
	 */
	public void add(UDChar operande) {
		int i, somme, retenue;
		retenue = 0;
		for(i = PRECISION - 1; i >= 0; i--){
			somme = chiffreToInt(chiffres[i])
			      + chiffreToInt(operande.chiffres[i])
				  + retenue;
			retenue = somme / 10;
			somme %= 10;
			chiffres[i] = intToChiffre(somme);
		}
	}

	/*
	 * Soustrait l'opérande UltraDecimal
	 */
	public void sub(UDChar operande) {
		int i, retrait, retenue;
		retenue = 0;
		for(i = PRECISION - 1; i >= 0; i--){
			retrait = 10 + chiffreToInt(chiffres[i])
			      - chiffreToInt(operande.chiffres[i])
				  - retenue;
			retenue = 1 - (retrait / 10);
			retrait %= 10;
			chiffres[i] = intToChiffre(retrait);
		}
		
	}

	/*
	 * multiplication d'un UD par un entier
	 * suffisant pour le calcul de séries simples
	 */
	public void mulEntiere(int operande){
		UDChar partiel, original;
		int i, chiffre, operandeSuiv;
		original = this.dupl();
		
		/* remise à zero de l'UD */
		for (i=0; i<PRECISION; i++) chiffres[i]='0';
		
		operandeSuiv = operande;
		for (i=0; i <10 ;i++){
			partiel = original.dupl();
			chiffre = operandeSuiv % 10;
			operandeSuiv /=10;
			partiel.mulUnit(chiffre);
			this.add(partiel);
			original.shiftLeft(1);
		}
	}

	/*
	 * multiplication d'un UD par un entier
	 * suffisant pour le calcul de séries simples
	 */
	public void divEntiere(int operande){
		int i, reste, result;
		reste = 0;
		result = 0;
		for (i = 0; i <PRECISION; i++){
			reste = reste * 10 + chiffreToInt(chiffres[i]);
			result = reste / operande;
			reste = reste % operande;
			chiffres[i] = intToChiffre(result);
		}
    	
	}

	/*
	 * Renvoie une chaine cotenant le nombre représenté 
	 */
	public String toString() {
		String str= new String("");
		int i;
		for(i=0; i<10; i++){
			str = str + chiffres[i];
		}
		str = str +",";
		for(i=10; i<PRECISION; i++){
			if (i%80==0)str=str+"\n";
			str = str + chiffres[i];
		}
		return str;
	}
	/*
	 * clone un UD
	 */
	private UDChar dupl(){
		UDChar ud = new UDChar(0);
		ud.add(this); /* pas optimal mais facile */
		return ud;
	}
	
	/*
	 * decalage à droite
	 * equivalent à multiplication par 10 puissance pas
	 */
	private void shiftRight(int pas){
		int i;
		if (pas <= 0) return;
		if (pas > PRECISION) pas = PRECISION;
		for (i = pas; i <PRECISION; i++){
			chiffres[i] = chiffres[i-pas];
		}
		for (i = 0; i < pas; i++){
			chiffres[i] = '0';
		}
	}
	
	/*
	 * decalage à gauche
	 * equivalent à multiplication par 10 puissance pas
	 */
	private void shiftLeft(int pas){
		int i;
		if (pas <= 0) return;
		if (pas > PRECISION) pas = PRECISION;
		for (i = 0; i < PRECISION-pas; i++){
			chiffres[i] = chiffres[i+pas];
		}
		for (i = PRECISION-pas; i <PRECISION; i++){
			chiffres[i] = '0';
		}
	}
	
	/*
	 * convertie un chiffre sous forme de char en entier
	 */
	private int chiffreToInt(char chiffre){
		switch(chiffre){
			case '0':
				return(0);
			case '1':
				return(1);
			case '2':
				return(2);
			case '3':
				return(3);
			case '4':
				return(4);
			case '5':
				return(5);
			case '6':
				return(6);
			case '7':
				return(7);
			case '8':
				return(8);
			case '9':
				return(9);
			default:
				return(0);
		}
	}
	
	/*
	 * Convertie un entier de 0 à 9 en char
	 */
	private char intToChiffre(int chiffre){
		switch(chiffre){
			case 0:
			  return('0');
			case 1:
				  return('1');
			case 2:
				  return('2');
			case 3:
				  return('3');
			case 4:
				  return('4');
			case 5:
				  return('5');
			case 6:
				  return('6');
			case 7:
				  return('7');
			case 8:
				  return('8');
			case 9:
				  return('9');
			default:
				return('0');
		}
	}
	
	/*
	 * multiplie par un chiffre de 0 à 9
	 */
	private void mulUnit(int operande){
		int chiffre, i, resultat, retenue;
		chiffre = operande%10;
		retenue = 0;
		for(i = PRECISION - 1; i >= 0; i--){
			resultat = chiffreToInt(chiffres[i])
				  * chiffre
				  + retenue;
			retenue = resultat / 10;
			resultat %= 10;
			chiffres[i] = intToChiffre(resultat);
		}
	}

}

 Conclusion

Pour le calcul sur plusieurs dizaine de décimales j'utilise un tableau de caractères.
Arbitrairement les 10 premiers sont avant la virgule et les suivants sont après.
Les UDChar représentent donc des nombres positifs < 10000000000.
Il n'y a aucun controle de dépassement.
Le code n'est absolument pas optimisé puisque je parcours l'intégralité les chiffres sans  chercher à éliminer les 0.
Je n'utilise aucun outils spécifique java puisque je veux avoir une implémentation proche en C pour comparer les vitesses.  (En plus ça m'arrange car je ne connais presque rien à Java).


 Historique

06 février 2005 12:00:34 :

 Sources de la même categorie

TEST DE PRIMALITÉ OPTIMISÉ par Julien39
Source avec Zip CLASSE MATRICE par frankladen11
Source avec Zip Source avec une capture RÉSOLUTION D'ÉQUATION GRÂCE AU CALCUL DES DÉTERMINANTS par frankladen11
Source avec Zip TYPE DE DONNÉES ABSTRAIT GRAPHE par smutsonberg
Source avec Zip Source avec une capture SIMPLEXE ET DUAL par MrRenaud

 Sources en rapport avec celle ci

Source avec Zip Source avec une capture CALCULER UNE ÉQUATION MATHÉMATIQUE COMPLEXE DANS UN STRING E... par DeadlyPredator
Source avec Zip TOUTES LES OPÉRATIONS POSSIBLES SUR UNE MATRICE par Captain Igloo
Source avec Zip CALCULATRICE par javajoum
Source avec Zip ADDITION DE DEUX NOMBRES RÉELS "COMME A LA MATERNELLE" par zinotron
CALCUL DE NOMBRES PREMIERS par niiizzzz

Commentaires et avis

Commentaire de Syruis le 20/09/2007 16:48:04

Bon ce code est magique... 10/10 mais une question qui me trote dans la tête pourquoi ...

On ne peut pas utiliser la fonctione Integer.valueOf pour aller d'un char vers une String ??

Merci,
Cordialement,
Syruis

 Ajouter un commentaire


Discussions en rapport avec ce code source dans le forum

calculer le nombre des objets identiques dans un vecteur [ par Meriem85 ] Bonsoir,J'ai un vecteur de Personne(name,pre,age),dans mon code je veux calculer le nombre des cas dans lesquelles il y a des personnes ds la vecteur calcul de nombre d'accolade [ par zeinebfki ] bonjour,j'ai un probléme qui consiste à délimiter les différents éléments d'un code en java.j'ai transformée le code en une chaine de caractére et je Automate java (calcul d'une expression mathématique simple) [ par Raph113 ] Bonjour à tous, je m'adresse à vous car j'ai un projet en java à réaliser et je suis un desesperé. Voila donc notre projet consiste à réaliser une app calcul d'integrale triples [ par soaaziz ] Bonjour,je suis un ingénieur en informatique et je suis intégré dans un projet purement mathématiquema mission consiste a réaliser un programme qui pe calcul d'integrale triples [ par soaaziz ] Bonjour, je suis un ingénieur en informatique et je suis intégré dans un projet purement mathématique ma mission consiste a réaliser un programme qui calcul d'une surface [ par mersai ] bonjour, dans un repère cartésien, j'ai un ensemble de points dont les coordonnées x et y sont connues. ils sont a l'intérieur d'un rectangle dont la Détecter le nombre de groups 1 dans une matrice donnée [ par thuyte ] Salut salut à toutDepuis quelques jours je cherche un algo en java quii retourn une liste de groupe de "1"  dans une matrice donnée. Par ex: Matrice d trouver le nombre d'arcs et de sommets dans un réseau [ par papillon2000 ] bonsoir à tous, je veux calculer et trouvez le nombre d'arcs et de sommets d'un réseau ou un graphe de contraintes. merci pour votre aide.


Nos sponsors


Sondage...

CalendriCode

Mars 2010
LMMJVSD
1234567
891011121314
15161718192021
22232425262728
293031    

Consulter la suite du CalendriCode

Photothèque

 
Développement réalisé par Nicolas SOREL (Nix) avec l'aide de : Cyril DURAND et Emmanuel (EBArtSoft), Merci à Vincent pour ses précieux conseils.
CodeS-SourceS.com© Toute reproduction même partielle est interdite sauf accord écrit du Webmaster
CodeS-SourceS.com© est une marque déposée tous droits réservés

Google Coop CodeS-SourceS Google Coop CodeS-SourceS
Temps d'éxécution de la page : 0,874 sec (3)

Nous contacter | Annoncer sur CodeS-SourceS | Mentions légales