source: soft/giet_vm/applications/rosenfeld/nrc2/src/nrarith2x.c @ 810

/* ------------------ */
/* --- nrarith2.c --- */
/* ------------------ */

/*
 * Copyright (c) 2000-2014, Lionel Lacassagne, All rights reserved
 * Univ Paris Sud XI, CNRS
 * 
 * Distributed under the Boost Software License, Version 1.0
 * see accompanying file LICENSE.txt or copy it at
 * http://www.boost.org/LICENSE_1_0.txt
 */

#include <stdio.h>
#include <stddef.h>
#include <stdlib.h>
#include <malloc.h>
#include <math.h> // fabs
// #include <memory.h> // memcpy

#include "mypredef.h"
#include "nrtype.h"
#include "nrdef.h"
#include "nrmacro.h"
#include "nrkernel.h"

/* ------------------------------------------------------------------------------------- */
IMAGE_EXPORT(void) addcnz_bmatrix(byte **src,long nrl,long nrh,long ncl, long nch, byte  cte, byte **dst)
/* ------------------------------------------------------------------------------------- */
{
        long i,j;
        byte *Xi, *Yi;
        
        for(i=nrl; i<=nrh; i++) {
                Xi = src[i];
                Yi = dst[i];
                for(j=ncl; j<=nch; j++) {
                        if(Xi[j])
                                Yi[j] = Xi[j] + cte;
                        else
                                Yi[j] = Xi[j];
                }
        }
}
/* ----------------------------------------------------------------------------------- */
IMAGE_EXPORT(void) addandc_bmatrix(byte **src,long nrl,long nrh,long ncl, long nch, byte  cte, byte **dst)
/* ----------------------------------------------------------------------------------- */
{
        long i,j;
        byte *Xi, *Yi;
        
        for(i=nrl; i<=nrh; i++) {
                Xi = src[i];
                Yi = dst[i];
                for(j=ncl; j<=nch; j++) {
                        if(Xi[j])
                                Yi[j] = Xi[j] + cte;
                        else
                                Yi[j] = Xi[i];
                }
        }
}
/* ---------------------------------------------------------------------------------------- */
IMAGE_EXPORT(void) addandc_si16matrix(sint16 **src,long nrl,long nrh,long ncl, long nch, short cte, sint16 **dst)
/* ---------------------------------------------------------------------------------------- */
{
        long i,j;
        sint16 *Xi, *Yi;
        
        for(i=nrl; i<=nrh; i++) {
                Xi = src[i];
                Yi = dst[i];
                for(j=ncl; j<=nch; j++) {
                        if(Xi[j]) Yi[j] = Xi[j] + cte;
                }
        }
}
/* ------------------------------------------------------------------------------------------- */
IMAGE_EXPORT(void) addandc_ui16matrix(uint16 **src,long nrl,long nrh,long ncl, long nch, short cte, uint16 **dst)
/* ------------------------------------------------------------------------------------------ */
{
        long i,j;
        uint16 *Xi, *Yi;
        
        for(i=nrl; i<=nrh; i++) {
                Xi = src[i];
                Yi = dst[i];
                for(j=ncl; j<=nch; j++) {
                        if(Xi[j]) Yi[j] = Xi[j] + cte;
                }
        }
}
/* ----------------------------------------------------------- */
IMAGE_EXPORT(int) count_bmatrix(byte **m, long nrl,long nrh,long ncl, long nch)
/* ----------------------------------------------------------- */
{
        long i, j;
        int s = 0;
        byte *Xi;
        
        for(i=nrl; i<=nrh; i++) {
                Xi = m[i];
                for(j=ncl; j<=nch; j++) {
                        s += Xi[j];
                }
    }
    return s;
}