IPP Software Navigation Tools IPP Links Communication Pan-STARRS Links

Changeset 5540


Ignore:
Timestamp:
Nov 17, 2005, 5:05:00 PM (21 years ago)
Author:
magnier
Message:

changed psImageSmooth from psF64 to psF32

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/eam_rel8_b0/psLib/src/imageops/psImageConvolve.c

    r5057 r5540  
    55 *  @author Robert DeSonia, MHPCC
    66 *
    7  *  @version $Revision: 1.26 $ $Name: not supported by cvs2svn $
    8  *  @date $Date: 2005-09-15 21:22:22 $
     7 *  @version $Revision: 1.26.4.1 $ $Name: not supported by cvs2svn $
     8 *  @date $Date: 2005-11-18 03:05:00 $
    99 *
    1010 *  Copyright 2004-2005 Maui High Performance Computing Center, University of Hawaii
     
    487487}
    488488
     489// XXX EAM : for now, image type must be F32
    489490void psImageSmooth (psImage *image,
    490491                    double  sigma,
     
    512513
    513514    // smooth in X direction
    514     temp = psVectorAlloc (Nx, PS_TYPE_F64);
     515    temp = psVectorAlloc (Nx, PS_TYPE_F32);
    515516    for (int j = 0; j < Ny; j++) {
    516         psF64 *vi = image->data.F64[j];
    517         psF64 *vo = temp->data.F64;
     517        psF32 *vi = image->data.F32[j];
     518        psF32 *vo = temp->data.F32;
    518519        for (int i = 0; i < Nx; i++) {
    519520            g = s = 0;
     
    528529            vo[i] = s / g;
    529530        }
    530         memcpy (image->data.F64[j], temp->data.F64, Nx*sizeof(psF64));
     531        memcpy (image->data.F32[j], temp->data.F32, Nx*sizeof(psF32));
    531532    }
    532533    psFree (temp);
    533534
    534535    // smooth in Y direction
    535     temp = psVectorAlloc (image->numRows, PS_TYPE_F64);
     536    temp = psVectorAlloc (image->numRows, PS_TYPE_F32);
    536537    for (int i = 0; i < Nx; i++) {
    537         psF64  *vo = temp->data.F64;
    538         psF64 **vi = image->data.F64;
     538        psF32  *vo = temp->data.F32;
     539        psF32 **vi = image->data.F32;
    539540        for (int j = 0; j < Ny; j++) {
    540541            g = s = 0;
Note: See TracChangeset for help on using the changeset viewer.