Changeset 682 in Sophya for trunk/SophyaLib/NTools
- Timestamp:
- Dec 10, 1999, 5:56:03 PM (26 years ago)
- Location:
- trunk/SophyaLib/NTools
- Files:
-
- 7 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/SophyaLib/NTools/datime.c
r517 r682 6 6 #include <math.h> 7 7 #ifdef __MWERKS__ 8 #include "mwerksmath.h" 9 // #include "unixmac.h" 8 #include "unixmac.h" 10 9 #endif 11 10 -
trunk/SophyaLib/NTools/generalfit.cc
r519 r682 5 5 #include <math.h> 6 6 #ifdef __MWERKS__ 7 #include "mwerksmath.h" 8 // #include "unixmac.h" 7 #include "unixmac.h" 9 8 #endif 10 9 #include <string.h> -
trunk/SophyaLib/NTools/histos.cc
r514 r682 1 1 // 2 // $Id: histos.cc,v 1. 6 1999-10-25 16:39:57ansari Exp $2 // $Id: histos.cc,v 1.7 1999-12-10 16:55:51 ansari Exp $ 3 3 // 4 4 … … 1781 1781 1782 1782 // Que faut-il ecrire? 1783 int errok = (dobj->err2) ? 1 : 0;1783 int_4 errok = (dobj->err2) ? 1 : 0; 1784 1784 1785 1785 // Ecriture entete pour identifier facilement -
trunk/SophyaLib/NTools/matxop.c
r220 r682 393 393 if (szv < 25) szv = 25; 394 394 395 if ( (FVeci = malloc(nvarmx*sizeof(float *))) == NULL )395 if ( (FVeci = (float**) malloc(nvarmx*sizeof(float *))) == NULL ) 396 396 { printf("InitRFitLin_Erreur: (Pb malloc(Veci)) \n"); 397 397 return(1); 398 398 } 399 399 400 if ( (FVSpace = malloc(nvarmx*szv*sizeof(float))) == NULL )400 if ( (FVSpace = (float*) malloc(nvarmx*szv*sizeof(float))) == NULL ) 401 401 { printf("InitRFitLin_Erreur: (Pb malloc(VSpace)) \n"); 402 402 free(FVeci); … … 407 407 408 408 409 if ( (FVecf = malloc(3*sizeof(float *))) == NULL )409 if ( (FVecf = (float**) malloc(3*sizeof(float *))) == NULL ) 410 410 { printf("InitRFitLin_Erreur: (Pb malloc(Vecf)) \n"); 411 411 free(FVeci); … … 414 414 } 415 415 416 if ( (FVSpacef = malloc(3*szv*sizeof(float))) == NULL )416 if ( (FVSpacef = (float*) malloc(3*szv*sizeof(float))) == NULL ) 417 417 { printf("InitRFitLin_Erreur: (Pb malloc(VSpacef)) \n"); 418 418 free(FVeci); … … 425 425 426 426 427 if ( (FFMtx = malloc(nvarmx*nvarmx*sizeof(double))) == NULL )427 if ( (FFMtx = (double*) malloc(nvarmx*nvarmx*sizeof(double))) == NULL ) 428 428 { printf("InitRFitLin_Erreur: (Pb malloc(FMtx)) \n"); 429 429 free(FVeci); … … 434 434 } 435 435 436 if ( (FFMtx2 = malloc(nvarmx*nvarmx*sizeof(double))) == NULL )436 if ( (FFMtx2 = (double*) malloc(nvarmx*nvarmx*sizeof(double))) == NULL ) 437 437 { printf("InitRFitLin_Erreur: (Pb malloc(FMtx2)) \n"); 438 438 free(FVeci); … … 444 444 } 445 445 446 if ( (FFSort = malloc((nvarmx+1)*nvarmx*sizeof(double))) == NULL )446 if ( (FFSort = (double*) malloc((nvarmx+1)*nvarmx*sizeof(double))) == NULL ) 447 447 { printf("InitRFitLin_Erreur: (Pb malloc(FFSort)) \n"); 448 448 free(FVeci); … … 705 705 if (szv < 25) szv = 25; 706 706 707 if ( (DVeci = malloc(nvarmx*sizeof(double *))) == NULL )707 if ( (DVeci = (double**) malloc(nvarmx*sizeof(double *))) == NULL ) 708 708 { printf("InitDFitLin_Erreur: (Pb malloc(Veci)) \n"); 709 709 return(1); 710 710 } 711 711 712 if ( (DVSpace = malloc(nvarmx*szv*sizeof(double))) == NULL )712 if ( (DVSpace = (double*) malloc(nvarmx*szv*sizeof(double))) == NULL ) 713 713 { printf("InitDFitLin_Erreur: (Pb malloc(VSpace)) \n"); 714 714 free(DVeci); … … 719 719 720 720 721 if ( (DVecf = malloc(3*sizeof(double *))) == NULL )721 if ( (DVecf = (double**) malloc(3*sizeof(double *))) == NULL ) 722 722 { printf("InitDFitLin_Erreur: (Pb malloc(Vecf)) \n"); 723 723 free(DVeci); … … 726 726 } 727 727 728 if ( (DVSpacef = malloc(3*szv*sizeof(double))) == NULL )728 if ( (DVSpacef = (double*) malloc(3*szv*sizeof(double))) == NULL ) 729 729 { printf("InitDFitLin_Erreur: (Pb malloc(VSpacef)) \n"); 730 730 free(DVeci); … … 737 737 738 738 739 if ( (DFMtx = malloc(nvarmx*nvarmx*sizeof(double))) == NULL )739 if ( (DFMtx = (double*) malloc(nvarmx*nvarmx*sizeof(double))) == NULL ) 740 740 { printf("InitDFitLin_Erreur: (Pb malloc(FMtx)) \n"); 741 741 free(DVeci); … … 746 746 } 747 747 748 if ( (DFMtx2 = malloc(nvarmx*nvarmx*sizeof(double))) == NULL )748 if ( (DFMtx2 = (double*) malloc(nvarmx*nvarmx*sizeof(double))) == NULL ) 749 749 { printf("InitDFitLin_Erreur: (Pb malloc(FMtx2)) \n"); 750 750 free(DVeci); -
trunk/SophyaLib/NTools/nbmath.c
r220 r682 1643 1643 { 1644 1644 int npt,k,i,nclass; 1645 double *clas s,aa1,c2,cut;1645 double *clas,aa1,c2,cut; 1646 1646 1647 1647 *a1 =0.; … … 1654 1654 /* printf("pass 1: %g*x c2=%g/%d\n",*a1,c2,npt); */ 1655 1655 1656 clas s =malloc( *n * sizeof(double) );1657 if( clas s== NULL ) {*n=npt; return(-3.);}1656 clas = (double*) malloc( *n * sizeof(double) ); 1657 if( clas == NULL ) {*n=npt; return(-3.);} 1658 1658 1659 1659 /* elimination des mauvais points */ … … 1662 1662 if(ey[i]<=0.) continue; 1663 1663 c2 = (y[i]-aa1*x[i])/ey[i]; 1664 clas s[nclass] = c2*c2;1664 clas[nclass] = c2*c2; 1665 1665 nclass++; 1666 1666 } 1667 qsort(clas s,(size_t) nclass,(size_t) sizeof(double),qSort_Dble);1667 qsort(clas,(size_t) nclass,(size_t) sizeof(double),qSort_Dble); 1668 1668 k = (int) ( (1. - per_clean ) * (double) nclass ); 1669 1669 if(k<0) k=0; 1670 1670 if(k>=nclass) k = nclass-1; 1671 cut = clas s[k];1671 cut = clas[k]; 1672 1672 k = 0; 1673 1673 for(i=0;i<*n;i++) { 1674 clas s[i] = ey[i];1674 clas[i] = ey[i]; 1675 1675 c2 = (y[i]-aa1*x[i])/ey[i]; 1676 1676 c2 *= c2; 1677 if(ey[i]>0. && c2>cut) {clas s[i] = -1.; k++;}1677 if(ey[i]>0. && c2>cut) {clas[i] = -1.; k++;} 1678 1678 } 1679 1679 /* printf("nombre pt tues %d cut=%g\n",k,cut); */ … … 1681 1681 /* 2sd passe */ 1682 1682 npt = *n; 1683 c2 = FitProp(x,y,clas s,&npt,&aa1);1684 if(c2<0.) {*n = npt; free(clas s); return(-2.);}1683 c2 = FitProp(x,y,clas,&npt,&aa1); 1684 if(c2<0.) {*n = npt; free(clas); return(-2.);} 1685 1685 *a1 = aa1; 1686 1686 *n = npt; 1687 1687 /* printf("pass 2: %g*x c2=%g/%d\n",*a1,c2,npt); */ 1688 1688 1689 free(clas s);1689 free(clas); 1690 1690 return(c2); 1691 1691 } … … 1730 1730 { 1731 1731 int npt,k,i,nclass; 1732 double *clas s,aa0,aa1,c2,cut;1732 double *clas,aa0,aa1,c2,cut; 1733 1733 1734 1734 *a0 = *a1 =0.; … … 1742 1742 /* printf("pass 1: %g + %g*x c2=%g/%d\n",*a0,*a1,c2,npt); */ 1743 1743 1744 clas s =malloc( *n * sizeof(double) );1745 if( clas s== NULL ) {*n=npt; return(-3.);}1744 clas = (double*) malloc( *n * sizeof(double) ); 1745 if( clas == NULL ) {*n=npt; return(-3.);} 1746 1746 1747 1747 /* elimination des mauvais points */ … … 1750 1750 if(ey[i]<=0.) continue; 1751 1751 c2 = (y[i]-(aa0+aa1*x[i]))/ey[i]; 1752 clas s[nclass] = c2*c2;1752 clas[nclass] = c2*c2; 1753 1753 nclass++; 1754 1754 } 1755 qsort(clas s,(size_t) nclass,(size_t) sizeof(double),qSort_Dble);1755 qsort(clas,(size_t) nclass,(size_t) sizeof(double),qSort_Dble); 1756 1756 k = (int) ( (1. - per_clean ) * (double) nclass ); 1757 1757 if(k<0) k=0; 1758 1758 if(k>=nclass) k = nclass-1; 1759 cut = clas s[k];1759 cut = clas[k]; 1760 1760 k = 0; 1761 1761 for(i=0;i<*n;i++) { 1762 clas s[i] = ey[i];1762 clas[i] = ey[i]; 1763 1763 c2 = (y[i]-(aa0+aa1*x[i]))/ey[i]; 1764 1764 c2 *= c2; 1765 if(ey[i]>0. && c2>cut) {clas s[i] = -1.; k++;}1765 if(ey[i]>0. && c2>cut) {clas[i] = -1.; k++;} 1766 1766 } 1767 1767 /* printf("nombre pt tues %d cut=%g\n",k,cut); */ … … 1769 1769 /* 2sd passe */ 1770 1770 npt = *n; 1771 c2 = FitLin(x,y,clas s,&npt,&aa0,&aa1);1772 if(c2<0.) {*n = npt; free(clas s); return(-2.);}1771 c2 = FitLin(x,y,clas,&npt,&aa0,&aa1); 1772 if(c2<0.) {*n = npt; free(clas); return(-2.);} 1773 1773 *a0 = aa0; 1774 1774 *a1 = aa1; … … 1776 1776 /* printf("pass 2: %g + %g*x c2=%g/%d\n",*a0,*a1,c2,npt); */ 1777 1777 1778 free(clas s);1778 free(clas); 1779 1779 return(c2); 1780 1780 } … … 1815 1815 { 1816 1816 int npt,k,i,nclass; 1817 double *clas s,aa0,aa1,aa2,c2,cut;1817 double *clas,aa0,aa1,aa2,c2,cut; 1818 1818 1819 1819 *a0 = *a1 = *a2 =0.; … … 1828 1828 /* printf("pass 1: %g + %g*x + %g*x**2 c2=%g/%d\n",*a0,*a1,*a2,c2,npt); */ 1829 1829 1830 clas s =malloc( *n * sizeof(double) );1831 if( clas s== NULL ) {*n=npt; return(-3.);}1830 clas = (double*) malloc( *n * sizeof(double) ); 1831 if( clas == NULL ) {*n=npt; return(-3.);} 1832 1832 1833 1833 /* elimination des mauvais points */ … … 1836 1836 if(ey[i]<=0.) continue; 1837 1837 c2 = (y[i]-(aa0+aa1*x[i]+aa2*x[i]*x[i]))/ey[i]; 1838 clas s[nclass] = c2*c2;1838 clas[nclass] = c2*c2; 1839 1839 nclass++; 1840 1840 } 1841 qsort(clas s,(size_t) nclass,(size_t) sizeof(double),qSort_Dble);1841 qsort(clas,(size_t) nclass,(size_t) sizeof(double),qSort_Dble); 1842 1842 k = (int) ( (1. - per_clean ) * (double) nclass ); 1843 1843 if(k<0) k=0; 1844 1844 if(k>=nclass) k = nclass-1; 1845 cut = clas s[k];1845 cut = clas[k]; 1846 1846 k = 0; 1847 1847 for(i=0;i<*n;i++) { 1848 clas s[i] = ey[i];1848 clas[i] = ey[i]; 1849 1849 c2 = (y[i]-(aa0+aa1*x[i]+aa2*x[i]*x[i]))/ey[i]; 1850 1850 c2 *= c2; 1851 if(ey[i]>0. && c2>cut) {clas s[i] = -1.; k++;}1851 if(ey[i]>0. && c2>cut) {clas[i] = -1.; k++;} 1852 1852 } 1853 1853 /* printf("nombre pt tues %d cut=%g\n",k,cut); */ … … 1855 1855 /* 2sd passe */ 1856 1856 npt = *n; 1857 c2 = FitPar(x,y,clas s,&npt,&aa0,&aa1,&aa2);1858 if(c2<0.) {*n = npt; free(clas s); return(-2.);}1857 c2 = FitPar(x,y,clas,&npt,&aa0,&aa1,&aa2); 1858 if(c2<0.) {*n = npt; free(clas); return(-2.);} 1859 1859 *a0 = aa0; 1860 1860 *a1 = aa1; … … 1863 1863 /* printf("pass 2: %g + %g*x + %g*x**2 c2=%g/%d\n",*a0,*a1,*a2,c2,npt); */ 1864 1864 1865 free(clas s);1865 free(clas); 1866 1866 return(c2); 1867 1867 } -
trunk/SophyaLib/NTools/nbrandom.c
r244 r682 495 495 496 496 sof = sizeof(struct tirage_alea); 497 if( (t = malloc(sof) ) == NULL ) {497 if( (t = (struct tirage_alea*)malloc(sof) ) == NULL ) { 498 498 printf("impossible d'allouer *tirage_alea par malloc \n"); 499 499 return(NULL); … … 503 503 504 504 sof = nbin * sizeof(double); 505 if( (t->Tab = malloc(sof) ) == NULL ) {505 if( (t->Tab = (double*)malloc(sof) ) == NULL ) { 506 506 printf("impossible d'allouer *tirage_alea.Tab par malloc \n"); 507 507 return(NULL); -
trunk/SophyaLib/NTools/nbrandom.h
r517 r682 8 8 #include <stdlib.h> 9 9 #ifdef __MWERKS__ 10 #include "mwerksmath.h" 11 #include "unixmac.h" 10 #include "unixmac.h" 12 11 #endif 13 12
Note:
See TracChangeset
for help on using the changeset viewer.