1

/*

     2

 * (c) 2008 Edouard TISSERANT

     3

 *

     4

 * Offered to the public under the terms of the GNU Lesser General Public

     5

 * License as published by the Free Software Foundation; either version 2

     6

 * of the License, or (at your option) any later version.

     7

 *

     8

 * This program is distributed in the hope that it will be useful, but

     9

 * WITHOUT ANY WARRANTY; without even the implied warranty of

    10

 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser

    11

 * General Public License for more details.

    12

 *

    13

 * This code is made available on the understanding that it will not be

    14

 * used in safety-critical situations without a full and competent review.

    15

 */

    16

    17

/****

    18

 * IEC 61131-3 standard function library

    19

 */

    20

    21

#include <limits.h>

    22

#include <float.h>

    23

#include <math.h>

    24

#include <time.h>

    25

#include <stdint.h>

    26

    27

#include <stdio.h>

    28

#include <stdlib.h>

    29

#include <string.h>

    30

#include <stdarg.h>

    31

    32

#ifdef DEBUG_IEC

    33

#define DBG(...) printf(__VA_ARGS__);

    34

#define DBG_TYPE(TYPENAME, name) __print_##TYPENAME(name);

    35

#else

    36

#define DBG(...)

    37

#define DBG_TYPE(TYPENAME, name)

    38

#endif

    39

    40

/*

    41

 * Include type defs.

    42

 */

    43

#include "iec_types_all.h"

    44

    45

extern TIME __CURRENT_TIME;

    46

extern BOOL __DEBUG;

    47

    48

/* TODO

    49

typedef struct {

    50

    __strlen_t len;

    51

    u_int16_t body[STR_MAX_LEN];

    52

} WSTRING;

    53

*/

    54

    55

# if __WORDSIZE == 64

    56

#define __32b_sufix

    57

#define __64b_sufix L

    58

#else

    59

#define __32b_sufix L

    60

#define __64b_sufix LL 

    61

#endif 

    62

    63

#define __lit(type,value,sfx) (type)value##sfx

    64

// Keep this macro expention step to let sfx change into L or LL 

    65

#define __literal(type,value,sfx) __lit(type,value,sfx)

    66

    67

#define __BOOL_LITERAL(value) __literal(BOOL,value,)

    68

#define __SINT_LITERAL(value) __literal(SINT,value,)

    69

#define __INT_LITERAL(value) __literal(INT,value,)

    70

#define __DINT_LITERAL(value) __literal(DINT,value,__32b_sufix)

    71

#define __LINT_LITERAL(value) __literal(LINT,value,__64b_sufix)

    72

#define __USINT_LITERAL(value) __literal(USINT,value,)

    73

#define __UINT_LITERAL(value) __literal(UINT,value,)

    74

#define __UDINT_LITERAL(value) __literal(UDINT,value,__32b_sufix)

    75

#define __ULINT_LITERAL(value) __literal(ULINT,value,__64b_sufix)

    76

#define __REAL_LITERAL(value) __literal(REAL,value,__32b_sufix)

    77

#define __LREAL_LITERAL(value) __literal(LREAL,value,__64b_sufix)

    78

#define __TIME_LITERAL(value) __literal(TIME,value,)

    79

#define __DATE_LITERAL(value) __literal(DATE,value,)

    80

#define __TOD_LITERAL(value) __literal(TOD,value,)

    81

#define __DT_LITERAL(value) __literal(DT,value,)

    82

#define __STRING_LITERAL(value) __literal(STRING,value,)

    83

#define __BYTE_LITERAL(value) __literal(BYTE,value,)

    84

#define __WORD_LITERAL(value) __literal(WORD,value,)

    85

#define __DWORD_LITERAL(value) __literal(DWORD,value,__32b_sufix)

    86

#define __LWORD_LITERAL(value) __literal(LWORD,value,__32b_sufix)

    87

    88

    89

    90

#define __INIT_REAL 0

    91

#define __INIT_LREAL 0

    92

#define __INIT_SINT 0

    93

#define __INIT_INT 0

    94

#define __INIT_DINT 0

    95

#define __INIT_LINT 0

    96

#define __INIT_USINT 0

    97

#define __INIT_UINT 0

    98

#define __INIT_UDINT 0

    99

#define __INIT_ULINT 0

   100

#define __INIT_TIME (TIME){0,0}

   101

#define __INIT_BOOL 0

   102

#define __INIT_BYTE 0

   103

#define __INIT_WORD 0

   104

#define __INIT_DWORD 0

   105

#define __INIT_LWORD 0

   106

#define __INIT_STRING (STRING){0,""}

   107

//#define __INIT_WSTRING

   108

#define __INIT_DATE (DATE){0,0}

   109

#define __INIT_TOD (TOD){0,0}

   110

#define __INIT_DT (DT){0,0}

   111

   112

typedef union __IL_DEFVAR_T {

   113

    BOOL    BOOLvar;

   114

   115

    SINT    SINTvar;

   116

    INT     INTvar;

   117

    DINT    DINTvar;

   118

    LINT    LINTvar;

   119

   120

    USINT   USINTvar;

   121

    UINT    UINTvar;

   122

    UDINT   UDINTvar;

   123

    ULINT   ULINTvar;

   124

   125

    BYTE    BYTEvar;

   126

    WORD    WORDvar;

   127

    DWORD   DWORDvar;

   128

    LWORD   LWORDvar;

   129

   130

    REAL    REALvar;

   131

    LREAL   LREALvar;

   132

   133

    TIME    TIMEvar;

   134

    TOD TODvar;

   135

    DT  DTvar;

   136

    DATE    DATEvar;

   137

} __IL_DEFVAR_T;

   138

   139

/*****************/

   140

/* Misc internal */

   141

/*****************/

   142

   143

/* function that generates an IEC runtime error */

   144

static inline void IEC_error(void) {

   145

  /* TODO... */

   146

  fprintf(stderr, "IEC 61131-3 runtime error.\n");

   147

  /*exit(1);*/

   148

}

   149

   150

   151

static inline void __normalize_timespec (struct timespec *ts) {

   152

  if( ts->tv_nsec < -1000000000 || (( ts->tv_sec > 0 ) && ( ts->tv_nsec < 0 ))){

   153

    ts->tv_sec--;

   154

    ts->tv_nsec += 1000000000;

   155

  }

   156

  if( ts->tv_nsec > 1000000000 || (( ts->tv_sec < 0 ) && ( ts->tv_nsec > 0 ))){

   157

    ts->tv_sec++;

   158

    ts->tv_nsec -= 1000000000;

   159

  }

   160

}

   161

   162

static inline struct timespec __time_to_timespec(int sign, double mseconds, double seconds, double minutes, double hours, double days) {

   163

  struct timespec ts;

   164

   165

  /* sign is 1 for positive values, -1 for negative time... */

   166

  long double total_sec = ((days*24 + hours)*60 + minutes)*60 + seconds + mseconds/1e3;

   167

  if (sign >= 0) sign = 1; else sign = -1;

   168

  ts.tv_sec = sign * (long int)total_sec;

   169

  ts.tv_nsec = sign * (long int)((total_sec - ts.tv_sec)*1e9);

   170

   171

  return ts;

   172

}

   173

   174

   175

static inline struct timespec __tod_to_timespec(double seconds, double minutes, double hours) {

   176

  struct timespec ts;

   177

   178

  long double total_sec = (hours*60 + minutes)*60 + seconds;

   179

  ts.tv_sec = (long int)total_sec;

   180

  ts.tv_nsec = (long int)((total_sec - ts.tv_sec)*1e9);

   181

   182

  return ts;

   183

}

   184

   185

static inline struct timespec __date_to_timespec(int day, int month, int year) {

   186

  struct timespec ts;

   187

  struct tm broken_down_time;

   188

   189

  broken_down_time.tm_sec = 0;

   190

  broken_down_time.tm_min = 0;

   191

  broken_down_time.tm_hour = 0;

   192

  broken_down_time.tm_mday = day;  /* day of month, from 1 to 31 */

   193

  broken_down_time.tm_mon = month - 1;   /* month since January, in the range 0 to 11 */

   194

  broken_down_time.tm_year = year - 1900;  /* number of years since 1900 */

   195

   196

  time_t epoch_seconds = mktime(&broken_down_time); /* determine number of seconds since the epoch, i.e. Jan 1st 1970 */

   197

   198

  if ((time_t)(-1) == epoch_seconds)

   199

    IEC_error();

   200

   201

  ts.tv_sec = epoch_seconds;

   202

  ts.tv_nsec = 0;

   203

   204

  return ts;

   205

}

   206

   207

static inline struct timespec __dt_to_timespec(double seconds,  double minutes, double hours, int day, int month, int year) {

   208

  struct timespec ts;

   209

   210

  long double total_sec = (hours*60 + minutes)*60 + seconds;

   211

  ts.tv_sec = (long int)total_sec;

   212

  ts.tv_nsec = (long int)((total_sec - ts.tv_sec)*1e9);

   213

   214

  struct tm broken_down_time;

   215

  broken_down_time.tm_sec = 0;

   216

  broken_down_time.tm_min = 0;

   217

  broken_down_time.tm_hour = 0;

   218

  broken_down_time.tm_mday = day;  /* day of month, from 1 to 31 */

   219

  broken_down_time.tm_mon = month - 1;   /* month since January, in the range 0 to 11 */

   220

  broken_down_time.tm_year = year - 1900;  /* number of years since 1900 */

   221

   222

  time_t epoch_seconds = mktime(&broken_down_time); /* determine number of seconds since the epoch, i.e. Jan 1st 1970 */

   223

  if ((time_t)(-1) == epoch_seconds)

   224

    IEC_error();

   225

   226

  ts.tv_sec += epoch_seconds;

   227

  if (ts.tv_sec < epoch_seconds)

   228

    /* since the TOD is always positive, if the above happens then we had an overflow */

   229

    IEC_error();

   230

   231

  return ts;

   232

}

   233

   234

/********************/

   235

/*   EN/ENO PARAMS  */

   236

/********************/

   237

   238

#define EN_ENO_PARAMS BOOL EN, BOOL *ENO

   239

#define TEST_EN(TYPENAME)\

   240

  if (!EN) {\

   241

    if (ENO != NULL)\

   242

      *ENO = __BOOL_LITERAL(FALSE);\

   243

    return __INIT_##TYPENAME;\

   244

  }\

   245

  else if (ENO != NULL)\

   246

    *ENO = __BOOL_LITERAL(TRUE);

   247

#define TEST_EN_COND(TYPENAME, COND)\

   248

  if (!EN || (COND)) {\

   249

    if (ENO != NULL)\

   250

      *ENO = __BOOL_LITERAL(FALSE);\

   251

    return __INIT_##TYPENAME;\

   252

  }\

   253

  else if (ENO != NULL)\

   254

    *ENO = __BOOL_LITERAL(TRUE);

   255

   256

/***************/

   257

/*   Time ops  */

   258

/***************/

   259

#define __TIME_CMP(t1, t2) (t2.tv_sec == t1.tv_sec ? t1.tv_nsec - t2.tv_nsec : t1.tv_sec - t2.tv_sec) 

   260

   261

static inline TIME __TIME_ADD(TIME IN1, TIME IN2){

   262

  TIME res ={IN1.tv_sec + IN2.tv_sec,

   263

             IN1.tv_nsec + IN2.tv_nsec };

   264

  __normalize_timespec(&res);

   265

  return res;

   266

}

   267

static inline TIME __TIME_SUB(TIME IN1, TIME IN2){

   268

  TIME res ={IN1.tv_sec - IN2.tv_sec,

   269

             IN1.tv_nsec - IN2.tv_nsec };

   270

  __normalize_timespec(&res);

   271

  return res;

   272

}

   273

static inline TIME __TIME_MUL(TIME IN1, LREAL IN2){

   274

  LREAL s_f = IN1.tv_sec * IN2;

   275

  time_t s = s_f;

   276

  div_t ns = div((LREAL)IN1.tv_nsec * IN2, 1000000000);

   277

  TIME res = {s + ns.quot,

   278

              ns.rem + (s_f - s) * 1000000000 };

   279

  __normalize_timespec(&res);

   280

  return res;

   281

}

   282

static inline TIME __TIME_DIV(TIME IN1, LREAL IN2){

   283

  LREAL s_f = IN1.tv_sec / IN2;

   284

  time_t s = s_f;

   285

  TIME res = {s,

   286

              IN1.tv_nsec / IN2 + (s_f - s) * 1000000000 };

   287

  __normalize_timespec(&res);

   288

  return res;

   289

}

   290

   291

static inline TIME __date_and_time_to_time_of_day(EN_ENO_PARAMS, TIME IN){

   292

  TEST_EN(TIME)

   293

  return (TIME){IN.tv_sec % 86400, IN.tv_nsec};

   294

}

   295

static inline TIME __date_and_time_to_date(EN_ENO_PARAMS, TIME IN){

   296

  TEST_EN(TIME)

   297

  return (TIME){IN.tv_sec - (IN.tv_sec % (24*60*60)), 0};

   298

}

   299

static inline TIME __time_add(EN_ENO_PARAMS, TIME IN1, TIME IN2){

   300

  TEST_EN(TIME)

   301

  return __TIME_ADD(IN1, IN2);

   302

}

   303

static inline TIME __time_sub(EN_ENO_PARAMS, TIME IN1, TIME IN2){

   304

  TEST_EN(TIME)

   305

  return __TIME_SUB(IN1, IN2);

   306

}

   307

static inline TIME __time_mul(EN_ENO_PARAMS, TIME IN1, LREAL IN2){

   308

  TEST_EN(TIME)

   309

  return __TIME_MUL(IN1, IN2);

   310

}

   311

static inline TIME __time_div(EN_ENO_PARAMS, TIME IN1, LREAL IN2){

   312

  TEST_EN(TIME)

   313

  return __TIME_DIV(IN1, IN2);

   314

}

   315

   316

/***************/

   317

/* String ops  */

   318

/***************/

   319

#define __STR_CMP(str1, str2) memcmp((char*)&str1.body,(char*)&str2.body, str1.len < str2.len ? str1.len : str2.len)

   320

   322

    TEST_EN(UINT)

   323

    return IN.len;

   324

}

   326

    TEST_EN_COND(STRING, L < 0)

   327

    STRING res = __INIT_STRING;

   328

    L = L < IN.len ? L : IN.len;

   329

    memcpy(&res.body, &IN.body, L);

   330

    res.len = L;

   331

    return res;

   332

}

   334

    TEST_EN_COND(STRING, L < 0)

   335

    STRING res = __INIT_STRING;

   336

    L = L < IN.len ? L : IN.len;

   337

    memcpy(&res.body, &IN.body[IN.len - L], L);

   338

    res.len = L;

   339

    return res;

   340

}

   342

    TEST_EN_COND(STRING, L < 0 || P < 0)

   343

    STRING res = __INIT_STRING;

   344

    if(P <= IN.len){

   345

	    P -= 1; /* now can be used as [index]*/

   346

	    L = L + P <= IN.len ? L : IN.len - P;

   347

	    memcpy(&res.body, &IN.body[P] , L);

   348

        res.len = L;

   349

    }

   350

    return res;

   351

}

   353

  TEST_EN(STRING)

   354

  va_list ap;

   355

  UINT i;

   356

  __strlen_t charcount = 0;

   357

  STRING res = __INIT_STRING;

   358

   359

  va_start (ap, param_count);         /* Initialize the argument list.  */

   360

   361

  for (i = 0; i < param_count && charcount < STR_MAX_LEN; i++)

   362

  {

   363

    STRING tmp = va_arg(ap, STRING);

   364

    __strlen_t charrem = STR_MAX_LEN - charcount;

   365

    __strlen_t to_write = tmp.len > charrem ? charrem : tmp.len;

   366

    memcpy(&res.body[charcount], &tmp.body , to_write);

   367

    charcount += to_write;

   368

  }

   369

   370

  res.len = charcount;

   371

   372

  va_end (ap);                  /* Clean up.  */

   373

  return res;

   374

}

   376

    TEST_EN_COND(STRING, P < 0)

   377

    STRING res = __INIT_STRING;

   378

    __strlen_t to_copy;

   379

   380

    to_copy = P > IN1.len ? IN1.len : P - 1;

   381

    memcpy(&res.body, &IN1.body , to_copy);

   382

    P = res.len = to_copy;

   383

   384

    to_copy = IN2.len + res.len > STR_MAX_LEN ? STR_MAX_LEN - res.len : IN2.len;

   385

    memcpy(&res.body[res.len], &IN2.body , to_copy);

   386

    res.len += to_copy;

   387

   388

    to_copy = IN1.len - P < STR_MAX_LEN - res.len ? IN1.len - P : STR_MAX_LEN - res.len ;

   389

    memcpy(&res.body[res.len], &IN1.body[P] , to_copy);

   390

    res.len += to_copy;

   391

   392

    return res;

   393

}

   395

    TEST_EN_COND(STRING, L < 0 || P < 0)

   396

    STRING res = __INIT_STRING;

   397

    __strlen_t to_copy;

   398

   399

    to_copy = P > IN.len ? IN.len : P-1;

   400

    memcpy(&res.body, &IN.body , to_copy);

   401

    P = res.len = to_copy;

   402

   403

    if( IN.len > P + L ){

   404

        to_copy = IN.len - P - L;

   405

        memcpy(&res.body[res.len], &IN.body[P + L], to_copy);

   406

        res.len += to_copy;

   407

    }

   408

   409

    return res;

   410

}

   412

    TEST_EN_COND(STRING, L < 0 || P < 0)

   413

    STRING res = __INIT_STRING;

   414

    __strlen_t to_copy;

   415

   416

    to_copy = P > IN1.len ? IN1.len : P-1;

   417

    memcpy(&res.body, &IN1.body , to_copy);

   418

    P = res.len = to_copy;

   419

   420

    to_copy = IN2.len < L ? IN2.len : L;

   421

   422

    if( to_copy + res.len > STR_MAX_LEN ) 

   423

       to_copy = STR_MAX_LEN - res.len;   

   424

   425

    memcpy(&res.body[res.len], &IN2.body , to_copy);

   426

    res.len += to_copy;

   427

   428

    P += L;

   429

    if( res.len <  STR_MAX_LEN && P < IN1.len)

   430

    {

   431

        to_copy = IN1.len - P;

   432

        memcpy(&res.body[res.len], &IN1.body[P] , to_copy);

   433

        res.len += to_copy;

   434

    }

   435

   436

    return res;

   437

}

   438

   439

   440

   442

    UINT count1 = 0; /* offset of first matching char in IN1 */

   443

    UINT count2 = 0; /* count of matching char */

   444

    while(count1 + count2 < IN1->len && count2 < IN2->len)

   445

    {

   446

        if(IN1->body[count1 + count2] != IN2->body[count2++]){

   447

            count1 += count2;

   448

            count2 = 0;

   449

        }        

   450

    }

   451

    return count2 == IN2->len -1 ? 0 : count1 + 1;

   452

}

   454

    TEST_EN(UINT)

   455

    return __pfind(&IN1, &IN2);

   456

}

   457

   458

/***************/

   459

/* Convertions */

   460

/***************/

   461

    /***************/

   462

    /*  TO_STRING  */

   463

    /***************/

   464

static inline STRING __bool_to_string(EN_ENO_PARAMS, BOOL IN)

   465

{

   466

    TEST_EN(STRING)

   467

    if(IN)

   468

        return (STRING){4, "TRUE"};

   469

    return (STRING){5,"FALSE"};

   470

}

   471

static inline STRING __bit_to_string(EN_ENO_PARAMS, LWORD IN){

   472

    TEST_EN(STRING)

   473

    STRING res = __INIT_STRING;

   474

    res.len = snprintf((char*)res.body, STR_MAX_LEN, "16#%llx", IN);

   475

    if(res.len > STR_MAX_LEN) res.len = STR_MAX_LEN;

   476

    return res;

   477

}

   478

static inline STRING __real_to_string(EN_ENO_PARAMS, LREAL IN){

   479

    TEST_EN(STRING)

   480

    STRING res = __INIT_STRING;

   481

    res.len = snprintf((char*)res.body, STR_MAX_LEN, "%.10g", IN);

   482

    if(res.len > STR_MAX_LEN) res.len = STR_MAX_LEN;

   483

    return res;

   484

}

   485

static inline STRING __sint_to_string(EN_ENO_PARAMS, LINT IN){

   486

    TEST_EN(STRING)

   487

    STRING res = __INIT_STRING;

   488

    res.len = snprintf((char*)res.body, STR_MAX_LEN, "%lld", IN);

   489

    if(res.len > STR_MAX_LEN) res.len = STR_MAX_LEN;

   490

    return res;

   491

}

   492

static inline STRING __uint_to_string(EN_ENO_PARAMS, ULINT IN){

   493

    TEST_EN(STRING)

   494

    STRING res = __INIT_STRING;

   495

    res.len = snprintf((char*)res.body, STR_MAX_LEN, "%llu", IN);

   496

    if(res.len > STR_MAX_LEN) res.len = STR_MAX_LEN;

   497

    return res;

   498

}

   499

    /***************/

   500

    /* FROM_STRING */

   501

    /***************/

   502

static inline BOOL __string_to_bool(EN_ENO_PARAMS, STRING IN){

   503

    TEST_EN(BOOL)

   504

    return IN.len == 5 ? !memcmp(&IN.body,"TRUE", IN.len) : 0;

   505

}

   506

   507

static inline LINT __pstring_to_sint(STRING* IN){