cximage/src/jasper/jpc/jpc_util.c - toolchain/benchmark - Git at Google

 /*
  * Copyright (c) 1999-2000 Image Power, Inc. and the University of
  *   British Columbia.
  * Copyright (c) 2001-2003 Michael David Adams.
  * All rights reserved.
  */

 /* __START_OF_JASPER_LICENSE__
  *
  * JasPer Software License
  *
  * IMAGE POWER JPEG-2000 PUBLIC LICENSE
  * ************************************
  *
  * GRANT:
  *
  * Permission is hereby granted, free of charge, to any person (the "User")
  * obtaining a copy of this software and associated documentation, to deal
  * in the JasPer Software without restriction, including without limitation
  * the right to use, copy, modify, merge, publish, distribute, sublicense,
  * and/or sell copies of the JasPer Software (in source and binary forms),
  * and to permit persons to whom the JasPer Software is furnished to do so,
  * provided further that the License Conditions below are met.
  *
  * License Conditions
  * ******************
  *
  * A.  Redistributions of source code must retain the above copyright notice,
  * and this list of conditions, and the following disclaimer.
  *
  * B.  Redistributions in binary form must reproduce the above copyright
  * notice, and this list of conditions, and the following disclaimer in
  * the documentation and/or other materials provided with the distribution.
  *
  * C.  Neither the name of Image Power, Inc. nor any other contributor
  * (including, but not limited to, the University of British Columbia and
  * Michael David Adams) may be used to endorse or promote products derived
  * from this software without specific prior written permission.
  *
  * D.  User agrees that it shall not commence any action against Image Power,
  * Inc., the University of British Columbia, Michael David Adams, or any
  * other contributors (collectively "Licensors") for infringement of any
  * intellectual property rights ("IPR") held by the User in respect of any
  * technology that User owns or has a right to license or sublicense and
  * which is an element required in order to claim compliance with ISO/IEC
  * 15444-1 (i.e., JPEG-2000 Part 1).  "IPR" means all intellectual property
  * rights worldwide arising under statutory or common law, and whether
  * or not perfected, including, without limitation, all (i) patents and
  * patent applications owned or licensable by User; (ii) rights associated
  * with works of authorship including copyrights, copyright applications,
  * copyright registrations, mask work rights, mask work applications,
  * mask work registrations; (iii) rights relating to the protection of
  * trade secrets and confidential information; (iv) any right analogous
  * to those set forth in subsections (i), (ii), or (iii) and any other
  * proprietary rights relating to intangible property (other than trademark,
  * trade dress, or service mark rights); and (v) divisions, continuations,
  * renewals, reissues and extensions of the foregoing (as and to the extent
  * applicable) now existing, hereafter filed, issued or acquired.
  *
  * E.  If User commences an infringement action against any Licensor(s) then
  * such Licensor(s) shall have the right to terminate User's license and
  * all sublicenses that have been granted hereunder by User to other parties.
  *
  * F.  This software is for use only in hardware or software products that
  * are compliant with ISO/IEC 15444-1 (i.e., JPEG-2000 Part 1).  No license
  * or right to this Software is granted for products that do not comply
  * with ISO/IEC 15444-1.  The JPEG-2000 Part 1 standard can be purchased
  * from the ISO.
  *
  * THIS DISCLAIMER OF WARRANTY CONSTITUTES AN ESSENTIAL PART OF THIS LICENSE.
  * NO USE OF THE JASPER SOFTWARE IS AUTHORIZED HEREUNDER EXCEPT UNDER
  * THIS DISCLAIMER.  THE JASPER SOFTWARE IS PROVIDED BY THE LICENSORS AND
  * CONTRIBUTORS UNDER THIS LICENSE ON AN ``AS-IS'' BASIS, WITHOUT WARRANTY
  * OF ANY KIND, EITHER EXPRESSED OR IMPLIED, INCLUDING, WITHOUT LIMITATION,
  * WARRANTIES THAT THE JASPER SOFTWARE IS FREE OF DEFECTS, IS MERCHANTABLE,
  * IS FIT FOR A PARTICULAR PURPOSE OR IS NON-INFRINGING.  THOSE INTENDING
  * TO USE THE JASPER SOFTWARE OR MODIFICATIONS THEREOF FOR USE IN HARDWARE
  * OR SOFTWARE PRODUCTS ARE ADVISED THAT THEIR USE MAY INFRINGE EXISTING
  * PATENTS, COPYRIGHTS, TRADEMARKS, OR OTHER INTELLECTUAL PROPERTY RIGHTS.
  * THE ENTIRE RISK AS TO THE QUALITY AND PERFORMANCE OF THE JASPER SOFTWARE
  * IS WITH THE USER.  SHOULD ANY PART OF THE JASPER SOFTWARE PROVE DEFECTIVE
  * IN ANY RESPECT, THE USER (AND NOT THE INITIAL DEVELOPERS, THE UNIVERSITY
  * OF BRITISH COLUMBIA, IMAGE POWER, INC., MICHAEL DAVID ADAMS, OR ANY
  * OTHER CONTRIBUTOR) SHALL ASSUME THE COST OF ANY NECESSARY SERVICING,
  * REPAIR OR CORRECTION.  UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY,
  * WHETHER TORT (INCLUDING NEGLIGENCE), CONTRACT, OR OTHERWISE, SHALL THE
  * INITIAL DEVELOPER, THE UNIVERSITY OF BRITISH COLUMBIA, IMAGE POWER, INC.,
  * MICHAEL DAVID ADAMS, ANY OTHER CONTRIBUTOR, OR ANY DISTRIBUTOR OF THE
  * JASPER SOFTWARE, OR ANY SUPPLIER OF ANY OF SUCH PARTIES, BE LIABLE TO
  * THE USER OR ANY OTHER PERSON FOR ANY INDIRECT, SPECIAL, INCIDENTAL, OR
  * CONSEQUENTIAL DAMAGES OF ANY CHARACTER INCLUDING, WITHOUT LIMITATION,
  * DAMAGES FOR LOSS OF GOODWILL, WORK STOPPAGE, COMPUTER FAILURE OR
  * MALFUNCTION, OR ANY AND ALL OTHER COMMERCIAL DAMAGES OR LOSSES, EVEN IF
  * SUCH PARTY HAD BEEN INFORMED, OR OUGHT TO HAVE KNOWN, OF THE POSSIBILITY
  * OF SUCH DAMAGES.  THE JASPER SOFTWARE AND UNDERLYING TECHNOLOGY ARE NOT
  * FAULT-TOLERANT AND ARE NOT DESIGNED, MANUFACTURED OR INTENDED FOR USE OR
  * RESALE AS ON-LINE CONTROL EQUIPMENT IN HAZARDOUS ENVIRONMENTS REQUIRING
  * FAIL-SAFE PERFORMANCE, SUCH AS IN THE OPERATION OF NUCLEAR FACILITIES,
  * AIRCRAFT NAVIGATION OR COMMUNICATION SYSTEMS, AIR TRAFFIC CONTROL, DIRECT
  * LIFE SUPPORT MACHINES, OR WEAPONS SYSTEMS, IN WHICH THE FAILURE OF THE
  * JASPER SOFTWARE OR UNDERLYING TECHNOLOGY OR PRODUCT COULD LEAD DIRECTLY
  * TO DEATH, PERSONAL INJURY, OR SEVERE PHYSICAL OR ENVIRONMENTAL DAMAGE
  * ("HIGH RISK ACTIVITIES").  LICENSOR SPECIFICALLY DISCLAIMS ANY EXPRESS
  * OR IMPLIED WARRANTY OF FITNESS FOR HIGH RISK ACTIVITIES.  USER WILL NOT
  * KNOWINGLY USE, DISTRIBUTE OR RESELL THE JASPER SOFTWARE OR UNDERLYING
  * TECHNOLOGY OR PRODUCTS FOR HIGH RISK ACTIVITIES AND WILL ENSURE THAT ITS
  * CUSTOMERS AND END-USERS OF ITS PRODUCTS ARE PROVIDED WITH A COPY OF THE
  * NOTICE SPECIFIED IN THIS SECTION.
  *
  * __END_OF_JASPER_LICENSE__
  */

 /*
  * $Id$
  */

 /******************************************************************************\
 * Includes
 \******************************************************************************/

 #include <assert.h>
 #include <stdio.h>
 #include <string.h>
 #include <math.h>
 #include <stdlib.h>
 #include <stdarg.h>

 #include "jasper/jas_math.h"
 #include "jasper/jas_malloc.h"

 #include "jpc_fix.h"
 #include "jpc_cs.h"
 #include "jpc_flt.h"
 #include "jpc_util.h"

 /******************************************************************************\
 * Miscellaneous Functions
 \******************************************************************************/

 int jpc_atoaf(char *s, int *numvalues, double **values)
 {
 	static char delim[] = ", \t\n";
 	char buf[4096];
 	int n;
 	double *vs;
 	char *cp;

 	strncpy(buf, s, sizeof(buf));
 	buf[sizeof(buf) - 1] = '\0';
 	n = 0;
 	if ((cp = strtok(buf, delim))) {
 		++n;
 		while ((cp = strtok(0, delim))) {
 			if (cp != '\0') {
 				++n;
 			}
 		}
 	}

 	if (n) {
 		if (!(vs = jas_malloc(n * sizeof(double)))) {
 			return -1;
 		}

 		strncpy(buf, s, sizeof(buf));
 		buf[sizeof(buf) - 1] = '\0';
 		n = 0;
 		if ((cp = strtok(buf, delim))) {
 			vs[n] = atof(cp);
 			++n;
 			while ((cp = strtok(0, delim))) {
 				if (cp != '\0') {
 					vs[n] = atof(cp);
 					++n;
 				}
 			}
 		}
 	} else {
 		vs = 0;
 	}

 	*numvalues = n;
 	*values = vs;

 	return 0;
 }

 jas_seq_t *jpc_seq_upsample(jas_seq_t *x, int m)
 {
 	jas_seq_t *z;
 	int i;

 	if (!(z = jas_seq_create(jas_seq_start(x) * m, (jas_seq_end(x) - 1) * m + 1)))
 		return 0;
 	for (i = jas_seq_start(z); i < jas_seq_end(z); i++) {
 		*jas_seq_getref(z, i) = (!JAS_MOD(i, m)) ? jas_seq_get(x, i / m) :
 		  jpc_inttofix(0);
 	}

 	return z;
 }

 jpc_fix_t jpc_seq_norm(jas_seq_t *x)
 {
 	jpc_fix_t s;
 	int i;

 	s = jpc_inttofix(0);
 	for (i = jas_seq_start(x); i < jas_seq_end(x); i++) {
 		s = jpc_fix_add(s, jpc_fix_mul(jas_seq_get(x, i), jas_seq_get(x, i)));
 	}

 	return jpc_dbltofix(sqrt(jpc_fixtodbl(s)));
 }

 jas_seq_t *jpc_seq_conv(jas_seq_t *x, jas_seq_t *y)
 {
 	int i;
 	int j;
 	int k;
 	jas_seq_t *z;
 	jpc_fix_t s;
 	jpc_fix_t v;

 	z = jas_seq_create(jas_seq_start(x) + jas_seq_start(y),
 	  jas_seq_end(x) + jas_seq_end(y) - 1);
 	assert(z);
 	for (i = jas_seq_start(z); i < jas_seq_end(z); i++) {
 		s = jpc_inttofix(0);
 		for (j = jas_seq_start(y); j < jas_seq_end(y); j++) {
 			k = i - j;
 			if (k < jas_seq_start(x) || k >= jas_seq_end(x)) {
 				v = JPC_FIX_ZERO;
 			} else {
 				v = jas_seq_get(x, k);
 			}
 			s = jpc_fix_add(s, jpc_fix_mul(jas_seq_get(y, j), v));
 		}
 		*jas_seq_getref(z, i) = s;
 	}

 	return z;
 }
	/*
	* Copyright (c) 1999-2000 Image Power, Inc. and the University of
	* British Columbia.
	* Copyright (c) 2001-2003 Michael David Adams.
	* All rights reserved.
	*/

	/* __START_OF_JASPER_LICENSE__
	*
	* JasPer Software License
	*
	* IMAGE POWER JPEG-2000 PUBLIC LICENSE
	* ************************************
	*
	* GRANT:
	*
	* Permission is hereby granted, free of charge, to any person (the "User")
	* obtaining a copy of this software and associated documentation, to deal
	* in the JasPer Software without restriction, including without limitation
	* the right to use, copy, modify, merge, publish, distribute, sublicense,
	* and/or sell copies of the JasPer Software (in source and binary forms),
	* and to permit persons to whom the JasPer Software is furnished to do so,
	* provided further that the License Conditions below are met.
	*
	* License Conditions
	* ******************
	*
	* A. Redistributions of source code must retain the above copyright notice,
	* and this list of conditions, and the following disclaimer.
	*
	* B. Redistributions in binary form must reproduce the above copyright
	* notice, and this list of conditions, and the following disclaimer in
	* the documentation and/or other materials provided with the distribution.
	*
	* C. Neither the name of Image Power, Inc. nor any other contributor
	* (including, but not limited to, the University of British Columbia and
	* Michael David Adams) may be used to endorse or promote products derived
	* from this software without specific prior written permission.
	*
	* D. User agrees that it shall not commence any action against Image Power,
	* Inc., the University of British Columbia, Michael David Adams, or any
	* other contributors (collectively "Licensors") for infringement of any
	* intellectual property rights ("IPR") held by the User in respect of any
	* technology that User owns or has a right to license or sublicense and
	* which is an element required in order to claim compliance with ISO/IEC
	* 15444-1 (i.e., JPEG-2000 Part 1). "IPR" means all intellectual property
	* rights worldwide arising under statutory or common law, and whether
	* or not perfected, including, without limitation, all (i) patents and
	* patent applications owned or licensable by User; (ii) rights associated
	* with works of authorship including copyrights, copyright applications,
	* copyright registrations, mask work rights, mask work applications,
	* mask work registrations; (iii) rights relating to the protection of
	* trade secrets and confidential information; (iv) any right analogous
	* to those set forth in subsections (i), (ii), or (iii) and any other
	* proprietary rights relating to intangible property (other than trademark,
	* trade dress, or service mark rights); and (v) divisions, continuations,
	* renewals, reissues and extensions of the foregoing (as and to the extent
	* applicable) now existing, hereafter filed, issued or acquired.
	*
	* E. If User commences an infringement action against any Licensor(s) then
	* such Licensor(s) shall have the right to terminate User's license and
	* all sublicenses that have been granted hereunder by User to other parties.
	*
	* F. This software is for use only in hardware or software products that
	* are compliant with ISO/IEC 15444-1 (i.e., JPEG-2000 Part 1). No license
	* or right to this Software is granted for products that do not comply
	* with ISO/IEC 15444-1. The JPEG-2000 Part 1 standard can be purchased
	* from the ISO.
	*
	* THIS DISCLAIMER OF WARRANTY CONSTITUTES AN ESSENTIAL PART OF THIS LICENSE.
	* NO USE OF THE JASPER SOFTWARE IS AUTHORIZED HEREUNDER EXCEPT UNDER
	* THIS DISCLAIMER. THE JASPER SOFTWARE IS PROVIDED BY THE LICENSORS AND
	* CONTRIBUTORS UNDER THIS LICENSE ON AN ``AS-IS'' BASIS, WITHOUT WARRANTY
	* OF ANY KIND, EITHER EXPRESSED OR IMPLIED, INCLUDING, WITHOUT LIMITATION,
	* WARRANTIES THAT THE JASPER SOFTWARE IS FREE OF DEFECTS, IS MERCHANTABLE,
	* IS FIT FOR A PARTICULAR PURPOSE OR IS NON-INFRINGING. THOSE INTENDING
	* TO USE THE JASPER SOFTWARE OR MODIFICATIONS THEREOF FOR USE IN HARDWARE
	* OR SOFTWARE PRODUCTS ARE ADVISED THAT THEIR USE MAY INFRINGE EXISTING
	* PATENTS, COPYRIGHTS, TRADEMARKS, OR OTHER INTELLECTUAL PROPERTY RIGHTS.
	* THE ENTIRE RISK AS TO THE QUALITY AND PERFORMANCE OF THE JASPER SOFTWARE
	* IS WITH THE USER. SHOULD ANY PART OF THE JASPER SOFTWARE PROVE DEFECTIVE
	* IN ANY RESPECT, THE USER (AND NOT THE INITIAL DEVELOPERS, THE UNIVERSITY
	* OF BRITISH COLUMBIA, IMAGE POWER, INC., MICHAEL DAVID ADAMS, OR ANY
	* OTHER CONTRIBUTOR) SHALL ASSUME THE COST OF ANY NECESSARY SERVICING,
	* REPAIR OR CORRECTION. UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY,
	* WHETHER TORT (INCLUDING NEGLIGENCE), CONTRACT, OR OTHERWISE, SHALL THE
	* INITIAL DEVELOPER, THE UNIVERSITY OF BRITISH COLUMBIA, IMAGE POWER, INC.,
	* MICHAEL DAVID ADAMS, ANY OTHER CONTRIBUTOR, OR ANY DISTRIBUTOR OF THE
	* JASPER SOFTWARE, OR ANY SUPPLIER OF ANY OF SUCH PARTIES, BE LIABLE TO
	* THE USER OR ANY OTHER PERSON FOR ANY INDIRECT, SPECIAL, INCIDENTAL, OR
	* CONSEQUENTIAL DAMAGES OF ANY CHARACTER INCLUDING, WITHOUT LIMITATION,
	* DAMAGES FOR LOSS OF GOODWILL, WORK STOPPAGE, COMPUTER FAILURE OR
	* MALFUNCTION, OR ANY AND ALL OTHER COMMERCIAL DAMAGES OR LOSSES, EVEN IF
	* SUCH PARTY HAD BEEN INFORMED, OR OUGHT TO HAVE KNOWN, OF THE POSSIBILITY
	* OF SUCH DAMAGES. THE JASPER SOFTWARE AND UNDERLYING TECHNOLOGY ARE NOT
	* FAULT-TOLERANT AND ARE NOT DESIGNED, MANUFACTURED OR INTENDED FOR USE OR
	* RESALE AS ON-LINE CONTROL EQUIPMENT IN HAZARDOUS ENVIRONMENTS REQUIRING
	* FAIL-SAFE PERFORMANCE, SUCH AS IN THE OPERATION OF NUCLEAR FACILITIES,
	* AIRCRAFT NAVIGATION OR COMMUNICATION SYSTEMS, AIR TRAFFIC CONTROL, DIRECT
	* LIFE SUPPORT MACHINES, OR WEAPONS SYSTEMS, IN WHICH THE FAILURE OF THE
	* JASPER SOFTWARE OR UNDERLYING TECHNOLOGY OR PRODUCT COULD LEAD DIRECTLY
	* TO DEATH, PERSONAL INJURY, OR SEVERE PHYSICAL OR ENVIRONMENTAL DAMAGE
	* ("HIGH RISK ACTIVITIES"). LICENSOR SPECIFICALLY DISCLAIMS ANY EXPRESS
	* OR IMPLIED WARRANTY OF FITNESS FOR HIGH RISK ACTIVITIES. USER WILL NOT
	* KNOWINGLY USE, DISTRIBUTE OR RESELL THE JASPER SOFTWARE OR UNDERLYING
	* TECHNOLOGY OR PRODUCTS FOR HIGH RISK ACTIVITIES AND WILL ENSURE THAT ITS
	* CUSTOMERS AND END-USERS OF ITS PRODUCTS ARE PROVIDED WITH A COPY OF THE
	* NOTICE SPECIFIED IN THIS SECTION.
	*
	* __END_OF_JASPER_LICENSE__
	*/

	/*
	* $Id$
	*/

	/******************************************************************************\
	* Includes
	\******************************************************************************/

	#include <assert.h>
	#include <stdio.h>
	#include <string.h>
	#include <math.h>
	#include <stdlib.h>
	#include <stdarg.h>

	#include "jasper/jas_math.h"
	#include "jasper/jas_malloc.h"

	#include "jpc_fix.h"
	#include "jpc_cs.h"
	#include "jpc_flt.h"
	#include "jpc_util.h"

	/******************************************************************************\
	* Miscellaneous Functions
	\******************************************************************************/

	int jpc_atoaf(char s, int numvalues, double **values)
	{
	static char delim[] = ", \t\n";
	char buf[4096];
	int n;
	double *vs;
	char *cp;

	strncpy(buf, s, sizeof(buf));
	buf[sizeof(buf) - 1] = '\0';
	n = 0;
	if ((cp = strtok(buf, delim))) {
	++n;
	while ((cp = strtok(0, delim))) {
	if (cp != '\0') {
	++n;
	}
	}
	}

	if (n) {
	if (!(vs = jas_malloc(n * sizeof(double)))) {
	return -1;
	}

	strncpy(buf, s, sizeof(buf));
	buf[sizeof(buf) - 1] = '\0';
	n = 0;
	if ((cp = strtok(buf, delim))) {
	vs[n] = atof(cp);
	++n;
	while ((cp = strtok(0, delim))) {
	if (cp != '\0') {
	vs[n] = atof(cp);
	++n;
	}
	}
	}
	} else {
	vs = 0;
	}

	*numvalues = n;
	*values = vs;

	return 0;
	}

	jas_seq_t jpc_seq_upsample(jas_seq_t x, int m)
	{
	jas_seq_t *z;
	int i;

	if (!(z = jas_seq_create(jas_seq_start(x) * m, (jas_seq_end(x) - 1) * m + 1)))
	return 0;
	for (i = jas_seq_start(z); i < jas_seq_end(z); i++) {
	*jas_seq_getref(z, i) = (!JAS_MOD(i, m)) ? jas_seq_get(x, i / m) :
	jpc_inttofix(0);
	}

	return z;
	}

	jpc_fix_t jpc_seq_norm(jas_seq_t *x)
	{
	jpc_fix_t s;
	int i;

	s = jpc_inttofix(0);
	for (i = jas_seq_start(x); i < jas_seq_end(x); i++) {
	s = jpc_fix_add(s, jpc_fix_mul(jas_seq_get(x, i), jas_seq_get(x, i)));
	}

	return jpc_dbltofix(sqrt(jpc_fixtodbl(s)));
	}

	jas_seq_t jpc_seq_conv(jas_seq_t x, jas_seq_t *y)
	{
	int i;
	int j;
	int k;
	jas_seq_t *z;
	jpc_fix_t s;
	jpc_fix_t v;

	z = jas_seq_create(jas_seq_start(x) + jas_seq_start(y),
	jas_seq_end(x) + jas_seq_end(y) - 1);
	assert(z);
	for (i = jas_seq_start(z); i < jas_seq_end(z); i++) {
	s = jpc_inttofix(0);
	for (j = jas_seq_start(y); j < jas_seq_end(y); j++) {
	k = i - j;
	if (k < jas_seq_start(x) \|\| k >= jas_seq_end(x)) {
	v = JPC_FIX_ZERO;
	} else {
	v = jas_seq_get(x, k);
	}
	s = jpc_fix_add(s, jpc_fix_mul(jas_seq_get(y, j), v));
	}
	*jas_seq_getref(z, i) = s;
	}

	return z;
	}